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1. Kapitel: Allgemeines 
§ 1. Einleitung i : x 


, a) Das Problem. Gegeben sei eine gewéhnliche Differentialgleichung erster _ 
Ordnung. Die Zeit ¢ sei die unabhangige, x die abhangige Variable. Die Glei- esa 


- chung kann in der Form a ; ie. 
j —— = «(x,t Lig 
We _ (x,t) | (1). 

1) Assistenz fiir Mathematik, ETH. 

' 2) Statistical Laboratory, The University. ; : ‘g 
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geschrieben werden. Wir betrachten nun den folgenden stochastischen Pro- 
zess x(t). s 

In jedem Zeitintervall (¢, ¢ + dt) besteht eine Wahrscheinlichkeit w(x, t) dt, 
dass x einen Sprung ausfiihrt. Wenn dieser zur Zeit ¢ stattfindet und wenn die 
abhangige Variable in (¢ — 0) den Wert y angenommen hat, so hat x unmittelbar — 
nachher die Verteilung G(y, x, ¢). Vom so erreichten Punkt x aus folgt der Pro- 
zess bis zum nachsten Sprung der durch (x, t) gehenden Lésung der Differential- 
gleichung. In Figur 1 sind einige Kurven der Lésungsschar und ein méglicher 


x 


Figur 1 
Moéglicher Verlauf eines Prozesses. 


Verlauf des Prozesses x(t) dargestellt. Der Anfang des Prozesses sei durch die 
Verteilung F(x, 0) gegeben. 
Wir kennen also 


die Differentialgleichung Kx Te 

die Dichte der Spriinge w(y,t) , ‘ 
die Verteilung nach diesen Gly, %,4) , 

die Anfangsverteilung £(%,,0) 


und wollen berechnen?) 
die Verteilung von x zu irgendeiner Zeit Els, by ; 
b) Methoden und Resultate. Die Lésungen werden auf zwei Arten erhalten. : 
Mit Integralgleichungen gewinnt man einen besondern Aspekt und allgemeinere — 
Ergebnisse des « Queueing »-Problems und ahnlicher Fragen (vgl. KENDALL[6])). 
Durch die Laplace-Transformation lasst sich unter anderem das «Gerdusch- 


. . 
8) Auf das Problem der Verteilung im Funktionenraum [im Lehrbuch von Doos [3] mit Q , 


bezeichnet: x,(@), w € Q] treten wir hier nicht ein. Auch die Autokorrelation des Prozesses 
wird nicht untersucht. 


4) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 303. 


i 


a 
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problem» behandeln. Hier wird die Stérung jeck oi als unstetig angenommen, 


wahrend in der Literatur das Hauptgewicht auf den kontinuierlichen Fall 
gelegt wird. 


§.2. Berspiele 


a) Wir betrachten ein Reservoir in einer wasserarmen Gegend (zum Beispiel 
Jura). Dieses werde direkt durch Regenwasser und nicht durch eine Quelle 
_ gespiesen. x(/) sei das Volumen des gespeicherten Wassers. Wir nehmen an, 
Regengiisse erfolgen mit einer Wahrscheinlichkeitsdichte w(t) und seien nur 
von kurzer Dauer. Wenn das Wasser sofort in das Reservoir geleitet wird, so 
nimmt x jeweils sprunghaft zu. —x(x, t) ist die pro Zeiteinheit gebrauchte 
Wassermenge. Bei leerem Reservoir ist x(0, t) = 0. Wenn mehr Wasser als die 
maximale Kapazitat A anfallt, so tiberlauft es sofort. Dies bedeutet x(x, 2) =—oo 
fiir x = A. Im Intervall 0 < x < A ist x negativ. 

Ahnlich verhalt es sich mit Stauseen fiir Kraftwerke. Die Bedingung, dass 
das Regenwasser sofort in den See gelangt, ist allerdings oft nicht erfiillt, 
besonders wenn Gletscher oder Walder im Einzugsgebiet liegen. Ein stetiger 
Zufluss kann jedoch beriicksichtigt werden, indem dieser zu x addiert wird. 
Fiir die folgenden Rechnungen muss aber x eine bestimmte Funktion von x 
und ¢ sein. Weder der allfallige stetige Teil des Zuflusses noch der fiir die Kraft- 
erzeugung bendtigte Abfluss diirfen ein stochastisches Element enthalten. 

b) Im Lagerungsproblem ist x(¢) die Menge der zur Zeit ¢ gelagerten Giiter. 
w ist die Wahrscheinlichkeitsdichte, dass eine neue Sendung eintrifft. G(y, x, ¢) 
ist die Verteilung von x nach deren Ankunft, wenn vorher die Warenmenge y 
betragen hat. Wir nehmen an, der Verkauf erfolge stetig. —% ist a Verkaufs- 

geschwindigkeit. Auch hier hat man die Bedingung 


4(0,4)=0, #(x,t)<0 fir «>0. 


c) Auch das Warteproblem (queueing problem) kann in dieser Form dar- 
gestellt werden. x(¢) ist jedoch nicht die Anzahl wartender Personen, sondern 
die durch diese fiir den Bedienenden aufgespeicherte Arbeit. w(y) ist die Wahr- 
scheinlichkeitsdichte, dass sich jemand anschliesst, wenn die Lange der 
Schlange, ausgedriickt in Arbeit, y betragt. G(x — y) ist die Verteilung der 
Arbeit, die ein Ankommender bringt. —x(x, | ist die Geschwindigkeit, mit 


- welcher der Bedienende arbeitet. ° 


d) Ein Geiger-Miller-Zahlrohr lost bei jedem registrierten Partikel einen 
kurzen Stromstoss aus. Dieser gelange auf einen C R-Kreis. Wie lautet die 
Wahrscheinlichkeitsverteilung der Spannung « am Kondensator? Bei jedem 
Stromstoss nimmt letztere sprunghaft zu und folgt zwischen den Unstetigkeiten 
der Differentialgleichung 


x 
eS Se 


Cc 
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Zahlrohr 


Figur 2 
Die Schaltung eines Zahlrohrs. 


Ahnlich sind die Verhaltnisse bei den St6rungen in einem Radio, wenn 
diese unstetig, und nicht, wie meist angenommen wird, kontinuierlich erfolgen. 


§ 3. Die Hauptformel 


Ableitungen seien im Sinn von SCHWARTZ [10] verstanden. Wir betrachten 
die Anderungen von F(x, ¢) langs einer Lésung der Differentialgleichung. Es ist © 
also dx = x dt. F(x + dx,t+ dt) setzt sich aus zwei Teilen zusammen: wenn — 


x 


Figur 3 


y < « ist und kein Sprung stattfindet, was mit der Wahrscheinlichkeit 
[B= doy.) aFly, 0 


erfolgt, bleibt F langs der Lésung der Differentialgleichung konstant. Dann 


ist noch der Einfluss eines Sprunges zu beriicksichtigen. Wir erhalten die 
Gleichung 


v 
— co 


E(x + dx, t+ dt) = fa — dt wy, t)] dF (y, t) + at | wly, t) G(y, x, t) dF(y, t) . 
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Nun fihren wir das totale Differential dF ein und dividieren durch dt: 


is?) 


sp % +5 = ~ | oly, t) dF ly, 2) Sid, io(y, t) G(y, x, t) dF (y, t) . 


Wir bezeichnen mit E(z) die Heavisidesche Funktion 
E(z)=1 fir <20, E(z)=0 fir z-<0, 
und definieren H(y, x, ¢) durch 
Hy, x,t) = E(x — y) — G(y, x, t). (2) 
Damit erhalten wir die fundamentale Formel 


oF 


ear - 
sp atsHr J w(y, t) H(y, x, t) dF(y, 2). (3) 


2. Kapitel: L6sung durch Integralgleichungen 


§ 4. Herleitung der Integralgleichung 


Diese Methode ist geeignet zur Lésung des stationaren Falles der Beispiele 
a), b) und c) von § 2. 
x sei eine nicht negative Variable. Wir setzen a(x) = — x und verlangen 


Oya ala Os fire 05 (4) 


Es sei f(x, t) = 0F/0x. Wenn wir die Abhangigkeit von ¢ nicht mehr explizit 
ausdriicken, erhalt (3) die Form 


f(2) a(x) — 5, = [ (0) Hy, ») dF iy) . 


Da x nicht negativ ist, geniigt es, die Integration von — 0 bis oo auszufiihren. 
Wir zerlegen den Bereich in (—0,+0) und (+0,00). Fiir alle positiven x 
diirfen wir durch a dividieren und erhalten nach einer kleinen Umstellung 


co 


Ay, (0) H(0, x) 1 OF 
f(x) — | POE” 19) dy = ay FO+ aay ae) 


+0 


Diese Gleichung hat die Form 
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mit 
Aly, (0) H(0, x) SS vem 
K(x, y) = HA JE SAU >, AX) = a(x) etek a(x) Ot ° 


a(*) 


Im stationdren Zustand fallt der letzte Term weg. Damit wird (6) eine 
Integralgleichung vom Typ von FrepHoim. Zur Lésung beniitzen wir die 
Neumannsche Reihe. Wir setzen 


K1(%, y) = K(x, 9), 
K(x, y) = ih K(x, u) K"-"1(u, y) du = fea, u) K(u, y) du, 


K*(x, yj K+ K? Ke 4 RA eee 


Ferner brauchen wir die Reziprozitatsformeln 
Sak + [ K(e, u) K*(u, y) du=K + [ K+, u) K(u, vy) du. (7) 
6 i) : 
Im stationaren Fall ist 
(x) = F(0) K(x, 0). (8) 


Somit erhalt man aus (6) 


Die Lésung davon ist 


(x) = F(0) K(x, 0) + Kee u) F(0) K(u, 0) du = F(0) K*(x, 0) . 
0 
Integriert man die Gleichung iiber x, so ergibt sich 
1 — F(0) = F(0) a K*(x, 0) dx, 
0 


F(0) = i ee CIB 


co , 


lee) 
1 + [K*(x, 0) dx 1+ [K*(z, 0) dx 
0 


0 


Dies ist die Lésung fiir den stationiren Fall. 
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Falls alle Spriinge positiv sind, was wir fiir dieses Kapitel im folgenden 
voraussetzen wollen, so ist G(y, x) = 0 fiir y > x. Diese Eigenschaft ubertragt 
sich auch auf A(y, x) und damit auf K(x, y). In diesem Fall ist das Integral 
in (6) nur bis x zu erstrecken; die Gleichung ist vom Typ von VOLTERRA. 


§ 5. Eigenschaften der Kerne 
und elementare Losungsmethoden fiir Volterra-Gleichungen 


a) K*(x, y) hangt von x, y und von der Funktion K(w, v) ab, jedoch nur von 
den Werten im Bereich y Sv Su XS x. 


Vv 


x 
Figur 4 
Der Bereich, in dem K(u, v) auf K*(x, y) einen Einfluss ausiibt. 


Beweis durch Induktion nach . Da es sich um Volterra-Gleichungen han- 
delt, ist K"(x, y) = 0 fiir y > x. Man kann daher den Integrationsbereich auf 
_ y <u S x beschranken und erhalt 


K(x, 9) = [ K(x, u) K(u, 9) du. 
J 


Nach der Induktionsannahme hangen nun die Faktoren des Integrals und 
die Grenzen nur von x, y und von K(w, v) im angegebenen Bereich ab, und 
somit gilt dies auch von der linken Seite. Fiir 1 = 1 ist die Behauptung erfiillt. 
- Durch die Summenbildung wird die Eigenschaft nicht zerstort und gilt daher 
auch fiir K*(x, y). | 

b) Wenn 

L(x, y) 2 K(x, y) 20, 
dann ist ) 
IEF (yyy eI (ay) 29. (10) 
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Denn bei der Bildung der iterierten und lésenden Kerne werden nur Pro- 
dukte, Integrale und Summen positiver Gréssen gebraucht. Somit bleiben die 
Ungleichungen erhalten. (Dieser Satz gilt allgemein fiir Fredholm-Gleichungen, 
wenn die Neumannsche Reihe konvergiert.) 

Es ist w > 0 und a > O. Da nach Voraussetzung alle Spriinge positiv sind, 
ist H > 0. Somit ist auch K(x, y) nicht negativ. 

c) K(x,y) = M(x, y) + L(x, 9). 

L sei eine kleine Korrektur, so dass die héhern Iterierten davon vernach- 
lassigt werden kénnen. Diejenigen von M sollen aber alle beriicksichtigt werden. 
Wir fiihren noch den Einheitskern I(x, y) ein durch die Definition 


I(x, y) = 6(x — y) (Diracsche 6-Funktion, vgl. Scuwartz [10)]). 


Wir schreiben hier die Faltung kurz als Produkt: 


[Ke u) L(u, y) du= KL. 


0 


Fir jeden Kern ist M= MI=IM. Ferner setzen wir M’=J + M* und 
behaupten 
K* = M*+M'LM’+M'LM'LM’+.---. (11) 


Wir verifizieren die Gleichung durch die Reziprozitatsformel. Es ist 
K+ KK*=M+L-+ (M+L) (M*4+ M’LM'+ M'LM'LM'+..-) 
=M+MM*+ MM'LM'+MMLM'LM'+... 
+L+2LM*+LM'LM’'+.... 
Nun ist nach (7) M+ M M* = M*, und aus der Definition von M’ folgt | 
L+LM*=LM' und 1+MM'=14+M+4+MM*=I1+M*t=M". 
Also wird 
K+ KK* = M*+M'LM’'+ M'LM'LM'+.--=K*, 


was zu zeigen war. (Auch dies gilt, falls die Reihe konvergiert, im Fredholm- 
schen Fall.) 
d) K sei ein Faltungskern: K(x, y) = K(x — y). Durch Induktion folgt, 


dass K” und K* auch nur von x — y abhangen. Wir fiihren die Laplace-Trans- 


formation ein: 
[o.e} 


M(s) = i e** K(z) dz. 


0 
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Die Transformierte von K” ist die n-te Potenz von M. So wird 


ra M(s) , 
M*(s) = M(s) + M%X(s) +--+ = T—- M@ fir |M(s)|<1. (12) 
Von M* erhalt man A*(x) durch Umkehrung der Transformation. 
Wenn K(x — y) nicht analytisch ist, so ist im allgemeinen auch K* nicht 
analytisch. In diesem Fall diirfte die Umkehrung der Laplace-Transformation 
schwierig sein. Hingegen erhalt man auch in diesem Fall leicht 


ie 2) 
- 


[ K*(x) dx = M*(0) = er gy wenn |M()|<1. (13) 
0 


Dadurch sind die Normierung von F und mit den Ableitungen von M*(s) 
in O auch die Momente gegeben. 

Diese Methode ist nicht identisch mit der von Kapitel 3. Hier wird die 
Laplace-Transformation zur Lésung der Integralgleichung (6), dort zur Lésung 
der Integro-Differentialgleichung (3) verwendet. 


A 


Figur 5 
Die Zerlegung des Gebietes 0 SS Srg 


e) K(x, y)= L(x, y)in «<A, 
K(x, y)=M(x,y)n y24, 
K(x, y) = N(x, y) in Meee Ni marl 


(vgl. Figur 5). Nach a) ist im Dreieck 0 Sy S ¥ <A K*=L* und im Bereich 


yon M K* = M*. Im Gebiet von N gehen wir von den Reziprozitatsformeln 


ee eh see | 
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aus. Es ist 


K*(x, y) = K(x, y) + Ne K(x, ) K*(u, y) du + ff K(x, u) K*(u, y) du 
Ad A 


A ta 
a= IN (A) + | N(x, u) L*(u, y) du + | M(x, u) K*(u, y) du. 
vi A 


Die ersten beiden Summanden sind bekannt. Nun halten wir y fest und be- 
nitzen die Abktirzungen 


A 
ful) = KC, 9), gy(2) = N(x, 9) + [NC u) E(u, 9) dee. 
y 
Damit erhalten wir die Integralgleichung 


fol) = gn(a) + [ M(x, ) fylu) de. 
Diese hat die Lésung 
ful) = K*(x, 9) = gy(a) + | MA (x, u) Gulu) de. (14) 


Damit ist das Problem gelést. 
Das Verfahren kann auch iteriert werden, indem die Definitionsbereiche 


von M oder L wieder unterteilt werden. 
Zwei Spezialfalle: 
al (a, 9) = 2h, y) ar eA 
=e (0) tir ex, 2 AS 


K*(x, y) = L*(x, y) fir x <A, 


= fire eA 

2 Ka, 9), = Le; 9). oan ae 
= N(%,.y). in 4 Aly <A, 
=0 Le BPs o78 4 


Hier ist M = 0. Also verschwindet auch M*, und wir haben 


A 
K*(x, y) = N(x, y) SNC: u) L¥*(u, y)\du in x =A, y< A 
aA 


f) Es sei 
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Im Zahler steht eine beliebige Funktion von x, im Nenner dieselbe Funktion 
von y. Mit Hilfe der Reziprozitatsformeln verifiziert man leicht 


K* f(x) *(y 
(x,y) = Hay We) - (15) 
g) Es sei 
K(x, y) = f(x) ely) fir x2y, K(x, y)=0 fir «< y. (16) 


GoursaT [5] hat auch die Verallgemeinerung auf 
= Sh) e(y) fir-x Sy, K(x,y)=0 fir x<y 


behandelt. Aus der Reziprozitatsformel erhalten wir 
K*(x, y) = f(x) fa (y) + fa Vs ya 


Nun leiten wir partiell nach x ab (auch hier im Sinn von SCHWARTZ): 


E(x, y) = f(x) Fay +f) ele) K*(e, 9) = K(x, ) (FS) + fe) ol), 
Se InKe(x, 9) = FE. + Ho) el) 


und analog 
Sy In K*(x, 9) = S15 — H9) elo) 


Das Minuszeichen kommt daher, dass das Integral hier nach der untern 
Grenze abgeleitet wird. Durch Integration erhalt man 


In K*(x, y) = In f(x) g(y) + ic i(u) (17) 


h) K(x, y) kénne im Dreieck 0 S$ y S$ x < A durch die Potenzreihe 
y= Om b;; xy? 
a] 
dargestellt werden. Fiir K*(x, y) setzen wir die Reihe 


= y anm * 
nm 
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an. Wir setzen diese in die Reziprozitatsformel 


x 
me iy cS 
y am x0 y= ST big xt yh + | D Dig xt DT ae y! du 
nm 4 y 1] kl 


yititk+ y? xt yl tirk +1 
= Shy w+ § hate eas jthk+1 


agkl 


ein. Durch Koeffizientenvergleich erhalt man 


1 
anm = One ene —i—j—1, m bi5 We ~ — Dy Ons A, m—5—Kk-1 7+ Re SS 2, +1° (18) 


1,7=0,1,2,...,44+7Sn—-1 4,k=0,1,2,....93+2 Sm-1 


Die Koeffizienten kénnen rekursiv berechnet werden. Insbesondere ist 


any = Ong = Png t 3) b:4 Oni-ina G— a4 GF 091, 2, 1 eee 


: (19) 


Da wir nur K*(x, 0) brauchen, fiihren wir den Konvergenzbeweis nur fiir 
> 4, x" durch. Fiir jedes B< A konvergiert die Reihe fiir K(x, y) in x SB 
absolut und gleichmassig. Also gilt in diesem Bereich 


L(x, y) =F [o3|#wsM. 
iy 


Wenn wir in (19) 6,; durch |6,;| ersetzen, erhalten wir statt der a, Koeffi- 
zienten, die wir mit c, bezeichnen. Es sind alle c, = 0. Der Konvergenzradius 
von 3} c, x” sei 9. In 0 <x< wird L*(x, 0) durch diese Reihe pete 
Wir miissen zeigen, dass 9 = B ist. 

Wir wahlen N so, dass 


See i 
b;;\"Bt Bi Se,) 6 
2 7 3B 
und setzen 
C, = Max (ec, B,¢, B* ..5,¢,-7 BD?) 
und 
b,= Max (|b, 0| Br", Lew Brr, aes | Yok s amet 5 


Es ist 


rhage) 


Cy BY = |Byg| BY + SY” |b,4| BY BBrt-i-1¢, 1 | 
aq 


n—1%~ 


SPr+ BC, 3" |b,,;| Be id ) 
aj 


a 


Vol. VII, 1956 Einige unstetige stochastische Prozesse 285 
Wir zerlegen die Summe ini +7 < N undi+j>N.Essein> N. 
i 
Cn Br = es -f- G. B ‘= oN M +. B Ce Sa 


Wir wahlen nun » so gross, dass 


nONS3BM und ARE 


Damit wird 
Get = 


Daraus folgt, dass die c,, B” beschrankt sind. Also konvergiert die Reihe in 
x < B. Da dies fiir jedes B < A gilt, ist der Konvergenzradius mindestens A. 


Nun ist aber |a,| <c,, man hat also eine konvergente Majorante der Reihe 
mx *(x, 0) = Dd) a, x". 


§6. Anwendungen 
a) Die Sprunghéhe sei exponentiell verteilt, das heisst 
Gly xh Shi 2 FOF often we Sy 
-(vgl. § 11). Es ist somit 
H(y,x)=e9-** fir x>y, H(y,x)=0 fir ~<y; a=a(x), w=v(y). 


a und w sind beliebig. Dies bedeutet im Beispiel des Schlangestehens, dass die 
Geschwindigkeit, mit welcher der Beamte arbeitet, und die Wahrscheinlich- 
keitsdichte, dass sich jemand anschliesst, beliebig von der Lange der Schlange 
abhangen kénnen. Es ist 


enka 


K(x, y) ev w(y) = f(x) e(y) fir x2y. 


a a(x) 


Somit erhalt man nach §5 g) 


In K*(x, y) = k(y — x) + In w(y) — In a(x) +f A du 
und damit insbesondere 
[ w(u)/a(u) du 
= 19)" ep a0 
K*(x, 0) = awe a , (20) 
b) Beschrankung der Variablen. Bei Reservoiren wird “% = —oo firx 2A, 


und somit K=0. Aus §5e) ist ersichtlich, dass K*(x, 0) =0 fir x 2 4, 
-wihrend fiir x <A K* nicht von A abhangt. Aus Gleichung (9) ist leicht der 

Einfluss auf F(x) abzulesen, der von einer Veranderung von A herrithrt: die 
‘Verteilung wird in A abgeschnitten (truncated) und neu normiert. 
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K = 0 fiir y => A kann im Schlangenproblem auftreten, wenn zum Beispiel 
ein Arzt keine neuen Anmeldungen entgegennimmt, falls er mit Arbeit tiber- 
lastet ist. Wenn also y = A ist, so wird w = 0. 

c) Ein von BLanc-Larrerre und Fortet [1], Seite 31, angegebenes Beispiel 
(Ausfluss von Wasser aus einem reat das unten eine Offnung hat und tropfen- 
weise aufgefiillt wird). Hier ist a(x —Vx. Wenn die Verteilung der Tropfen- 
grésse (also der Sprunghdhe) ea ist, so wird nach dem obigen, mit 
w= 1 

falvu du 1 : q 
K*(x, 0) = —— e-*e 60 oe ee (21) 
yx yx 

Das zur Normierung notwendige Integral kann durch die Substitution 
v =\Vx auf ein Fehlerintegral zuriickgefiihrt werden. 

Wenn die Sprunghdhe konstant = c ist, so wird K(x, y) = w( (y) Vx im 
Halbstreifen 0 < x — y<c, y =O, sonst 0. Mit einer ahnlichen Methode wie 
in §5e) kann K*(x, y) sukzessivein OS %¥—y<c, cS x—y<2e,.:. bee 
rechnet werden. Die Rechnungen werden jedoch bald umstandlich. 

Wenn die Sprunghohe klein ist im Verhaltnis zum mittleren Wasserstand, 
so kann % naherungsweise linearisiert werden. Falls ausserdem w in bezug auf 
y linearisiert werden darf und G nur von x — y und ¢ abhangt, so lasst sich die 
Lésung auch im nichtstationaéren Fall mit der Methode des § 8 berechnen. 

d) Ein Beispiel, wie die Methode von § 5 e) weitergefiihrt werden kann. Die 
Sprunghohe sei konstant = 1, a und w seien konstant. Wir setzen w/a = c. So 
wird, 0< y S x stillschweigend vorausgesetzt, 


K(x, y)=c¢ in OSx%—y< 1, K(x, 9)=0 In Lae gee 


lni07Ss 4S ist 
K*(x, 9) 0 ee ; 
lo lS 4 92 wird 


1 1 j 
g(x) = N(x, 0) + [ N(x, u) L*(u, 0) du = 0 + fe cetdu=ce—ce@-v, 
0 


Wegen K = K(x —y)- istin 1 Su Ss 4 <2 
M*(x, u) = K*(x — u) = c ee) , 
Setzt man dies in (14) ein, so erhalt man durch Integration 


: 
K*(x, 0) = e@-2 (c & + ¢2— ¢ — ¢2 x). 
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Dieses Resultat legt den Ansatz nahe 
Kas Ol P(ayet* ga 2h. 


Der Index von K* gibt an, dass die Formel fiir » < x<n-+ 1 gilt. P, ist 


ein Polynom u-ten Grades. Wir verifizieren die Gleichung durch die Reziprozi- 
tatsformel. Fiir 1 = 1 ist 


x-—1 n x 
K*(x, 0) = K(x, 0) + [ 0.K* du + | c K* (u, 0) du a c K*(u, 0) du 
0 x—1 n 


1 z 


=0+0 +f 6 Pant) au +e P,,{u).e% du. 
2 0 


Nun setzen wir 


also 


Damit wird 
Pola) e* = ¢ Kf) — ¢ Ry, (2) e* + oR, (2) 6 — cK, (0). 
Die Glieder mit e°* ergeben 
F(z) = ¢ R,(2) SOR, Y) « 
Wenn wir noch den Zusammenhang von P, R und R’ beriicksichtigen: 
Ri(2) = ¢ R,,_.4(2)« 
Ber Vergleich der Konstanten ergibt 


; R,,(0) = Recah hi} : 
So erhalt man rekursiv 


R,(z) = —cz+& 


Pz) FB ce’ =e, 


R,(2) = 5 8 -cez—ce +e ie 
Pye) = at c8a(@—1) bee — (0+ ee, 


ape s = = . ter ara, Gea Sel ditey Stan) ter Se etimy et es. as aay cate oe er 
ce) 6 . o- . . 
a . . . . ae 2. 


Auf diese Art kann K* etappenweise berechnet werden. 


= 
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e) In gewissen Fallen ist auch die Kombination von § 5 d) und e) vorteil- 
haft. So kénnte im vorigen Beispiel, da K = K (x — y) ist, die Laplace-Trans- 
formation gebraucht werden. Bereits durch M(0) = ¢ erhalt man nach (13) die 
Normierung 
Cc ak 


ap SO Sse ee 


/ K*(x, 0) dx = } 
0 

Kombiniert man die Methoden, so erhalt man schrittweise die exakte 
Dichtefunktion /(x). 

Der Vorteil der Laplace-Transformation liegt darin, dass die Normierung 
und die Momente leicht exakt erhalten werden. Mit der Methode von § 5 e) 
erhalt man die Normierung — wenn «x nicht beschrankt ist — erst asymptotisch, 
dafiir die (allerdings noch nicht normierte) Dichtefunktion exakt. Der Haupt- 
vorteil liegt jedoch darin, dass sie auch fiir den Fall gilt, wo K nicht Funktion 
von x—y ist. Dies kommt besonders zur Geltung, falls a(x) sprunghaft 
variiert. Dies tritt ein, wenn in einem Kraftwerk die Anzahl der in Betrieb 
gesetzten Turbinen vom Wasserstand des Stausees abhangt. Auch bei einem 
Schalter kann dies vorkommen: ein Beamter bedient, wenn die Schlange kurz, 
zwei, wenn sie lang ist. 


3. Kapitel: L6sung durch die Laplace-Transformation 
§ 7. Herleitung der Differentialgleichung 


Wir tiben auf (3) die zweiseitige Laplace-Transformation aus. 


Ji a+ ) ods = — [es dx [ues t) H(y, x,t) dF(y,t). (24) 


w, H und die Lésung F seien so, dass in einem Bereich B der s-Ebene die verwen- 
deten Integrale existieren. Nun setzen wir voraus, dass H eine Funktion von 
x — y und ¢ ist, dass also die Héhen der Spriinge nicht davon abhangen, von 
wo aus dieselben erfolgen. Denn unter dieser Voraussetzung erhalt man rechts 
eine Faltung. Nicht angegebene Integrationsgrenzen sind in diesem Kapitel 
(ausser § 11) — co und +o. Die rechte Seite von (24) wird, abgesehen vom 
Vorzeichen, 
i ie H(x — y,t) e@-» 84 w(y, t) dF(y, t) dx = “ H(z, t) 8 dz i‘ w(y, t) 6” dF (y, t) . 
Damit wir die Transformation der Produkte bilden k6nnen, setzen wir zu- 
ndachst voraus, dass 

N 


—#% = DS) aelt) x, wy, 1) = DP a). 


of 
| 
: 
7 
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Wir definieren nun 


M(s, t) = / e** f(x, t) dx, h(s,t) = wi dH (x,t) , (26) 
und setzen voraus, dass in einem Teilbereich B’ von B fiir alle ¢ 
lim e**F(x)=0 und lim e A(x) =0. (27) 
X—>+ co X¥-+> -£ 00 


Somit erhalt man durch partielle Integration 
if e® F(x) dx = —~ M(s) , if eH (x) dx = —~ 2(s) . 
M bedeute die partielle Ableitung von M nach ¢, M(”) die n-te nach s. So ist 
M(s, t) = [x Bea t(wer) dx é 
Damit erhalt man die Gleichung 


~¥ a, Ms) 2 M == Yp, Me 
_ oder 
N 3 
Dd, M's, t) [an(t) s + Pal) Ws, t)] + M(s, t) = 0. (28) 
n=0 
Man hat also hier eine partielle Differentialgleichung N-ter Ordnung fiir 


M mit der Anfangsbedingung, dass M(s, 0) die Laplace-Transformierte von 
f(x, 0) ist. 


§ 8. Lésung der Gleichungen erster Ordnung 


Im Fall N = 1 setzen wir — x(x, t) = a(t) + d(¢) x und w(y, t) = A(t) + g(t) y. 
- Somit lautet die Differentialgleichung 
[a(t) s + p(t) hs, t)] M + [B(t) s + gf) Als, t)] M+ Ms 0. (29) 


Fir K = In M, die erzeugende Funktion der Kumulanten (Semiinvarianten), 
lautet die Gleichung 


3 (ast+ph)+ (bst+qh)K'+K=0. (30) 
Wir fiihren statt s eine neue Variable wu = u(s, ¢) ein. Es ist 
dK OK , . 0K / du OK 
: bia OG bx ten gO (aK). > OK* ou 
lawl see dak me Rame ere du ds * 


ee ad 


- ZAMP. VII/19 © 


* 


ieee Bill | 
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In die Gleichung fiir K eingesetzt 


OK 


Ou Ou OK ae 
stants tae (bs+ 9H + 3p) + 


Or ae RS 
Der Faktor von 0K/du kann zum Verschwinden gebracht werden, wenn fur 
u = const die Lésungen der gewohnlichen Differentialgleichung 


ds 
oe eed (31) 


beniitzt werden. Man erhalt somit mit s = s(u, ¢) und wu =s fir ¢=0 
a 
K(s(u, T), T) = al [a s(u, t) + ph(s(u, t),#)]dé+K(w,0). (32) 
0 


Die Integration ist langs w = const auszufiihren. Am Schluss ist die Koordi-_ 
natentransformation von s nach ~ wieder riickgangig zu machen. 

Zuniachst einige Spezialfalle. . 

a) Wenn q fiir alle ¢ verschwindet, w also nicht von y abhingt, so ist die 
Koordinatentransformation von h, und damit von G unabhangig. Man erhalt 


ds _ 
dina 
Wir setzen 


F 
[ole ax = i b(x) dt 
S=ue° , also. S.se'S- ee" 


Somit erhalten wir 


T _ usar 
K(S, T) =—[plt)h See aie ie 


T T 
fo — [0(x) ar — [o(c)ac (33) | 
ai) se! dt KN. S en ° ,0 | 


=k, + kK, +K,. 


Ore 


kK ist eine Summe von erzeugenden Funktionen von Kumulanten. Also ist 
F(x, t) eine Faltung von drei unabhangigen Verteilungen. 


EG hangt von S,T,b,pundh, aber nicht von a und F (x, 0) ab. 
is hangt von S, T, a und b, aber nicht von #, h und F(x, 0) ab. 
K, hangt von S, T, db und F(x, 0), aber nicht von a, bp und h ab. 
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Wenn wir p = 0 setzen, so bleiben K, und K, wngedndert. Diese Funktionen 
geben also nur an, wie die Aifahecvertethene langs den Lésungen der Differen- 
tialgleichung (1) gefiihrt wird. Der Einfluss der-Spriinge wird durch kK, dar- 
gestellt (vgl. Lévy [7], § 12,1°). 

b) Die Differentialgleichung (31) lasst sich auch noch leicht lésen, wenn 


6, q und h nicht von ¢ abhangen. In diesem Fall kann die Gleichung separiert 
werden. 


ds is ds 
PT tea 1S cara g eeiater (34) 


Durch Umkehrung dieser Gleichung erhalt man s(u, ¢). Dies setzt man in (32) 
ein, integriert und driickt dann wieder w durch s und ¢ aus. 
c) Die dussern Bedingungen seien konstant. Wenn eine stationare Verteilung 
existiert, so ist AK = 0, und man erhalt 
zre\ _ 4S + ph 
K's) = ee (35) 
Bei der Lésung tritt eine Integrationskonstante auf. Diese kann nicht aus 
der Anfangsbedingung erhalten werden, wie bei der Losung der partiellen 
Differentialgleichung, da F(x, 0) im stationaren Fall nicht bekannt ist. Da 


aber K die erzeugende Funktion der Kumulanten einer Wahrscheinlichkeits- 


verteilung ist, muss K(0) = 0 sein. Somit wird 


; a2 tp bls) 
K(s) = — b= tigke az). (36) 


Nun zuriick zum allgemeinen Fall. Wir lésen das Problem ohne den Umweg 


_iiber die Variablentransformation mit Hilfe der Taylor-Reihen in bezug auf s. 


aA 


ln eed ee S| 


Wir bilden 
(as+ ph) = 3" 4,(t) 8", (bs+qh) = Sy Ball | 
AE (37) 
M(s, t) = oe m,(t) s” , sr) a 1 Yala ide a4 
n=0 n=0 
n!m, (=M,,) ist das n-te Moment von F(x, ¢), insbesondere ist also m)= 1. 


pa! ky an die entsprechenden Semiinvarianten (Kumulanten). Da H die Diffe- 
renz der Verteilungen E und G ist, wird 


h(0) =0, n! h()(0) = —n-tes Moment von G. 


Also 
A,=0, A,=4-— p-mittlere Sprunghohe, 


n! A, = — p: n-tes Moment der Sprunghdhenverteilung G. 
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Analog mit den B,,. Man kann diese in (29) einsetzen und erhalt sukzessive 
lineare, inhomogene Differentialgleichungen fiir die Momente MN, (t) , Malt) »<tald 
Da in (30) (as +h) als Summand und nicht als Faktor auftritt, werden die 
Differentialgleichungen fiir die Kumulanten etwas einfacher. Durch Koeffi- 
zientenvergleich erhalt man 


Aq) + 3! Ball) hea (+1) + hy =0. (38) 


n = 1 ergibt die Gleichung 
ak k, By +A,=0. 


Dies ist eine gewodhnliche, lineare, inhomogene Differentialgleichung erster 
Ordnung fiir 2, (t). n = 2 gibt 


hy 2h By ey, Ba eon 


Der Ausdruck in der Klammer ist nach dem vorigen bekannt. Somit kann 
auch ,(t) berechnet werden. Solange die Momente von G existieren, erhalt man 
lineare Differentialgleichungen fiir k,, ky, Rs, .... Die Anfangswerte f,,(2) sind 
durch die Semiinvarianten von F(x, 0) gegeben. 


§9. Die Gleichungen hoherer Ordnung im stationaren Fall 


Wir betrachten zunachst den Fall N = 2 von (28). Es sei 
wy)=p+qyt+try?, —a(*x)=a+dxn%+cx?. 
Somit erhalt man im stationadren Fall 
M" (cs+rh)+M’' (bs+qh)+M (as+ph)=0. (39) 


In der Lésung treten zwei Integrationskonstanten auf. Diese miissen so 
gewahlt werden, dass M die Laplace-Transformierte einer Verteilung ist. Die 
Bedingung M(0) = 1 liefert sofort eine der Konstanten, die ausdriickt, dass die 
gesamte Wahrscheinlichkeit = 1 ist. Die andere muss zum Ausdruck bringen, 
dass die Dichtefunktion /(*) = 0 ist. Notwendige und hinreichende Bedin- 
gungen, dass M(i s) eine charakteristische Funktion ist, sind bekannt. Die- 
jenigen, welche sich auf globale Eigenschaften beziehen, kénnen nur dann_ 
ohne weiteres angewandt werden, wenn (39) geschlossen integriert werden kann. 
Im Falle, dass die Differentialgleichung durch eine MacLaurin-Reihe integriert 
wird, lassen sich diese Kriterien nicht gut anwenden. Hingegen kénnen unter 


dieser Annahme wenigstens Ungleichungen fiir die Koeffizienten dieser Reihe 
hergeleitet werden. 
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Wir beniitzen wieder den Ansatz (37) und die analoge Gleichung 


cstrh= S"C, 9 
n=0 
Wir fiihren dies in (39) ein und setzen den Koeffizienten von s” = 0: 


oan [@ + 2) @ +1) miza Cre + (G+ 1) mir Brit im Ani=0. (40) 
t=0 
Dabei ist m) = 1, wahrend m, = m noch bestimmt werden muss. Wir gewinnen 
Ungleichungen fiir m, wenn wir diejenigen fiir die absoluten Momente (siehe 
zum Beispiel CRAMER [2], Seite 176) beniitzen. Wenn wir folgern kénnen, dass 


_ x sicher nicht negativ ist, so stimmen die gewdhnlichen Momente mit den 


absoluten tiberein. Es muss also, wenn 


iM, =n! om, 


_ die Momente von F(x) sind, gelten 


M; = Ms M,,-1 . (41) 
Diese Ungleichung ist scharfer als die 6fters gebrauchte 
Mun = Mieth) 


und fiihrt in unserm Fall iiberdies auf einfachere Rechnungen. Beispiel: 


2 he i sah Me) = 


Wenn alle Spriinge die Héhe 1 haben, ist 


33 


= pea 


2 
—h=s+ oa I 
Man erhalt also aus (39), 
sM"+(2M+M’)h=0, 


mM'=—(2M+mM)4, 
s2 
M, + Mys+ MS +++ = (2+m+(2m+M,)s+(2M,+M) > +-'] 
Ss Ga 
x{1+ 54+ 3°"), 
M,=2 +m, 
m ib ot 

My = 2 1 My eed oe ae re NS, 

2 5 


+8—m, 


2 
M,= + 77 +2m+M,+2M,+M=93 6 


Nacht he! wires nate. o) ah. 610th Lame) cele. “oh B.sele Re UMeh Chis. Rappe: SSH ers she ieuie Ra eee ee ee ae 


294 ANDREAS DALCHER ZAMP 
Nun bilden wir die Ungleichungen 
M,—m?=2+m—m20 (mS2), 
M,m — M3= 1m? +3m—4—4m—m?20, 


5m*—2m—820 (m=1,48). 


M,M,— M2 =0 und M,M,— Mf 20 geben keine Verscharfung der Un- 
gleichungen fiir m, hingegen eine Relation zwischen den absoluten Zentral- 
momenten: 4, = o4. Allerdings ist dies eine Gleichung vierten Grades in m. 


fy — 04 = M,—4mM,+ 6m? M,—3m*—(2+m—m’*)? 20, 
52-7 m+ m+ 8mi—4m*>0 (m<1,83), 
also 1,48 < m < 1,83. 
Wenn / rational ist, so werden nach Multiplikation mit dem Nenner die 


Klammern in (39) Polynome in s. Man erhalt so eine Rekursionsformel fiir die 
Momente M,,. 


Beispiel. Die Sprunghéhe sei exponentiell verteilt. Also ist H = e~** fiir © 
x = 0, und es wird 
1 ra —1 
eee Br we (s—k) a. 
5 As) [e a hi ara 
0 
Ss 
ils) = (42) 


Wir dividieren (39) durch / und erhalten 


M"(cs—ckh+n+M'(bs—bk+q)+M (as—ak+})=0. | 


Auf die Integration dieser Gleichung durch die Methode von LAPLACE 


kommen wir in § 11 zuriick. Hier verwenden wir die Momente. Wir leiten | 
n-mal ab und setzen s = 0, 


M,+2(7—¢k) + M,41(9—bk+ cn) +My(p—ak+5n) + M,,an=0. 
Es sei nun 


k=1, —#=5x%+x2; w(y)=10+y. 
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Setzen wir diese Werte ein, so erhalten wir 


Mn+2 = Mai (1 — 4) + M, (Sn + 10) , 


M,=1, 
M,=m, 
M,= 10-4 m, 


M, = —304+ 27m, 
M, = 260 — 134m, 


anew a) (Gee) Gee, gab ke SP. Gh.e terete ee (wie: kh: ah ene” we 


Da hier x = 0 ist, kénnen wir (41) anwenden. 
m=1: 10— 4m —m* 20 (a S 1,742), 
m= 2: 11m? +50m—100 20 (m = 1,503), 
nw = 3: 1700 — 760 m— 193 m2? =O (m <1,593). 
Die folgenden Ungleichungen ergeben 


m = 1,544, m=<1,578 und m= 1,547, 
also 
1547 <m <1,578. 


Die Ungleichung py, = o* liefert hier keine Verscharfung. Dass die Ein- 

schrankung von m hier viel besser ist als im ersten Beispiel, liegt daran, dass 
_ hier x(x) weniger von einer Geraden abweicht als dort. 
Anhand eines Beispiels wollen wir noch den Fall diskutieren, dass % ein 
- Polynom dritten Grades ist. Dies fiihrt auf eine Differentialgleichung dritter 
- Ordnung. Wir beschranken uns wieder auf den stationaéren Fall. Das Problem 
_ wird mit der Momentenmethode behandelt. Es sei 


—#=2%—2x7+ x3, hasty Sten, w=1. 
- Wir setzen dies in Gleichung (28) ein und erhalten 
—M(s) (} +24...) 42M) — 2M") + MM) =0. 


| Eine Integrationskonstante ist M(0) = 1, die beiden andern x = M ‘(0) und 
_y=M"(0), sind, wie m in den vorherigen Beispielen, noch unbekannt. Wir 
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bilden die Taylorentwicklung der obigen Gleichung und erhalten 

1 
—(l+sx+---) (> eg 2 a ++) +2(x+ sy +++) 


—2y+sM,+---)+(Mat+sM,+-:)=9, 


1 
—~+2%-2y+M=0, M, = 29-24 ti (43) 
be ae 7 7h 
sip aame 2M,+M,=0, Mase 2.) wed ois 


Nun betrachten wir die Ungleichungen (41) fiir die Momente. 
M,M,—M?20, y—-#x20, 
M,M,—-MZ20, 2xy—2074+5-—y?20, 
M,M,— M320, 2y? —~fLaytty—(2y-224 5) iy Bee. 
1 


—2y? 408+ 2 wy—2y+2x— 


C0: ete: 08 (ee aye ee) ois oa © © 06 ee 0 (elle ‘e @ ce 6 6\'6 10 © 6 is) «© 6) 4) (6) elle je: e)e) susuele io) «) hele a) ames) 


i 


In Figur 6 sind die Parabel und die beiden Ellipsen eingezeichnet. Bereits 
aus der Tatsache, dass die mittlere SprunghGhe 1/2 ist, folgt, dass (x, y) in der 
grdssern Ellipse liegt. Man erhalt die Ungleichung M, <1, unabhangig von der 
Streuung of der Sprungverteilung: diese braucht nicht einmal zu existieren! 
Analog folgt aus of = 2/6 — (1/2)2, dass (x, y) in der kleinern Ellipse liegt und 
dass M, < 0,6 ist. Auch dies gilt unabhangig von den héhern Momenten von G. 
Dass F mehr Momente als G hat, hangt mit dem Umstand zusammen, dass x(t) 
eine senkrechte Asymptote hat. Dadurch werden extrem grosse x-Werte rasch 
verkleinert. 

Wenn noch mehr Momente von G bekannt sind, erhalt man andlor weitere 
Kegelschnitte. So wird das Gebiet, in dem (x, y) liegen kann, weiter verkleinert. 
Dieses schraénkt sich jedoch im allgemeinen nicht auf einen Punkt zusammen. 
Dies riihrt daher, dass die Ungleichungen (41) nicht hinreichend sind, damit 
M(s) die Laplace-Transformierte einer Wahrscheinlichkeitsverteilung ist. Man 


kann die Beziehung noch verallgemeinern, indem man statt » ein beliebiges 
reelles z setzt. So erhalt man 
M,,,M,_,—MZ=0. 


+e 


Der Beweis erfolgt genau gleich wie bei ganzen (siche CRAMER [2], Seite 176). 
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Yea 


Figur 6 


Die durch Gleichung (44) dargestellten Kegelschnitte. 


Auch diese Ungleichungen geniigen nicht, dass die M, die Momente einer 
Wahrscheinlichkeitsverteilung sind. Wir zeigen dies durch ein Gegenbeispiel. 


Es sei 
Mla 1 tars = 2 <2. Mls) =2e-") far 22:2). 


Wenn ein F mit diesen absoluten Momenten existiert, so ist insbesondere 
[aE (x) =1, [|x|aF()=1, [x2 aF(x) =1. 


Die Verteilung von |x| hat die gewéhnlichen Momente M)= M,= M,=1, 
und somit ist deren Streuung = 0, also | «| = 1 mit W = 1. Dies steht aber im 
Widerspruch zu den angegebenen héhern absoluten Momenten. 


§ 10. Partielle Gleichungen zwetter Ordnung 
Im Fall 


— (x, t) = a(t) + b) x + lt) x7, w(y,t) = de) + 9) y+ (f° 
erhalt man nach (28) die partielle Differentialgleichung 


M'(cs+rh)+M'(bs+ qh) +M(ast+ph)+M=0. 
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Noch einen andern Fall kénnen wir auf eine Differentialgleichung zweiter 
Ordnung zuriickfiihren. Es sei 


w(y, t) = pt) + at) y, 


x sei aber kein Polynom, sondern von der Form 


— (x,t) = a(t) — 2, (45) 
Wahrend die vorherigen Probleme in diesem Kapitel dem Charakter des 
Beispiels d) von § 2 entsprechen, hat dieses den einer «Queue» (Schlange): 
x = 0 fiir ein festes x und annahernd konstant fiir grosse x. 
Um in diesem Fall die Differentialgleichung zu erhalten, multiplizieren wir 
(3) mit x e** und integrieren von —oo bis +00. 


—ff (ax —b)etdx+ rf Fe) est dx 


I 


s / dx /@ + y) x H(x— y) e f(y) dy 


— ff + ay) eee ((x — 9) + y) Hx — 9) 10) dy ae, 
Somit erhalt man 
, 0 (-1 
—aM'+bM+— (=m) 


a ai 


Wenn wir noch mit s multiplizieren und zusammenfassen 


. h) M— —* hh M’ (= h) M’—q——hM". 


Ss 


M'+qhM" + [as+ph+ (x — +) | M’ 
fon i (46) 
-—M+ E (2 -=) —bs|M=0. 
Im stationaren Fall lasst sich diese Gleichung nach der Methode von § 9 


losen. Im allgemeinen jedoch ist dies auch bei einfachen Verteilungen G eine 
schwierige Differentialgleichung. 


§ 11. Stationarer Zustand bei exponentiell verteilter Sprunghohe 
Es sei 
Gaya l= eke fir a4 20) 1Gia) sai ieee 
Damit wird nach (42) 
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Wir setzen dies in (28) ein und dividieren durch h. Im stationiren Fall 
erhalt man somit 


N 
$5, M) (s) (a, S— ay k+ Pn) = 2 a) 
n=0 


Dies ist eine Laplacesche Differentialgleichung. Sie wird gelést durch den 
Ansatz 


M(s) = / 8? f(x) dx. (48) 


Dabei sind /(x) und die Integrationsgrenzen geeignet zu wihlen. Setzt man den 
Ansatz in die Differentialgleichung ein, so erhalt man 


a 


fi sD) ay x” f(x) dx + | Be DN (Py ok a) #° xy dx = 0. 


Die Summen sind die Polynome — x(x) und w(x) + k x(x). Das erste Integral 
wird nun partiell umgeformt. 


— [esx fa) dx = =e" & f(x) + [er (a f(x))' dx. 
Somit erhalten wir 
soe Hh | 4 i= (a f(x)’ + (w+ hx) Ma) dx =0. (49) 
Wir wahlen nun f(x) so, dass die eckige Klammer verschwindet : 


x f(x) + f(*) (@ +w+ks«)=0. 


f’(*) ye wtkx 
eal ares 
Sh +w(u)/x(u) du 
fe mr hs (50) 


Nun werden die Integrationsgrenzen von (48) [nicht die im Exponenten von 
(50)] so bestimmt, dass e*” « f(x) in diesen verschwindet. Also muss in den 
Grenzen X gelten rs 
= g k+w/sdu 
4(X) f(X)=ce* == (0) (51) 


Die Gleichung gelte in X (= —oo) und in +00, aber nicht im Intervall 

X < x < oo. So ist (47) erfiillt durch (48) und (50). Die Integrationskonstante 

‘e wird durch M (0) = / f(x) dx =1 erhalten. Da M(s) die Laplace-Transfor- 
= . 

-mierte der Verteilung ist, haben wir mit f(x) direkt die gesuchte Wahrschein- 


> 
7 
, 


‘ 


: 
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lichkeitsdichte gefunden. Der Umweg tiber die Laplace-Transformation konnte 
also mit der Gleichung (50) abgeschnitten werden. 

Bemerkenswert ist, dass in der Formel (50) gar nicht zum Ausdruck 
kommt, dass w und % Polynome in x sind. Wir lassen nun beliebige Funktionen 
zu. Wir bilden f(x) nach (50) und miissen priifen, ob durch diesen Ansatz 
Gleichung (3) erfiillt ist. Es ist also zu zeigen, dass 


Die linke Seite ist 


x 
~ [k+wlidu 
nx =le" 
Rechts erhalten wir 
y 
x 4% — fk+wjzdu 
Cc 
— [w(y) e-984(y) dy = — fol) ew-* dy 
fi 
a — | w/xdu 
ee =o f eva) e(% —y) k oe e * dy 
x 
ae ~ fojiau 
= —¢ el—“)k poe dy 
y x 
— fojxau 


Wegen (51) verschwindet der Ausdruck an der untern Grenze. Somit stimmen 
die beiden Seiten der Gleichung iiberein. 

Da alle Spriinge positiv sind, kann x nicht kleiner als z werden, falls x(z) 0. 
Eine Untermenge von x < z kénnte zwar ergodische Klasse sein (siehe DooB [3]). 
Da aber beliebig hohe Spriinge méglich sind, ist dies (ausser x = —oo) nur 
moglich, wenn darin w = 0. Sonst ist im stationaren Fall « > z. 

Ferner beachten wir, dass (50) mit (20) iibereinstimmt bis auf die untere 
Grenze des Integrals und die damit zusammenhangende Normierung. In (20) 
ist die Grenze 0, also die untere Schranke der x-Werte. In (50) wiirde nach (51) 


das Integral mit X als unterer Grenze divergieren. Deshalb wahlten wir eine 
andere: %,. 


Die Divergenz des Integrals 


x xX 
w w 
[yu | wan dx = | wat =00 
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bedeutet, dass mit W=1 x >X bleibt, wahrend nach (9) x = 0 mit positiver 
Wahrscheinlichkeit angenommen wird. Die Lésunmgen (20) und (50) unter- 
scheiden sich also durch die Verhaltnisse an der untern Grenze. 


Wenn * und w in x, verschwinden, so bleibt der Prozess in x, konstant, 
sobald dieser Punkt erreicht ist. Wenn 


*y 


IF dx =[w dt 


konvergiert, so geschieht dies mit positiver Wahrscheinlichkeit ; x, ist in diesem 
Fall absorbierende Schranke. 


4. Kapitel: Diverses 
§ 12. Ansatz fiir den mehrdimensionalen Fall 
x ist ein Vektor (x,, %, ..., %»)- Die entsprechende Differentialgleichung 
lautet 


oder, skalar geschrieben, 


. _ ax; 


x; = di SM, Bon te ge Xe the 4 Ne Dees eS 


Dies umfasst auch den Fall einer Differentialgleichung u-ter Ordnung. 
Wenn diese skalar ist, so setzen wir 


% = (%,%,X,..., x! 
« FA o < 5 eee Fi 
4i= X41 fiir t<n, Xn = K(X, x, aaa . t) . 


Auch die Anfangsverteilung, die Dichte der Spriinge und die Verteilung 
nach diesen sind Funktionen mehrerer Variablen: 


F(x,t), wly,t), GY, x,t). 
‘Formal andert sich bei der Ableitung der Hauptformel (3) nichts. 


; § 13. Berechnung von Naherungen 

Auch im nichtstationaren Fall lasst sich f(x, t) naherungsweise nach der 
folgenden Methode berechnen. Wir beschranken uns zundichst auf den Fall, 
‘dass « nicht von ¢ abhangt. x = 0 sei untere Grenze. Wir teilen die Zeit in 
gleiche Intervalle At. Dann werden die in den Momenten ¢; = 7 At null werden- 
den Lésungen der Differentialgleichung (1) eingetragen (siehe Figur 7). Durch 
diese und die Senkrechten ¢ = ¢; erhalt man ein Netz von Schnittpunkten, und 
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Figur 7 
damit die Ordinaten x) = 0, %4, %_)--+, %n, +... Wir setzen 4; = ¢;+ Az/2 und 
AF, inner a F(x, ti) a F(%, 4 t;) , 
é 1, Reale 
Gg ng = G(%p) %pr Gi) — G(s %q 1b) 


(Fiir die Definition von %;, siehe Figur 7.) 
So wird 


F(0, ¢;,,) = (1 — w(0, t}) At] F(x, ¢) + F(x, §) w(0, 3) AtG(O, Gast) 
a ay AF, ¢3,3 w(x,, t;) At G(x, Xt) , 
R 


AF 


nt+1 


= [1 — w(x, t;) At] AF, .1 + F(a, &) w(0, tj) Ae AG ni 
nie DS: ARS aS W(X, t;) At AG, 4, : 
h 


So erhalt man fiir diese 4F, ;,, lineare Kombinationen der AF, ;, die sukzes- 
sive gelést werden kénnen (zum Beispiel mit Rechenautomaten). 

Wenn «x mit ¢ variiert, so werden die x; von ¢ abhangig. Im iibrigen verandert 
sich jedoch die Methode nicht. Statt w(x,, ¢) kann auch das Mittel von w im 
Viereck 

(Xe sae 4) (%ps ti) (Xpoas tga) (Mp ts 41) 
genommen werden. 
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Wenn das Problem auf eine Integralgleichung gebracht werden kann, so 
lasst sich diese durch ein System linearer Gleichungen annihern. Da es sich 
um Volterra-Gleichungen handelt, weist die Matrix Dreiecksform auf. Das 
(zunachst noch nicht normierte) f(x) kann rekursiv erhalten werden. Diese 
Methode gibt weniger Rechenarbeit als die erste, ist aber nur fiir den stationaren 
Fall brauchbar. Auch ist nicht ersichtlich, wie schnell sich die Anfangsvertei- 
lung der asymptotischen niahert. 


§ 14. Anmerkungen tiber Existenzsatze 


Zunachst muss die Differentialgleichung (1) in jedem Punkt (x, ¢), der vom 
Prozess angenommen werden kann, eine eindeutige Losung haben. Dies hat 
wenigstens fiir wachsende ¢ zu gelten; fiir abnehmende ¢ ist die Bedingung bei 
den Problemen des 2. Kapitels in x = 0 nicht erfiillt. 

Nun wenden wir uns der Frage zu, ob das Hauptproblem eine eindeutige 
Lésung habe. Wenn x = x(c, t) die Lésungen von (1) sind, so kénnen wir eine Ko- 
ordinatentransformation durchfiihren mit ¢ und der Integrationskonstanten c als 
neue Variable. Damit sind die einzigen Veranderungen von c die Spriinge. So ist 
das Problem auf den «absolut unstetigen Fall» von W. FELLER [4] zuriickgefihrt. 

Aber auch wenn die Existenz einer Lésung bewiesen ist, braucht es keine 
stationare zu geben. Im «Queue»-Problem mit —x =a und w= const ist 
bewiesen worden (zum Beispiel LINDLEY [8]), dass eine stationdre Lésung 
existiert, wenn ~G < a ist. In der Integralgleichung dussert sich dies darin, 


dass [BA(x, 0) dx konvergiert. Dies ist nicht auf den Spezialfall = const 
‘ 0 


und a = const beschrankt. Denn falls dieses Integral konvergiert, so erhalt man 
mit (9) ein F(x), das die Integralgleichung und damit auch (3) befriedigt. 

_ Analog ist es mit der Laplace-Transformation. Wenn wir eine stationare 
oder variable Lésung der Differentialgleichung (28) gefunden haben, so ist die 
Umkehrung der Laplace-Transformation eine Lésung von (3). Wenn diese eine 
Wahrscheinlichkeitsverteilung ist, so ist mit der Lésung des Problems auch 
deren Existenz bewiesen. Weiteres iiber Existenzsatze findet sich in BLANC- 
‘LAPIERRE und Forrtet [1] und in Doos [3]. 


Es ist meine angenehme Pflicht, Herrn Prof. Dr. W. SAXER, Ziirich, sowie 
‘auch Herrn Prof. Dr. A. PFLUGER, Ziirich, und Herrn D. V. LInDLEy, Lecturer, 
Cambridge, fiir die wohlwollende Unterstiitzung meinen aufrichtigen Dank aus- 
zusprechen. 

LITERATURVERZEICHNIS 
[1] A. Branc-Laprierre und Fortet, Théorie des fonctions aléatoives (Masson Cie., 


Paris 1953). I ; 
[2] Haratp Cramér, Mathematical Methods of Statistics (Princeton University 


Press 1946). 


roe AEN 


ie el 


304 WitHetm KLINGENBERG ZAMP 


[3] J. L. Doon, Stochastic Processes (John Wiley and Sons, New York; Chap- 
man and Hall, London 1953). 
[4] Witty Fetter, Zur Theorie dey stochastischen Prozesse (Existenz- und Ein- 
 deutigkeitssatze), Math. Ann, 773, 113-160 (1936). ii 
[5] E. Goursat, Détermination de la résolvente d’une classe ad’ équations intégrales, 
"” Bull. Sci. math, 68, 114-150 (1933), 
[6] Davip G. KENDALL, Stochastic Processes in the Theory of Queues, Ann. Math, 
Statist. 24, 338 (1953). 
[7] Pau Livy, Processus stochastiques et mouvement Brownien (Gauthier- 
Villars, Paris 1948). 
[8] D. V. Linprey, Theory of Queues with a Single Server, Proc. Camb. phil. 
Soc. 48, 277 (1952). 
[9] WERNER SCHMEIDLER, Integralgleichungen mit Anwendungen in Physik und 
Technik (Geest und Portig, Leipzig 1950). 
[10] Laurent Scuwartz, Théorie des distributions (Hermann, Paris 1950). 


' 


Summary 


Stochastic processes of the following type are considered. At random time 
points ¢, the variable *(t) jumps from y to x, say. The heights + — y of the jumps 
have a given distribution G*(+ — y) that may depend on y or ¢. Between the 
jumps, +(f) is a solution to a given differential equation dx/dt = a(x, t). We look 
for the distribution F(¥, t) of x at time ¢> 0, F(x, 0) being given. In the stationary 
case, stable distributions are investigated. 

If there is a lower boundary ¥, and if F(*%,) > 0, the problem is similar to the 
queueing problem. We solve it in the stationary case with integral equations of 
the Volterra type. Other problems can be transformed to differential equations 
for the moment generating functions. These equations are partial in the non 
stationary and ordinary in the stationary case. 


(Eingegangen: 16. September 1955.) 


Die Anzahl det Nullstellen eines Polynoms 
in Gebieten mit stiickweise rationalen Randkurven 


Von WILHELM KLINGENBERG, Hamburg?) : 


Wir fecacmicee ein Polynom w(z) mit komplexen Koeffizienten. Fiir viele 
Anwendungen ist es niitzlich, wenn man nach einem einfachen Verfahren die 
Anzahl der Nullstellen in gewissen Gebieten der komplexen Zahlenebene be- 
stimmen kann, ohne die genaue Lage der Nullstellen berechnen zu miissen. 

Falls das Gebiet eine Halbebene oder eine Kreisscheibe ist, so kennt man 
seit HERMITE, Hurwitz, Routu, FROBENIUS u. a. 2) mehr oder weniger vollstan- 


1) Mathematisches Seminar der Universitat, Institut fiir angewandte Mathematik. 
*) Fir ausfiihrliche Literaturangaben verweisen wir auf: M. MarvEen, The Geometry of thd 
Zeros of a Polynomial in a Complex Variable (New York 1949). | 
‘ 
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dige Lésungen des vorliegenden Problems; ein fiir numerische Zwecke besonders 
geeignetes Verfahren gibt WriELanpr’) an. Von allgemeineren Gebieten der 
komplexen Ebene scheint bisher nur der Winkeiraum oder das Kreiszweieck 
in einer Arbeit von SHERMAN‘) behandelt zu sein. 

Wir wollen in dieser Note zeigen (vgl. Satz 2 in Abschnitt 5 und die Er- 
ganzung in Abschnitt 6), dass man auch fiir eine bedeutend allgemeinere Klasse 
von Gebieten die Anzahl der Nullstellen eines Polynoms in einem solchen 
Gebiet in endlich vielen Schritten durch rationale Operationen und Grdssen- 


_vergleichung bestimmen kann; fiir den Spezialfall des Winkelraums und der 


Halbebene werden wir dabei auf das Verfahren von SHERMAN bzw. WIELANDT 
gefiihrt. Und zwar brauchen wir bei unserem Verfahren von dem Gebiet nur 
vorauszusetzen, dass es einfachzusammenhiangend ist und dass sein Rand sich 
aus endlich vielen Stiicken von rationalen Kurven zusammensetzt. 

Das Kernstiick unserer Methode besteht in einem Verfahren, die Anderung 


2x I= | dargw (1) 
© 
des Arguments von w zu berechnen, wenn w = w(¢) ein Stiick € einer rationalen 
Kurve durchlauft (siehe Abschnitte 3 und 4). Wir benutzen dabei ein Lemma 
(vgl. Abschnitt 4) tiber verallgemeinerte Sturmsche Ketten, das so formuliert 
ist, dass es vielleicht auch bei der Einsetzung noch allgemeinerer Kurvenstiicke 
€ in (1) von Nutzen sein kann. 


1. Funktionentheoretische Vorbereitungen 


Bei unserem Verfahren stiitzen wir uns auf den bekannten Cauchyschen 


Satz von der Windungszahl’). 


Es sei ¢ eine geschlossene, orientierte und rektifizierbare Kurve in der 
komplexen z-Ebene. 2, liege nicht auf c. Dann heisst die Zahl 


nz ¢) = (221) / (¢ — a))-1 dz (2) 


c 


“die Windungszahl von ¢ beziiglich 2p. 


Sei nun 
w(z) = Td (z— 2)" (m; ganz und +0) (3) 


7 


eine rationale Funktion, so dass keine Nullstelle oder kein Pol auf c liegt. 


3) H. WiELANDT, Beitrdge zur mathematischen Behandlung komplexer Eigenwertprobleme. 1: Ab- 


_azthlung der Eigenwerte komplexer Matrizen, Deutsche Luftfahrtforschung, Forschungsbericht 
Nr. 1806/1. 13 S. (1943). 


c 


- 


4) S. SHERMAN, Generalized Routh-Hurwitz Discriminant, Phil. Mag. [7] 37, 587-551 (1946). 
5) Vel. fiir das Folgende etwa: L. V. AHLFORS, Complex Analysis (New York 1953). 


ZAMP VII/20- 


306 WILHELM KLINGENBERG ZAMP 


© bezeichne das Bild w(c) von ¢ in der komplexen w-Ebene. Dann folgt aus 


(2) und (3) 
w- dw = pe m,(z — 2;)~1 dz 
sae 
oder 


n(0, ©) = yah mM; n(z;, ¢) - (4) 


Das heisst: Aus der Windungszahl n(0, €) der Kurve € = w(c) in der w-Ebene — 


beziiglich des Punktes w = 0 bestimmt sich die in (4) rechtsstehende Summe. 

Ist nun speziell ¢ der positiv durchlaufene Rand eines einfachzusammen- 
hangenden Gebietes J" der z-Ebene und z; ein nicht auf ¢ gelegener Punkt, so 
ist die Windungszahl n(z;, c) gleich 1 oder gleich 0, je nachdem z; zu J” gehért 
oder nicht. Wenn nun also w(z) eine rationale Funktion (3) ist, fiir die keines der 


2;2u ¢ gehdrt, so folgt aus (4): 2(0, €) mit C= w/(c) ist gleich der Differenz N—P — 


aus der Anzahl N der Nullstellen und der Anzahl P der Pole von w(z), die zu 
I’ geh6ren. 
Damit ist das Problem, die Anzahl der Nullstellen eines Polynoms w(z) in 


einem einfachzusammenhangenden Gebiet zu bestimmen, auf die Berechnung 


der Windungszahl einer geschlossenen Kurve zuriickgefiihrt. 


2. Die Berechnung der Windungszahl einer geschlossenen Kurve, die 


sich aus endlich vielen Stiicken von rationalen Kurven zusammensetzt 


Wir betrachten eine geschlossene orientierte Kurve € in der komplexen 
w-Ebene, die nicht die Punkte w = 0 und w =co enthalt und die sich aus endlich 
vielen Stiicken ©, (0 = 1,...,7) von rationalen Kurven zusammensetzt. Das 
bedeutet, dass jedes €, eine Parameterdarstellung 


w,(@) = “20 (a, StS, (5) 


gestattet, wo P,(t), Q,(t), R,(t) reelle Polynome sind, von denen wir annehmen — 


kénnen, dass sie keinen gemeinsamen Teiler + const besitzen. Die Vorausset- 
zung, dass © weder den Punkt w = 0 noch den Punkt w = oo enthilt, bedeutet: 


Ro) +0 und Pid)+ Qi) +0 fir te [a,, ig (6) 
Die Windungszahl (0, €) ist wegen 


n(0, ©) = (2 %1)-3 Des dw = (2x)-1 Sfa argw (7) 
’ & ° &, 


6 . 5 . . . 
) Wir benutzen [a, b] zur Bezeichnung eines abgeschlossenen und Ja, b[ zur Bezeichnung eines ~ 


offenen Intervalls mit den Endpunkten @ und b. Die zugehérigen links- bzw. rechtshalboffenen 
Intervalle werden mit Ja, b] und [@, bf bezeichnet. 
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durch die Anderung ; 
2x1, =| dargw (8) 


& 


des Arguments von w lings des rationalen Kurvenstiicks €, bestimmt. 
Offenbar hat das Integral J, in (8) denselben Wert, wenn wir die Kurve C, 
durch die Kurve €¥ mit der Parameterdarstellung 


. wet) = P,Q) +70, (t) (a, S¢Sd,) (9) 
ersetzen. 


3. Die Berechnung der Anderung von arg w langs eines orientierten 
Stiicks einer rationalen Kurve’) 


Wir wollen den Wert des Integrals (1) berechnen, wo € ein orientiertes Stiick 
einer rationalen Kurve mit der Parameterdarstellung 


w(it)= P()+7Q0(t) (@aStsb) (10) 
ist; P(t) und Q(é) sind hier reelle Polynome, nicht beide = const, und mit der 
Eigenschaft : ; 
P(t) + Q2(¢) +0 fiir te [a,b]. (11) 


Wir unterscheiden die beiden Méglichkeiten§) 


Fall 1: Grad P(t) = Grad Q@), | 
Fall 2: Grad P(t) <Grad g(t). | e 
Hiermit definieren wir 
im Fall 1: fo) = PW), fi) =—00), | a 
im Fall 2: 2) =), A()= PO, | 
“so dass also in jedem Fall gilt: 
: Grad f,é) = Grad f,(é) und Grad f,(t) > 0. (14) 


Wir teilen die Punkte w +0 und +oo in vier Klassen ein, je nachdem 
argw zu einem der Intervalle (0, 2/2[; [x/2, a[; [7,3 2/2[; [3 2/2, 2a[ gehért. 
Diese Klassen nennen wir beziehungsweise ersten, zweiten, dritten und vierten 
- Quadranten und bezeichnen sie mit Ox, Or, Qi und Qjy. 

7) In dieser und der folgenden Nummer unterdriicken wir den Index @. Insbesondere bezeichnet 
jetzt also € nicht, wie in Abschnitt 2, notwendig eine geschlossene Kurve, sondern nur ein Sttick 


_ einer rationalen Kurve. 
8) Wir vereinbaren, dass das Nullpolynom einen Grad <0 hat. 
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Die Anderung (1) von argw lings © wird davon abhangen, wie oft beim 
Durchlaufen von € Quadranteniibergainge im zunehmenden und im abneh- 


menden Sinn stattfinden. 
Beschrinken wir uns zundchst auf den Fall 1. Hier bedeutet eine Nullstelle ¢, 


von fo(t) [= P(t)], dass € die imaginare Achse trifft. Die Bedingung 


(A) Es gibt « > 0 so, dass ay < ty — e und dass fiir alle ¢ mit 4; -e<¢< 
gilt: fo(¢) filo) < 0 

besagt gerade: Es gibt eine dem auf der imaginaéren Achse gelegenen Punkt 

w(t)) vorangehende Umgebung auf ©, die in Q; oder Qy liegt; da die Punkte 

der imaginaren Achse zu Qj; oder Qyy gehéren, bedeutet also eine Nullstelle f, 

mit der Eigenschaft (A) einen Quadrantenwechsel im zunehmenden Sinn. 
Ganz ebenso erkennt man, dass eine Nullstelle ¢) von f(t) mit der Eigen- 

schaft 


(B) Es gibt ein e > 0 so, dass ty + ¢ < bg und dass fiir alle ¢ mit 
pete Ye. ote yey oO 


einen Quadrantenwechsel im abnehmenden Sinne bedeutet. ) 
Auf die gleiche Weise sind auch im Fall 2 die Bedingungen (A) und (B) 
gekennzeichnet, wobei nur an die Stelle der imaginaren die reelle Achse tritt. 
Wenn wir annehmen, dass A-mal eine Nullstelle von /,(¢) mit der Eigen-— 
schaft (A) und S-mal eine Nullstelle von f,(¢) mit der Eigenschaft (B) auftritt, so 
liefert offenbar die Zahl 


es S (15). 


eine erste Naherung fiir die Anderung von (2)-! argw langs . : 
Eine genauere Naherung erhalten wir mit Hilfe der Funktion N(é), die wir 
folgendermassen fiir die Werte ¢ aus [a, b] definieren: 


0, falls arg z(t) zu oder zu ehort, 

im Fall 1: N(t) = | BR ae ae 
| 1, falls argz(t) zu Q,, oder zu Q,y, gehért; { 
(16) 
O, falls arg z(t) zu d hort ; 
im Fall 2: N(t) = Se ie ae | 
1, falls arg2(t) zu Q,; oder zu Q,,, gehért. ; 
4 : 


pee (12), (13) ist diese Definition offenbar gleichbedeutend mit der fol-— 
genden: 


| 0, falls f,(t) =O oder f,(¢) 4,(t) <0, 
\ee (17) _ 


1, falls f(t) =0 oder f(t) A(t) > 0. 


| 
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Hiermit erkennt man dann, dass die Zahl 


RS ON 


NiO) 
-. = - 


. (18) 
gleich einem Viertel der Anzahl der gewonnenen oder verlorenen Quadranten 
ist, wenn man die Kurve € durchlauft. (18) stellt also, bis auf einen Fehler vom 
Betrag <2/2, die Anderung von (2 7)-! argw langs © dar. 

Den genauen Wert von (1) bendtigen wir fiir das in Abschnitt 2 aufgewor- 
fene Problem nicht; wir wollen ihn aber doch der Vollstaéndigkeit halber an- 
geben. Dazu bestimmen wir aus 


“a=argw(a) und f= argw(b) (19) 
die Werte x und f durch 
= t Bi J — Pp L 
od =amod7>, B = B mod (0 =a, p< 3). (20) 
Damit erhalten wir dann fiir (1) den Wert 
A= (270)-* / dargw = — ot) ACB ne oa Ae (21) 


4. Die Berechnung der Differenz R-S; ein Satz iiber 
verallgemeinerte Sturmsche Ketten 


Es bleibt fiir die Lésung der in Abschnitt 3 gestellten Aufgabe noch die 
Differenz R — S aus der Anzahl R der Nullstellen von /,(¢) mit der Eigenschaft 
_ (A) und der Anzahl S der Nullstellen von /,(¢) mit der Eigenschaft (B) zu 
bestimmen. 

Dazu erklaren wir mit Hilfe der Polynome f,(¢), /,(¢) eine Folge 


F(t) = (fo(t), ---» fit) 2 9) (22) 


von Polynomen nach folgender Vorschrift: 
Falls /,(¢) = 0, lassen wir F(t) mit f,(t) abbrechen. Andernfalls bestimmen 
_ wir nach dem euklidischen Algorithmus ein Polynom /,(¢) durch 


—f iss) = fal) + gilt) FAC) (23) 
mit z= 1; g;,(t) ist also ein Polynom, und es ist 
Grad f;,, = Gradf; und Grad /,,, < Grad /,_,. (24) 


_ Falls f,(¢) = 0, so lassen wir F(t) mit /,(¢) abbrechen ; andernfalls bestimmen wir 
- gemiss (23) mit 7 = 2 ein Polynom f,(¢). So fortfahrend gelangen wir wegen (24) 
nach endlich vielen Schritten zu einem Polynom /,,,(t) = 0; die Folge F’(¢) lassen 
_ wir dann mit dem vorhergehenden Polynom /,(t) abbrechen. 
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Mit W(t) bezeichnen wir die Anzahl der Vorzeichenwechsel in der Folge 
F(t) (22) an der Stelle t. Hiermit gilt dann die Beziehung 


R—S=W(a) — W(b). (25) 


Den Beweis fiihren wir mit Hilfe des folgenden Lemmas’) : 

Lemma. Es sei F(t) = (fold), --- » frlt)) (1 = 0) eine Folge von /+ 1 auf dem 
Intervall [a, 6] definierten und dort stetigen Funktionen mit folgenden Eigen- 
schaften: 

1, fit) = 0 fur f¢[a,0)- 
II. Die Anzahl der Nullstellen von f(¢) auf [@, b] ist endlich. 

Ill. Falls f,(¢,) = 0 fiir t)€ [a, b] und 0 <7 <1, so gilt f,4(to) fisalto) < 0. 

IV. Falls f,(¢y) = 0 fiir t) € [a, 0], so gilt f,(¢) + 0. 

Es bezeichne R bzw. S die Anzahl der Nullstellen von /,(¢) auf [a, 0} mit 


der Eigenschaft (A) bzw. (B) aus Abschnitt 3. W(f) sei die Anzahl der Vorzei- : 


chenwechsel der Folge F(¢) an der Stelle ¢. Dann gilt die Beziehung (25). 


Beweis des Lemmas. Wir verfolgen, wie sich W(t) als Funktion von ¢ beim ~ 


Durchlaufen des Intervalls [a, 6] andert. 


Offenbar bleibt W(t) in der Umgebung solcher Werte ¢ ungedndert, fiir die 


keines der /;(¢) aus F(t) verschwindet. Sei jetzt t) eine Nullstelle von /,(¢) mit 


a 


0<i< 1. Nach III besitzt dann die Folge (/,_,, /;, f;,1) an der Stelle ¢) genau ~ 
einen Vorzeichenwechsel. Da wegen der Stetigkeit /,_, und /;,, in einer Um- — 


gebung von ¢, ihr Vorzeichen nicht andern, erfahrt W(t) in dieser Umgebung — 


durch den Teil (/,;_,, ;, f::1) keine Anderung. 
Es bleibt wegen I jetzt nur noch der Einfluss einer Nullstelle ¢) von /o(?) 


auf W(t) zu untersuchen. Diese Nullstellen sind isoliert, und /,(¢) und /,(é) sind — 


stetig. Fiir t= 7%) besitzt die Folge (fo, /;) keinen Vorzeichenwechsel; in einer 


hinreichend kleinen Umgebung Jf) — é, t) + e[ von fy (soweit sie zu [a, 6] gehort), 


andert /,(¢) sein Vorzeichen nicht, und f,(¢) hat ausser t = ¢, keine Nullstelle. 


Die Bedingung (A) besagt nun nichts anderes, als dass die Folge (fp, /,) in 


der linksseitigen Umgebung ]¢) — ¢, fo[ von ¢, (soweit sie zu [a, 6] gehort) einen — 
Vorzeichenwechsel hat, das heisst: W(t) = W(t) — 1 fiir t€ ]éy — €, to[. Ebenso | 
besagt die Bedingung (B), dass die Folge (fp, /,) in der rechtsseitigen Umge- 


bung Jfy, 4) + e[ von ¢,( soweit sie zu [a, 6] gehort) einen Vorzeichenwechsel hat, | 
also W(t) = W(t.) + 1 fiir £€ Jt, tp + e[. Aus der Definition von R und S ergibt 


sich damit: W(b) = W(a) — R + S, was zu beweisen war. 


Um dieses Lemma auf unser Problem anwenden zu k6nnen, brauchen wir 


nur noch zu zeigen, dass unsere oben definierte Folge F(t) (22) von Polynomen 
die Eigenschaften I bis IV des Lemmas besitzt. 

Eigenschaft I folgt aus der Bemerkung, dass das letzte Polynom /,(¢) in 
(22) der grésste gemeinsame Teiler von f,(t) und f,(é) ist und f,(¢) und 7,(¢) keine 


- 
) Vel. zu dem Lemma: F. A, WILLERS, Methoden der praktischen Analysis, Kap. IV (Berlin 


1928), Ferner: G. Vatiron, Théorie des fonctions, 2° éd., n° 53 (Paris 1948) 
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Nullstelle auf [a,b] gemeinsam haben. II ist erfiillt, da /,(é) ein Polynom ist, 
III und IV folgen aus (23) und daraus, dass /,(t) und f; wi(4),0 S14 < J, nicht 
gemeinsam verschwinden auf [a, 0]. 


5. Zusammenfassung der Ergebnisse 


Saiz 7. Wir betrachten eine geschlossen orientierte Kurve € der komplexen 
w-Ebene, die nicht die Punkte w= 0 und w= oo enthilt und die sich aus 
y Stiicken ©, von rationalen Kurven zusammensetzt. Die Windungszahl 
n(0, ©) ist gemass (7) gleich der Summe tiber die 7 Integrale J » (8) und damit 
nach den Betrachtungen von Abschnitt 3 und 4 durch 


n(0, ©) = oy Wo(4e) — Wo(be) i. No(4e) — No(be) | (26) 
A 


2 a J 


gegeben (9 = 1,...,7). Die hier auftretenden Zahlen W,(a,), W,(b,), No(4), 
N,(0,) bestimmen sich folgendermassen: 

(5) und (6) sei eine Parameterdarstellung w,(¢) des Kurvenstiickes €,. Zu 
den Zahlerpolynomen P,(t) und Q,(¢) bestimmt man gemass (12) und (13) die 
Polynome /,9(¢) und f,,(¢) (in den Abschnitten 3 und 4 ist der Index @ unter- 
_driickt). Durch (17) ist dann N,(t) und damit speziell N,(a,) und N,(0,) de- 
finiert. Zu /,o(¢) und f,,(¢) erklart man dann mit Hilfe des euklidischen Algorith- 

mus gemass (23), (24) eine Folge 


Fo) = Foolt) » «+» foi,(4)) , = 9) (27) 


von Polynomen, wobei foi, (t) der grésste gemeinsame Teiler von f/,) und f,; 
bzw. P,(t) und Q,(¢) ist. W,(t) bezeichnet dann die Anzahl der Vorzeichen- 
-wechsel von F,(t) in (27) an der Stelle ¢; damit sind dann auch W,(a,) und 
W,(0,) in (26) bestimmt. 

Satz 2 (Eine Anwendung von Satz 1). Es sei J” ein einfach zusammenhan- 
-gendes endliches Gebiet der komplexen z-Ebene. Der Rand ¢ von I’ setze sich 
aus 7 Stiicken ¢, von rationalen Kurven zusammen. Ferner sei w(z) ein Poly- 
nom, das keine Nullstelle auf dem Rand ¢ besitze. Dann erfiillt die Bildkurve 
€ = w(c) von c die Voraussetzungen von Satz 1, und zwar setzt sich © aus den 
y rationalen Kurvenstiicken €, = w(c,) zusammen, die Bilder von den c¢, unter 
der Abbildung z > w(z) sind. Die Anzahl N der Nullstellen von w(z) in "ist 
durch die Windungszahl n(0, €) gegeben, die man nach dem in Satz 1 beschrie- 
_benen Verfahren bestimmen kann. 
Bemerkungen. 1. Wenn man also nach dem in Satz 2 beschriebenen Verfah- 
ren die Anzahl der Nullstellen eines Polynoms w(z) im Gebiet J” bestimmen 
will, so wahlt man zuniichst eine gewisse Parameterdarstellung 


z,(t) ue oO 5 qo(t) (a, oe RE b,) (28) 
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von ¢,; hier sind £,(t), Yo(t), 7o(¢) reelle Polynome ohne gemeinsamen Faktor . 
+ const, und 7,(¢) verschwindet nicht auf [a,, |. Ist dann w(z) das zu unter- 
suchende Polynom, so erhalt man aus 


w,(t) = w(z,()) (29) 


mit (28) eine Parameterdarstellung vom Typ (5), (6) der Bildkurve €, = w(c,). 
Auf diese Darstellungen (5), (6) von ©, wird dann das in Satz 1 beschriebene 
Verfahren zur Bestimmung von n(0,€) angewandt, womit man dann die 
Anzahl N der Nullstellen von w(z) in J” gefunden hat. 

2. In Satz 2 wird vorausgesetzt, dass das Polynom w(z) keine Nullstelle auf 
dem Rand ¢ des Gebiets J" haben soll. Diese Voraussetzung lasst sich im Ver- — 
laufe der Rechnung nachpriifen, und zwar ist sie offenbar damit gleichbedeu- — 
tend, dass keine der Bildkurven €, = w(c,) den Punkt w = 0 enthalt. Anders ~ 
ausgedriickt: In der Parameterdarstellung (5) von ©, diirfen P,(¢) und Q,(t) 
nicht gemeinsam in einem Punkt von [a,, b,] verschwinden. Nun bestimmt man 
aber bei der Aufstellung der Folge F,(¢) in (27) durch den euklidischen Algorith- — 
mus den gréssten gemeinsamen Teiler foi, (4) von P,(t) und Q,(t). Die in Rede 
stehende Voraussetzung ist also gleichbedeutend damit, dass das Polynom 
Ie ip(4) keine Nullstelle auf [a,, 5,] besitzt, und dies kann man ja bekanntlich 


mit Hilfe der gewohnlichen Sturmschen Kette nachpriifen. 


6. Eine Erweiterung von Satz 2 fiir unendliche Gebiete 


Wir erganzen die in Satz 2 beschriebene Methode fiir die Berechnung der — 
Anzahl der Nullstellen in gewissen Gebieten noch dahingehend, dass wir fiir 
diese Gebiete auch den Punkt z = co als Randpunkt zulassen. Auf diese Weise — 
wird zum Beispiel ein durch eine Halbebene oder einen Winkelraum gegebenes — 
Gebiet unserer Methode zuganglich. 

Wir betrachten also ein einfachzusammenhangendes Gebiet J” der kom- 
plexen z-Ebene, das von y rationalen Kurvenstiicken c, begrenzt ist, die ganz 
im Endlichen liegen, mit Ausnahme des Endpunktes von c, und Anfangspunktes 
von ¢,, die mit z = oo zusammenfallen mégen. Wir fragen nach der Anzahl der 
Nullstellen eines Polynoms w(z) in I", wobei keine Nullstelle von w(z) auf dem 
Rand c von I liegen soll. 

Wir vermindern nun zunachst J” um die zu einer Kreisumgebung K von 
z= oo gehérenden Punkte, wobei wir K so wahlen, dass sie keine Nullstellen — 
von w(z) enthalt. Das Restgebiet bezeichnen wir mit I”. I” ist ganz im End- | 
lichen gelegen und, wenn wir nur K hinreichend nahe an co wahlen, einfach- | 
zusammenhdngend. Der Rand ¢’ von J” setzt sich aus 7 + 1 rationalen Kurven- 


stiicken zusammen, Zu denen auch derjenige Kreisbogen f gehért, den der 
Randkreis von K mit J" gemeinsam hat. . 


: 
: 
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Wir kénnen also auf J” die in Satz 2 beschriebene Methode zur Bestimmung 
der Anzahl N der Nullstellen von w(z) in J” anwenden. Diese Anzahl stimmt 
mit der auf das Gebiet J” beziiglichen iiberein. N andert sich nicht, wenn wir 
die Kreisumgebung K sich auf den Punkt z = co zusammenziehen lassen. Dies 
lauft auf folgendes Verfahren fiir die Bestimmung von N hinaus: 

Ergdnzung zu Satz 2. Das einfachzusammenhangende Gebiet J” der kom- 
plexen z-Ebene werde von » Stiicken ¢, von rationalen Kurven begrenzt, die 
ganz im Endlichen liegen mégen mit Ausnahme des Endpunktes von ¢, und 
Anfangspunktes von ¢,, die mit z= co zusammenfallen mégen. w(z) sei ein 
Polynom n-ten Grades, das keine Nullstellen auf dem Rand ¢ von J" besitzt. 
Dann ist die Anzahl N der Nullstellen von w(z) in J” durch 


{ Wo(a@e) — Wo(be) No(4e) — No(%o) | N(%, — ,) 

vm 3 F saemh wenhl ok Saeiae 
gegeben. Hier sind «, und /, mit Hilfe der Parameterdarstellung z,(t), a, St <b 
und 2z,(¢), 4, St S by, von ¢, und ¢, definiert durch 

By = lim argz,(#), a= lim argz,(¢) . (31) 
t—+>b,— t—da,+ 


Die Zahlen W,(a,), W,(bo), No(a), N.(b,) ergeben sich ebenso wie in Satz 2 
beschrieben, soweit es sich nicht um die auf den unendlichen Punkt 


24(0) = 29(4,) = 00 


fiihrenden Parameterwerte 6, und a, handelt. In diesen Fallen definieren wir 
W,(6,) und W,(a,) als Grenzwert durch 


W,(b,) = lim W(t), We(ag) = lim W,() (32) 
t—>b,— t—>a, + 
wahrend wir fiir die Definition von N,(b,) und N,(a,) auf die Formel (16) 
zuriickgreifen und hier die Werte 
argw,(b,) = lim argw,(t), argw,(a,) = lim arg #,() (33) 
t—b,— ay 

einsetzen; hier sind w,(é) = w(z,(¢)) und w,(t) = w(z,(¢)) die Parameterdar- 

stellungen von €, = w(c,) und ©, = w(c,). 
Bemerkung. Im allgemeinen wird man im vorstehenden Fall die Parame- 
terdarstellung z,() von ¢, so wahlen, dass das Parameterintervall die Form 
{a,, co] hat. In diesem Fall bestimmt sich der Wert W,(6,) = W,(co) = lim W,(t) 


einfach als die Anzahl der Vorzeichenwechsel derjenigen Folge, die man erhalt, 
wenn man in der Polynomfolge F,(¢) die Polynome f, ,(¢) durch ihre héchsten 
-Koeffizienten ersetzt. 
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Entsprechend wird man in der vorstehenden Erganzung zu Satz 2 fiir ¢, 
die Parameterdarstellung 2,(¢) so wahlen, dass das Parameterintervall die 
Form [—oo, },] hat. Der Wert W,(a,) = W,(—oo) we lim W,(t) ist dann einfach 


gleich der Anzahl der Zeichenwechsel in derjenigen Folge, die man erhalt, 
wenn man in der Folge F(t) die Polynome f,,(¢) durch ihre héchsten Koeffi- 
zienten, multipliziert mit (—1)%:, ersetzt, wobei g; der Grad von fg; ist. 

Insbesondere kann auch ¢, = ¢, sein, das heisst: J’ wird durch ein einziges 
rationales Kurvenstiick, das durch den Punkt z = co geht, begrenzt. In diesem 
Falle wird man eine Parameterdarstellung z,(¢) von c, so wahlen, dass die im 
Endlichen gelegenen Punkte von ¢, dem Parameterintervall |— oo, cof ent- 
sprechen. So wird man zum Beispiel fiir den Rand ¢ der Halbebene Re(z) < 0, 
also fiir die imaginare Achse, die Parameterdarstellung 2(¢) =7¢ wahlen; ver- 
gleiche auch das folgende Beispiel. 


7. Ein Beispiel 
Wir fragen nach der Anzahl der Nullstellen des Polynoms 
w(z) = 224+ 2—3 _ (34) 
in dem Gebiet J’, das durch die Parabel ¢ = ¢,: 
z,() = #2 —1t -(— co Sf Soo) (35) 
begrenzt wird. | 
Nach dem in Satz 2 und der Erganzung zu Satz 2 beschriebenen Verfahren 


bestimmen wir zundchst die Parameterdarstellung w(t) = w(z,(¢)) der Bild-— 
kurve @, = w(¢): 


w(t) = — 3744+ 2—3471(-38+28—=P(t)+710,0). (36)8 
Wir befinden uns also gemass (12) im Fall 1 und erhalten damit durch Anwen-— 
dung des euklidischen Algorithmus gemass (23), (24) fiir die Folge K(é) von 
Polynomen: 


fo(t) = #8 — 314 + £2 3 

f(t) = 33° ake 4 ¢ | 
fall) = su Be +3 
f(t) = ze ae (37), 
f(t) = 26 t? 23 | 
f(t) = 125 , . 
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Zur Bestimmung der Anzahlen W,(—oo) und iV,(co) der Vorzeichenwechsel 
in der vorstehenden Polynomfolge F,(¢) (36) fiir grosse negative ¢ bzw. grosse 
positive ¢ schreiben wir die Folgen 


fs eye tira eee ged Las ee 

La oe, sgpualagte ited erE : Sian 

oe 5) a, 3” 4 20, 52 — > Solis 3 4 26, 52 re (38) 
der héchsten Koeffizienten der Polynome /,,() auf, einmal multipliziert mit 
(—1)%:, g; = Grad /, ,(¢), das andere Mal ohne diesen Faktor. Hieraus ergibt sich: 


W,(—co) =4,- W,(co) = 2. (39) 
Ferner finden wir 
tat lim argz,(#)=0, £#, =lim argz,()=0. (40) 
—>-—00 t—oo 


Endlich bemerken wir fiir den (auch in unserem Beispiel vorliegenden) Fall, 
dass der Rand ¢ von /’sich aus einem einzigen rationalen Kurvenstiick ¢, mit 
dem Anfangs- und Endpunkt z = co zusammensetzt, dass gilt: 


N,(—oo) — Ny(oo) = 0. (41) 


Denn wenn wir fiir c, eine Parameterdarstellung (28), 9 = 1, mit dem Parame- 
terintervall [—co, co] wahlen, so erhalten wir wegen w,(¢) = w(z,(¢)), w(z) ein 
Polynom: 

lim argw(z,(t)) — lim argw(z,(¢)) = lim arg Ul(— 9) =Ooderz. (42) 
t—>— 00 : t—>0o t—>0o w(z,(t)) 

Hieraus folgt aber mit der Definition (16) von N,(¢) Gleichung (41). 

Fiir die Anzahl N der Nullstellen in unserem Beispiel finden wir nun nach 

Formel (30) mit (39), (41), (40) den Wert: 

West) Wald) 4g. Nal 90) = Naloo) 


S(t Bs). 
5 : oa SALW tas) 


Mer ZG 


Vereinfachungen fiir den Fall der Halbebene. Wir betrachten als Gebiet J" 
speziell eine Halbebene, deren Rand ¢ = ¢, durch 


z(t) =ct+d; —co<t<oo (c,d komplex und c + 0) (44) 
gegeben sei. In diesem Fall vereinfacht sich die Formel (30) fiir die Anzahl NV 


der Nullstellen eines Polynoms w(z) zu 


W,(— 0) — W, (ov) n (45) 


N= 2 2? 


hier sind W,(—co) und W,(co) ebenso wie frither definiert, und 1 ist der Grad 
des Polynoms w(z) (vgl. die Fussnoten 3 und 9). 
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Zum Beweis bemerken wir, dass auf die Halbebene J’ die Voraussetzungen 
zutreffen, die oben zu der Gleichung (41) gefiihrt haben; der mittlere Term in 
(30) fallt also fort. Ausserdem folgt aus der Definition (31) von f, und a, (die 
Indizes 1 und 2 bezeichnen in unserem Fall ja dasselbe!) mit w,(¢) = w(ct+ d): 
w(—ct-+ ad) 


= lim ae era) 


a, — B, = lim arg w(ct + d) eee arg w(ct + a) 


t{—>— 00 


Summary 


We are dealing with the problem of calculating the number N of zeros of a 
polynomial w(z) in certain regions J’ of the complex plane. If I" is a halfplane 
or a disc, there exist many well-known solutions of the problem. If J’is of a more 
general type, there is only a paper of SHERMAN‘) for the case, that J” is the inter- 
section of two halfplanes. In this paper we give the solution for the case that the 
region J’ is simply connected and the boundary of J’ is composed by a finite 
number of pieces of rational curves. Our method generalizes the method of 
WIELANDT®) for the halfplane and the method of SHERMAN‘). The core of our 
paper is the way of computing the variation of arg w, if w runs along a piece of 
a rational curve. This computation is based on a lemma about generalized Sturm 
sequences, 
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Uber eine anschauliche Darstellung der Vorginge bei der 
Fortpflanzung von unstetigen Signalen in Wellenleitern 


Von ADALBERT RUBINOWICcz, Warschau!) 


1. Wellenleiter mit beliebigem Querschnitt 


Zur Darstellung der Vorgange in elektromagnetischen oder akustischen 
Wellenleitern im Falle der Fortpflanzung von abgehackten, zeitlich periodi- 
schen Wellensignalen ist es sehr bequem, die gesamte Wellenbewegung in. 
einzelne Wellentypen aufzulésen, von denen sich jede mit einer bestimmten 
Phasen- und Gruppengeschwindigkeit ausbreitet. Solche Wellentypen ergeben 
sich, falls man die gesamte Wellenbewegung, ebenso wie im stationaren Falle, 
nach den Eigenschwingungen des Querschnitts des Wellenleiters entwickelt 
[5]?). Insbesondere kann man durch Superposition von solchen Wellentypen die 
Vorgange darstellen, die in einem Wellenleiter im akustischen Falle von eine 1 
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punktformigen Strahler nullter Ordnung und im elektromagnetischen Falle von 
einer punktférmigen Dipolquelle angeregt werden [2]. Speziell fiir die Anwen- 
dung in der Radartechnik ist eine solche Entwickiung von Bedeutung, weil hier 
die Wellenlange der Dipolstrahlung und der Querschnitt des Wellenleiters in 
der Regel aneinander derart angepasst.werden, dass im stationiren Falle in 
weiterer Entfernung von der Strahlungsquelle nur ein einziger Wellentypus 
in Frage kommt. 

Eine solche Zerlegung des Wellenvorganges hat aber den schwerwiegenden 
‘Nachteil, dass sie uns tiber den Entstehungsmechanismus der Vorlaufer und 
ihre Eigenschaften gar keine Auskunft gibt. Um einen Einblick in das Wellen- 
feld der Vorlaufer zu erhalten, muss man namlich die Felder der Vorlaufer aller 
Wellentypen superponieren, ohne Riicksicht darauf, ob fiir sie der Wellenleiter 
durchlassig oder undurchlassig ist. Dies fiihrt notwendigerweise auf einen ganz 
uniibersichtlichen Ausdruck fiir die gesamte Wellenbewegung. 

Um uns eine Vorstellung von der Entstehung der Vorlaufer zu verschaffen, 
-miissen wir daher den Wellenvorgang als ein einheitliches Ganzes betrachten. 
Wir wollen dabei unseren Uberlegungen den praktisch wichtigeren Fall der 
elektromagnetischen Wellen in einem von vollkommen leitenden Wanden 
begrenzten Wellenleiter zugrunde legen. Sie gelten aber im grossen und ganzen 
auch in dem einfacheren akustischen Falle. 

Wir setzen voraus, dass die elektromagnetische Strahlung durch einen 
infinitesimalen elektrischen Dipol angeregt wird, der etwa zur Zeit t = 0 zu 
schwingen beginnt. Der Beginn der zeitlichen Abhangigkeit seines Dipol- 
momentes soll dabei nur so beschaffen sein, dass die Wellenfront der Strahlung 
bei der freien Ausbreitung den Bedingungen geniigt, die an einer sich mit 
Lichtgeschwindigkeit bewegenden elektromagnetischen Sprungflache erfiillt 
sein miissen (vgl. [6] und [7]}). Der weitere Verlauf der Dipolschwingungen soll 
zunachst an keine Voraussetzungen gebunden sein. 

Bis die Strahlung an die Wande des Wellenleiters gelangt, pflanzt sie sich 
zunichst als freie Strahlung fort, so als ob der Wellenleiter gar nicht vorhanden 
ware. Dies ist intuitiv klar, kann aber auch ganz streng mit Hilfe des Eindeu- 
tigkeitssatzes fiir das gemischte Problem des elektromagnetischen Feldes [4], [6] 
bewiesen werden. Eine solche freie Strahlung besitzt ja die fiir unser Problem 
vorgeschriebene Singularitat, geniigt den Anfangsbedingungen und erfiillt bis 
zu ihrem Auftreffen an die Wande des Wellenleiters auch die Randbedingungen. 

Weiter kann man behaupten, dass auch spaterhin in der ganzen Wellenfront 
selbst und knapp dahinter die elektromagnetische Strahlung sich von der 
freien Strahlung des Dipols nicht unterscheidet. Der Grund dafiir liegt darin, 
dass die Strahlung in der Wellenfront nur durch das von den Wanden des 
‘Wellenleiters reflektierte und abgebeugte elektromagnetische Feld abgeandert 
werden kann. Der Weg, den diese letztere Strahlung zuriicklegen miisste, um 
zur Wellenfront zu gelangen, ist jedoch linger als der direkte Weg vom strah- 
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lenden Dipol. Auch diese Tatsache lisst sich mit Hilfe des erwahnten Eindeu- 
tigkeitssatzes streng beweisen [7]. Wir sehen somit, dass selbst im Falle eines 
punktférmigen Dipols die Wellenfront gar nicht eben ist, sondern durch eine 
Kugelflache gegeben wird. Je weiter jedoch diese Wellenfront im Wellenleiter 
fortschreitet, um so mehr niahert sie sich einer Ebene, die senkrecht auf der 
Fortschreitungsrichtung steht. Dabei nehmen in der Wellenfront die Feld- 
stirken ebenso wie bei der freien Ausbreitung ab. Hier macht sich also noch 
nicht die energiesammelnde Wirkung des Wellenleiters bemerkbar. 

Die Wellenfront ist eine elektromagnetische Sprungflache, da die Feld- 
starken an ihrer Vorderseite verschwinden und nur an ihrer Riickseite im all- 
gemeinen von Null verschieden sind. Eine solche Sprungflache pflanzt sich 
aber bekanntlich (vgl. etwa [6], [7]) wie eine Huygenssche Hiillflache fort. Ihr 
Fortschreiten kann daher mit Hilfe des Huygensschen Enveloppenprinzips 
verfolgt werden. Bei der Reflexion und Beugung einer solchen Sprungflache 
entsteht wieder eine Sprungflache. Die Lage der reflektierten Wellenfront kann 
daher ebenfalls mit Hilfe des Huygensschen Enveloppenprinzips festgelegt 
werden. Jeder Punkt der, den Wellenleiter begrenzenden, Wand bildet dabei 
den Ausgangspunkt von Huygensschen Sekundarwellen. Die Lage der reflek- 
tierten Hiillflache wird daher durch die Hiillflache aller Kugeln gegeben, deren 
Mittelpunkte auf der reflektierenden Wand legen, und die auch die urspriing- 
liche einhiillende Hiillflache nach ihrem unbehinderten Durchgange durch die 
Wand einhiillen. Daraus folgt, dass die einfallende und die reflektierte Hiill- 
flache sich in einer auf der reflektierenden Wand liegenden Kurve schneiden 
und dass die Fortschreitungsrichtungen der beiden Einhtillenden an den Wan- 
den des Wellenleiters mit dessen Langsrichtung die gleichen Winkel ein- 
schliessen. 

Die durch die Reflexion einer Wellenfront entstehende Sprungflache stellt 
wieder eine Wellenfront dar, da auf ihrer Vorderseite die elektromagnetischen 
Feldstarken ebenso verschwinden wie bei der einfallenden Wellenfront, aus der 
sie entsteht. Durch Reflexion der zweiten Wellenfront entsteht eine dritte usw. j 
Auf diese Weise wird das im Wellenleiter befindliche elektromagnetische Feld 
in jedem Zeitmoment ¢ > 0 von einer endlichen, im Laufe der Zeit aber un- 
unterbrochen anwachsenden Anzahl von elektromagnetischen Sprungflachen’ 
durchsetzt. Alle die angefiihrten Tatsachen lassen sich mit Hilfe des Eindeutig- 
keitssatzes streng beweisen (vgl. [7], Kap. VI). In allen diesen elektromagne- 
tischen Sprungflachen miissen selbstverstandlich die hier geltenden Bedingun-. 
gen (vgl. zum Beispiel (6]) erfiillt sein. Jede dieser Sprungflachen signalisiert 
die Ankunft eines neuen Wellenfeldes, das durch Reflexion und Beugung aus 
dem gerade ihm vorangehenden entsteht. Da die bereits vorhandene Wellen- 
front sich mit der aus ihr entstehenden neuen in einer langs der Wand verlau- 


fenden Kurve schneidet, hangen alle die aufeinanderfolgenden Wellenfronten 
miteinander zusammen. 
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Bezeichnet # den Winkel zwischen der Normalen in einem gegebenen Punkte 
der Wellenfront und der Normalen an den Querschnitt des Wellenleiters, so 
bewegt sich die Wellenfront langs einer zur Langsrichtung des Wellenleiters 
parallelen Geraden mit der Geschwindigkeit c/cos . Man kann dies leicht fest- 
stellen, wenn man sich in den beiden Zeitmomenten ¢ und ¢ + dt die Lagen der 
Wellenfronten vergegenwartigt und beachtet, dass sich diese mit der Licht- 
geschwindigkeit c in der Richtung ihrer Normalen bewegen. Je grésser somit 
der Winkel # ist, um so schneller schreitet die Wellenfront in der Langsrichtung 
des Wellenleiters fort. Nehmen wir an, dass die Vorderseite der Wellenfront in 
jedem ihrer Punkte nach aussen konvex ist, das heisst, dass der Winkel  mono- 
ton wachst, wenn wir ausgehend von dem vordersten Punkte der Wellenfront 
zu ihren riickwarts gelegenen Punkten fortschreiten. Es werden sich dann die 
weiter zuriickliegenden Partien der Wellenfront schneller bewegen als die vor- 
deren, und die Wellenfront wird daher bei ihrem Fortschreiten sich immer mehr 
einer senkrecht zur Langsrichtung des Wellenleiters liegenden Ebene nahern. 
Ist hingegen die Vorderseite der Wellenfront nach aussen konkay, so wird sie 
sich bei ihrem Fortschreiten immer schiefer gegentiber der Langsrichtung ein- 
stellen. Wie man sich aus der oben angegebenen Konstruktion der reflektierten 
Hiillflachen klarmachen kann, entsteht aus einer Hiillflache, die an ihrer 
Vorderseite nach aussen konvex ist, eine ebensolche reflektierte Hiillflache. Da 
im Falle eines Wellenleiters die vom Dipol bei der freien Ausstrahlung emittierte 
Wellenfront nach aussen konvex ist, so miissen auch alle aus ihr durch sukzes- 
sive Reflexionen und Beugungen entstehenden Wellenfronten an ihrer Vorder- 
seite nach aussen konvex sein. Sie werden sich daher bei ihrem Fortschreiten in 
der Langsrichtung des Wellenleiters immer mehr Ebenen nahern, die senkrecht 
auf der Langsrichtung des Wellenleiters stehen. Da alle sukzessiven Wellen- 
fronten langs der Wande des Wellenleiters miteinander zusammenhangen, 
bedeutet dies, dass sie sich auch einander immer mehr nahern miissen. Die 
sukzessiven Wellenfronten werden sich daher immer mehr an der vordersten 
Wellenfront zusammendrangen. 

Es entsteht nun die Frage nach der Intensitatsverteilung in den auf diese 
Weise entstehenden Vorliufern. Um streng giiltige Aussagen machen zu kénnen, 
miisste man die Intensitatsverteilung in den durch wiederholte Reflexion und 
Beugung an den Wanden des Wellenleiters entstehenden Wellen zur Verfiigung 
haben, die einzig in dem im nachsten Abschnitt zu besprechenden Spezialfall 
(wo nur Reflexionen an Ebenen stattfinden) bekannt ist. 

- Um die Entstehung der den verschiedenen Wellentypen entsprechenden 
Hauptwellen plausibel zu machen, betrachten wir zunachst die Gesamtheit der 
iiber einen Beobachtungspunkt P hinwegstreichenden Wellenfronten. Da sich 
diese in der Nahe der vordersten Wellenfront haufen, werden sie in dem Punkte 
P nach dem Eintreffen der vordersten Wellenfront zunachst in sehr rascher, 
mit der Zeit aber immer mehr abnehmender Aufeinanderfolge eintreffen, Eine 
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Hauptwelle mit einer gegebenen Frequenz kann sich in einem Punkte P erst 
dann bemerkbar machen, wenn die Frequenz der Aufeinanderfolge zweier Wel- 
lenfronten mit sich verstirkenden Feldstarkeverteilungen einer Eigenfrequenz 
des Querschnittes des Wellenleiters gleich wird. Um einzusehen, dass eine solche 
Hauptwelle auch spaterhin bestehen bleibt, muss man noch annehmen, dass die 
Phasenflichen der nach den Wellenfronten ankommenden Teile der reflektierten 
und gebeugten Wellen im Takte der Wellenfronten den Punkt P passieren. Wir 
wollen jedoch hier auf dieses Problem nicht weiter eingehen. Im Abschnitt 2 
wird sich dieser Tatbestand im Falle eines rechteckigen Querschnittes auf 
einem ganz elementaren Wege ergeben. 

Offenbar kann die Zeit, die zwischen der Aufeinanderfolge zweier Wellen- 
fronten verstreicht, nicht langer sein als die Periode der tiefsten Eigenschwin- 
gung des Querschnittes. Auf diese Weise versteht man, dass der Wellenleiter 
eine Sperrwirkung aufweisen muss fiir die Frequenzen, die kleiner sind als die, 
die dieser Schwingung entspricht. 


2. Die Wellenbewegung in den Vorlaufern in Wellenleitern mit Quer- 
schnitten, die sich mit Hilfe des Thomsonschen Spiegelungsverfahrens 
behandeln lassen 


Um uns ein vollstandigeres Bild von den Vorgangen in Wellenleitern zu ver- 
schaffen, wollen wir den Spezialfall in Betracht ziehen, bei dem die Lésung der 
Maxwellschen Gleichungen im stationaren Falle mit Hilfe des Thomsonschen 
Spiegelungsverfahrens hergestellt werden kann. Der Querschnitt eines solchen 
Wellenleiters wird durch ein Dreieck oder Rechteck gegeben, und der urspriing- 
liche Dipol und die Gesamtheit seiner Bilder liegen in einer durch einen senk- 
rechten Querschnitt hindurchgehenden Ebene in den Punkten eines zweidimen- 
sionalen Gitters. Figur 1 zeigt ein solches Gitter fiir den Fall eines rechteckigen 
Querschnittes, wobei der Einfachheit halber angenommen wird, dass das Dipol- 
moment in der Querschnittsebene liegt. 

Um fiir einen Wellenleiter, der sich nach dieser Methode behandeln lasst, die 
Lésung im Falle eines abgehackten Wellensignals zu erhalten, muss man nur 
den urspriinglichen Dipol sowie seine sémtlichen Spiegelbilder synchron und 
nach ein und demselben, sonst aber beliebigen Zeitgesetze, etwa vom Zeit- 
moment ¢ = 0 an, schwingen lassen. Zu einer gegebenen Zeit ¢ > 0 besteht dann 
die Wellenbewegung in dem Wellenleiter aus der Welle des urspriinglichen 
Dipols sowie den Wellen aller seiner Bilder, von denen die Wellenbewegung 
in dem Zeitabschnitt (0, ¢) zu den Raumpunkten des Wellenleiters gelangen 
konnte. In jedem Raumpunkte P des Wellenleiters setzt sich somit das elektro- 
magnetische Feld aus der Wellenbewegung zusammen, die von dem urspring- 
lichen Dipol sowie allen seinen Bildern ausgestrahlt wurde, die sich in einer 
Kugel befinden, deren Mittelpunkt in P liegt und deren Radius gleich c¢ ist 
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fe = Lichtgeschwindigkeit). Die Wellenbewegung in P wird somit stets durch 
die Superposition einer endlichen Anzahl von Kugelwellen erhalten. Man kann 
also in dem betrachteten Spezialfalle ganz einfach und elementar den sukzes- 
siven Aufbau der Wellenbewegung verfolgen. Dies wird durch die Tatsache 
bedingt, dass hier (wie aus der Sommerfeldschen Theorie der Beugung an einem 
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Keil oder Winkelspiegel bekannt ist) nur Reflexionen, aber keine Beugungser- 
scheinungen auftreten. 

Die angegebene Entstehungsweise der Wellenbewegung bestatigt zunachst 
im betrachteten Spezialfall das Ergebnis unserer allgemeinen Uberlegungen in 
Abschnitt 1. Die von den einzelnen Dipolbildern emittierten Wellen entsprechen 
den an den Wanden des Wellenleiters reflektierten Strahlungen. In Figur 1 
wurden in jedem zu einem Dipolbild gehérigen Rechteck die Wande Aj, A», 
B,, B, eingetragen, an denen die urspriingliche Strahlung reflektiert werden 
musste, um die von dem betreffenden Bildpunkte ausgehende Strahlung zu 
ergeben. Der Reihenfolge der Wande ist dabei keine besondere Bedeutung bei- 
zumessen. Die Strahlung, die von dem rechts oben gelegenen, mit A, A, B, 
bezeichneten Dipolbild ausgeht, kénnen wir uns durch Reflexion an der Ebene 
A,, sodann an der Ebene B, (es ergibt sich auf diese Weise eine Strahlung, die 
dem Dipolbilde A, B, ulrspticht) und schliesslich an der Ebene A, entstanden 
denken. Wir penne uns aber auch vorstellen, dass diese Welle hist sukzessive 
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Reflexionen an den Ebenen B,, A, und A, zustande kommt. Durch eine 
Reflexion an den gleichen Ebenen (in der Reihenfolge A,, A, B, oder By, Ay, 
A,) entsteht aber auch eine Strahlung, die wir in der Figur 1 dem links 
oben ebenfalls mit A, A, By bezeichneten Dipolbilde zuordnen. Unter der An- 
nahme, dass der strahlende Dipol in einer Querschnittsebene liegt, sind in der 
Figur 1 die Lagen der Dipolbilder durch entsprechende Pfeile angedeutet. 
Man erkennt, dass das ganze Dipolgitter in vier ineinandergestellte einfache 
Dipolgitter zerfallt, von denen jedes einer bestimmten, durch Pfeile angedeu- 
teten Dipolrichtung entspricht. Der urspriingliche Dipol liegt in dem Felde, 
in dem keine Buchstaben eingetragen sind. Die Langen der Seiten des Periodizi- 
titsrechteckes werden durch die doppelten Seitenlangen 2a und 20 des Recht- 
eckes gegeben, das den Querschnitt des Wellenleiters darstellt. 

Die Wellenfronten der von den einzelnen Dipolbildern ausgehenden Kugel- 
wellen stellen die Wellenfronten der einzelnen, an den Wanden des Wellenlei- 
ters reflektierten Wellen dar. Da zur Zeit t diese Wellenfronten auf Kugeln mit 
den Radien ct liegen, deren Mittelpunkte sich in den einzelnen Dipolbildern 
befinden, bestatigt es sich, dass die vorderste Wellenfront, die vom urspriing- 
lichen Dipol stammt, kugelf6rmig ist und dass das elektromagnetische Feld 
auf ihr und knapp hinter ihr durch die freie Strahlung dieses Dipols gegeben 
wird. Auch alle nachfolgenden Wellenfronten sind Kugelflachen. Dass alle 
diese Wellenfronten in ihren Schnittkurven mit den Wanden des Wellenleiters 
miteinander zusammenhangen, folgt aus ganz einfachen elementargeometrischen 
Betrachtungen. 

Einen Begriff von der gegenseitigen Lage der von den einzelnen Dipolbildern 
stammenden Wellenfronten vermittelt Figur 2. Sie stellt die Wellenfronten 
einer linearen Dipolreihe (etwa ...A?A,,.A,A,, 4,, urspriinglicher Dipol, 
Ay, A, Ay, A, AG, ...) in einer durch sie hindurchgehenden Ebene dar. Die 
Wellenfronten tragen die aus der Figur 2 zu entnehmenden Nummern der 
Dipolbilder, von denen sie ausgestrahlt wurden. Man sieht, dass alle sukzessiven 
Wellenfronten miteinander zusammenhdngen und dass sie sich an der vorder- 
sten Wellenfront zusammendrangen. Es ist nicht schwer, sich klarzumachen, 
dass, je weiter die vorderste Wellenfront vorriickt, eine stets wachsende Anzahl 
nachfolgender Wellenfronten sich immer enger an sie anschliesst. : 

Um sich in dem betrachteten Spezialfall einen quantitativen Einblick in die 
Beschaffenheit des gesamten Feldes zu verschaffen, muss man nur von einer 
entsprechenden Darstellung der im rein periodischen Fall herrschenden Ver- 
haltnisse ausgehen. In diesem Fall kann man namlich das in einem Beobach- 
tungspunkte P des Wellenleiters auftretende Feld im wesentlichen mit Hilfe 
der Methode der stationaéren Phase berechnen. Dabei muss jedoch vorausgesetzt 
werden, dass die Entfernung des Punktes P von der Gitterebene im Vergleich 
mit den Seitenléngen aund b des rechteckigen Querschnittes hinreichend gTOss | 
ist. Man hat die Rechnung zunichst fiir jedes einzelne der vier einfachen Dipolg 
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gitter durchzufiihren, aus denen das gesamte in F igur 1 dargestellte Dipol- 
gitter besteht [8], [1]. Das Ergebnis der Rechnung kann man dahin zusammen- 
fassen, dass die in P eintreffende Strahlung in erster Naherung von gewissen 
um «wirksame Punkte» gelegenen « wirksamen Bereichen» zu kommen scheint, 
von denen jeder nur eine relativ kleine Gruppe von Dipolbildern umfasst. Vom 
Punkte P aus gesehen, liegen die wirksamen Punkte in Richtungen, in denen 
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Die Wellenfronten der Wellen des urspriinglichen Dipols und seiner Bilder. 


das ebene Dipolgitter gemass der Fraunhoferschen Beugungstheorie die Spek- 
tren der verschiedenen Ordnungen ausstrahlt. Dies findet statt, trotzdem wir 
es hier nicht mit der Fraunhoferschen, sondern mit der Fresnelschen Beugung 
gu tun haben. Wie aus der Figur 1 zu ersehen ist, haben von den vier in- 
einandergestellten einfachen Dipolgittern je zwei antiparallele Dipolmomente. 
Dies bedeutet, dass in erster Naherung das Spektrum nullter Ordnung, das 
senkrecht zur Dipolebene ausgestrahlt wird, nicht auftritt. Dieser Umstand hat 
aber nicht zur Folge, dass auch die Spektren anderer Ordnungen verschwinden. 
Bei schiefer Ausstrahlung treten namlich zwischen den Strahlungen der vier 
ineinandergestellten einfachen Dipolgitter, da sie in der Gitterebene gegen- 
‘einander verschoben sind, noch zusatzliche Phasendifferenzen auf. Dass das 
Spektrum nullter Ordnung verschwindet, wird sozusagen durch das Ver- 
-schwinden des Strukturfaktors bedingt. 
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Die Wellenfelder, die von den einzelnen wirksamen Punkten ausgestrahlt 
werden, kénnen in der ersten, das heisst Fraunhoferschen Naherung als eine 
Superposition eines elektrischen und magnetischen Wellentypus von der glei- 
chen Ordnung angesehen werden. Im Falle eines rechteckigen Querschnitts 
fallen niimlich die Phasengeschwindigkeiten dieser beiden Wellentypen im all- 
gemeinen miteinander zusammen. . 

Aus der obigen Darstellung des Mechanismus des Zustandekommens der 
Wellenbewegung im periodischen Falle kann man bei abgehackten pericdischen 
Signalen unmittelbar einige Folgerungen ziehen. Dabei soll angenommen 
werden, dass die vom urspriinglichen Dipol ausgesandte Strahlung nach einer 
sehr kurzen Anlaufzeit (die notwendig ist, damit die Strahlung an der Wellen- 
front die Bedingungen fiir die bewegten elektromagnetischen Sprungflachen 


erfiillt) praktisch periodisch wird. Ist der Radius c ¢ der Kugel um P so klein, — 


dass in dem Kreise, die diese Kugel aus der Dipolebene herausschneidet, nur 


der urspriingliche Dipol liegt, so besteht im Punkte P die gleiche elektromagne- 


tische Strahlung wie bei der freien Ausbreitung. Sobald aber dieser Schnitt- 


kreis bereits so gross geworden ist, dass er den ganzen wirksamen Bereich um — 


den urspriinglichen Dipol enthalt, so verschwindet in P die Strahlung in erster 
Naherung, da die von den vier verschiedenen einfachen Dipolgittern stammen- 
den Wellen einander durch Interferenz zerstéren®). 


Die Intensitaét der Strahlung steigt erst wieder an, sobald der Schnittkreis 


einen der iibrigen wirksamen Bereiche erreicht. Dieser Zeitpunkt markiert den 
Beginn der Ankunft der Hauptwelle des elektrischen und magnetischen Wellen- 


typus, denen der betreffende wirksame Punkt entspricht. Man kann dabei er- 
warten, dass dieses Anwachsen durch ein Fresnelsches Integral gegeben wird. — 


Wie aus der Kirchhoffschen Theorie der Fresnelschen Beugungserscheinungen 
bekannt ist (vgl. [7]), ergibt sich namlich immer ein solches Integral, sobald nur 


uber ein Teilgebiet des wirksamen Bereiches integriert wird. Dass dies tat-— 


sdchlich zutrifft, zeigt eine von PEARSON [3] im akustischen Falle abgeleitete, 
hier aber unverandert geltende Naherungsformel!). 

Um die Fortpflanzungsgeschwindigkeit der von einem ins Auge gefassten 
wirksamen Punkt kommenden Hauptwelle festzulegen, bezeichnen wir mit @ 
den Winkel, den im Beobachtungspunkte P die Blickrichtung nach dem ent- 
sprechenden wirksamen Punkte mit der Normalen an die Gitterebene ein- 


3) Auch in dem von SomMERFELD [9] behandelten Falle eines in einem dispergierenden Medium 
sich fortpflanzenden abgehackten Wellensignals ist die Intensitat der Feldstarken in den Vorlaufern 
dusserst gering. Anders liegen die Verhaltnisse im Falle eines akustischen Wellenleiters mit schall- 
harten Wanden. Wegen der Grenzbedingung Oy/On = 0 fiir das Geschwindigkeitspotential gp haben 
die vier ineinandergestellten einfachen Schallquellengitter alle die gleichen Phasen. Hier findet also 
kein Absinken der Wellenamplitude in den Vorliufern statt. 

4) Im akustischen Falle sind die Beitrage der iibrigen wirksamen Bereiche kleiner als die des 
zentralen, in dem die urspriingliche Schallquelle liegt. Infolge der zwischen den vier ineinander- 


gestellten einfachen Gittern auftretenden Phasendifferenzen wird ja die Gesamtstrahlung eines 
solchen wirksamen Bereiches geschwacht. 


= 
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schliesst. Sei ¢,, der Zeitpunkt, in dem die Hauptwelle in dem Punkte P ein- 
trifft. Die Entfernung des Punktes P von der Gitterebene wird dann durch 
z= Ct, cos® gegeben. Die Geschwindigkeit der Fortpflanzung der Haupt- 
welle wird somit gleich vg = 2/ty = ¢ cos@. 

Der Winkel # wird durch die gleiche Bedingung festgelegt, die im Falle der 
Fraunhoferschen Beugungserscheinungen fiir ein ebenes rechteckiges Gitter die 
Richtung der Fortpflanzung der Spektren der verschiedenen Ordnungen be- 
stimmt. Fiir ein Spektrum, dessen Ordnung durch die ganzen Zahlen m und 
n festgelegt wird, wird cos ?,,,, durch den Ausdruck 


cost... = 1 (52) ~ (Fy (m,n = 0, +1, +2, :..) 


gegeben. Fiir die Geschwindigkeit vg ergibt sich somit der bekannte Ausdruck 
fiir die Gruppengeschwindigkeit in einem Wellenleiter vom rechteckigen Quer- 


schnitt oe = 
Whe > m 2 dai (nr 2 
ug =0|/1 (>= (Fe) 
Auch der Ausdruck fiir die Phasengeschwindigkeit vp; = c/cos#,,, lasst sich 
aus dem verwendeten Bilde einfach ermitteln. 
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Summary 


The electromagnetic radiation in a wave guide due to an electric dipol may 
be considered as consisting of the incident primary radiation and the waves 
generated by successive reflections from the walls of the pipe. This picture was 
applied to study the propagation of a cut-off periodic wave, first generally for 
any section of the wave guide and then more in detail for a rectangular section. 
The field consists in any given moment after the beginning of the dipol radiation 
of a finite number of such waves, the fronts of which are generally electromagnetic 
- wave surfaces. In the foremost wave front and the points just behind it the 
electromagnetic field is given by the primary radiation. But very soon the field 
subsides in the forerunner till the arrival of the first main wave. 


(Eingegangen: 21. Juli 1955.) 
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The Rolling-Up of a Vortex Sheet 


By Marvin STERN, San Diego, California’) 


1. Introduction 


A surface of discontinuity or vortex sheet is known to exist behind a wing 
of finite span [1]2). Further, this sheet, which is shed planar, is known to roll 
in from its tips. Tests have indicated a behavior as shown below. 

This, and related problems such as the rolled vortex sheet occurring in flow 
around a sharp corner, were first discussed by PRANDTL. In lectures published 
in Berlin in 1924, he presented some of this work where he had represented the 
rolled sheet in the form of a logarithmic spiral. 

In 1931, KapEN [2] investigated this problem both experimentally and 
theoretically. As an approximation in the theoretical treatment, a two-dimen- 
sional time dependent flow was considered rather than the three-dimensional 
steady state problem that corresponds to Figure 1. ; 


——*< +4 


Figure 1 Figure 2 
Vortex sheet behind wing of finite span. Time dependent rolling of vortex sheet. ~ 


A Ay segment of the vortex sheet is chosen. A coordinate system is fixed 
to this segment while the wing is allowed to move relative to it with a velocity 
—v. The influence of the entire sheet on this segment is assumed to differ little 
froma cylindrical extension of it. The resulting two-dimensional time dependent 
flow is then considered to respresent the real flow with the correspondence y= vt. 


z Convair. A Division of General Dynamics Corporation. 
) Numbers in brackets refer to References, page 342. 
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A further simplification is effected by assuming that the distance between 
vortex cores in Figure 2 is great enough so that their influence upon each 
other is small. Half of the sheet can then be studied alone, thus finally reducing 
the problem'to the analogous one of the semi-infinite vortex sheet. 


= 


Se 


Figure 3 


Rolling of semi-infinite vortex sheet. 


In this latter problem, KADEN recognizes the fact that no fundamental 
length appears. It is concluded, then, that the resulting motion of the vortex 
sheet must consist merely of homogeneous expansions, thus continually yield- 
ing geometrically similar figures with time. 

In 1935, WESTWATER |{3] approximated the continuous vorticity distribution 
along the sheet by a set of discrete vortices and performed a step-by-step 
numerical calculation giving the form of the sheet at various increments of 
time. His results showed fairly good agreement with KADEN’s. 

We will consider the problem of the two-dimensional time dependent flow 
associated with the semi-infinite vortex sheet. Although KADEN obtains an 
approximation to the form of the spiral curve, it would be difficult to obtain 
the associated flow field from KADEN’s results. We will work with the complex 
potential of the flow, thus yielding information on the entire flow field. 

In contrast with KADEN, we obtain differential equations describing the 
flow field that holds over the whole plane, rather than just conditions that are 
to be satisfied locally at the spiral core. We thus have means of determining 
series expansions for our solutions, whereas KADEN can properly only give the 
leading term in an asymptotic behavior at the spiral core. 

At the spiral core, our solution agrees with that of KADEN. We find the 
nature of the essential singularity there. In addition, we find several more terms 
and then show means for obtaining successive terms uniquely. 

Although we were unable to develop an expansion at infinity that satisfied 
our differential equations, we did show consistency of our assumed asymptotic 
behavior there. 


2. Formulation of the Problem 


We consider the problem of the rolling up of a semi-infinite vortex sheet. 
Specifically, we take the two-dimensional potential flow about a semi-infinite 


plate, W = 2 a(z)'”. 


es 
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Consider the plate to be instantaneously withdrawn, leaving the semi-infinite 
vortex sheet with vorticity distribution equal to the velocity jump that existed 
across the plate y(x) = 2 a/(x)'”. The flow about this vortex sheet is identical 
now with the flow that took place about the semi-infinite plate. 


| 
| 
| 
| 


Figure 4 
Flow around semi-infinite plate. 


Since the velocity at any point adjacent to the sheet directly above is the © 
exact negative of that directly below, there is no unbalanced pressure across the ~ 
sheet and therefore no motion of the sheet might be expected. This condition 


is violated, however, at the origin. Here the infinite velocity around the origin 
results in a force in the horizontal direction at this corner (the so-called leading 
edge suction). Due to the inability of the sheet to resist any force, the static 


equilibrium is expected to be upset, and a motion or curling of the sheet emanat- 


ing from this free edge is expected. 


Our problem, then, is the two-dimensional time dependent motion of the | 


curling of this sheet and the determination of its associated flow field. 


3 


a 
Figure 5 


Rolled sheet at time#=1. |§,| =1, 
o = arc length (is time independent) 


Figure 6 
Similarity transformation in time. 


The vortex sheet is represented 
by x = {%,,%.} wherex=6&. 
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Modeling after KADEN, we limit ourselves to the study of solutions to this 
problem that are homogeneous in time. We assume one time independent 
curve, §(c), representing the curled vortex sheet at time ¢ = 1. 

In the course of time, the curve shown in Figure 5 representing the curled 
vortex sheet is assumed to magnify according to its homogeneous dependence 
upon time, giving geometrically similar curves generated from rays through 
the origin. Thus, on the (x, ¥,)-plane, the time dependent curves would appear 
as shown in Figure 6. ; 


Conditions to Be Satisfied 


On the (*,, %,)-plane, the hydrodynamic equations of conservation of mass 
and irrotationality of the flow field will be satisfied by finding an analytic 
function that characterizes the complex potential, y, of the flow field. The bound- 
ary conditions to be satisfied are such that the flow at any time adjacent to the 
vortex sheet must be tangential to it. In addition, it must be required that there 
be no pressure differential across the vortex sheet. These conditions will be 
satisfied on the (*,, x,)-plane; the homogeneity in time will lead to equations 
to be satisfied on the (&,, &,)-plane. 


Potential Flow 


We assume that there exists a time dependent velocity potential, g(x), 
which is the real part of an analytic function, such that the velocities in the 
flow field are given by the gradient u=V,q. Associated with this p(x) there 
will exist a time independent velocity potential, @(€), also the real part of an 
analytic function. These potentials are related homogeneously in time, 


g(x) = 1 P(E) . 


Boundary Condition 


The boundary condition to be satisfied is the condition of sliding on the 
(x1, ¥2)-plane. For all time ¢, points adjacent to or on the vortex sheet can have 
velocity parallel to the curve only. Since points of the geometric curve have 
velocity x, due to the homogeneous expansion itself, the sliding condition on 
the resulting velocity is 


u—x,||§. 
for all ¢;i. e., there is a scalar F such that for alli 
Vip — X;= — F(a, t) §,- 
Since x, = 2 ¢"-1 &, we attain homogeneity in ¢ by requiring that 


F(o,t) =i" Fo), Vap= Ve. 
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To satisfy the latter equation, we require 
g(x) = Pr * O(5) . 


Dividing through by the homogeneity factor, ¢"~1, the sliding condition is 
written on the (&,, €,)-plane as 


VD = n& — F(o) &,. (1) 


The vortex sheet represents a curve across which the potential is discon- 
tinuous. If we define @+ as the value of ® on the concave side of the curve, and 
@- on the convex side, then the mean and jump of potential across the curve 
can be written as 


mean Cee (++ @-), jump f-G+-O@. 


a 


Equation (1) can then be written for @ and @® 


Vb=nE—FE,, <a 


Vb= —Fé. @ 


Pressure Condition 
The time dependent Bernoulli equation © 


1 
Wess: (Vo)? + 2 = const 
can be written on either side of the time dependent discontinuity curve. 
Since pressure is continuous across this curve, we find 


G:+VeVp=0.- 


Transcribing to the time independent plane, we find, to within a homogeneity 
factor of 12"—-2, 

(2n—-1) 64 (-nE+76)pé=0.. ; 
q 


Insertion of (1) yields ; 
(2n—-1)6-FO@ =0. 


~ 


3. Determination of the Degree of Homogeneity 


The degree of homogeneity in time is obtained by satisfaction of th 
prescribed asymptotic behavior at infinity on both the x- and &-planes. 


c 
: 
: 
e 
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Behavior for Large ¢ 
On the §-plane we require that the curve becomes asymptotic to the §,-axis. 
§~ { 0,0 } ‘ 


For the potential, we require that its behavior tends to that of flow around the 
slit ; i.e., we require asymptotically, as ¢ > co 


P~+2a(o)"?, Sw4a(o)?. 
The asymptotic velocity behavior would then be, as required, 
VO ~ +a(o)-*"£1,0}, VO~O, VO~2a(o)-*?{1, 0}. 
From equation (1) we obtain 
Fwauc, ewe ala} #*. 
Insertion into (2) yields 


: (2” — 1) 4a(o)'? —no2a(c)“*27~0 
which gives 


w| bs 


fo 


Further, even on the x-plane we require the potential to behave as for flow 
around a slit regardless of time. 


g~ +12a(x,)"* independent of time . 
But 
y= f2n-1 ®@D 


~ +#2"-12 a(o)'” 
mF f20A 17g (E, JF Ae 2 1M? (ae). 


Toeattain the correct asymptotic behavior of » for all time, we should take 
w= 2/3. 


Motion of Vortex Particle 


In considering the motion of the vortex particles, it is realized that this 
‘motion is a composite of the geometric magnification of the curve homo- 
geneously in time and the movement of the vortex particles along the curve. 
The value of o affixed to the vortex particle is now time dependent, o = a(t). 
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The motion of the particle is given by 
x=? &(6), where o=o(t). 


The velocity of the vortex particle, given by the total derivative dx/dt, can 
alternately be described as the mean velocity across the discontinuity surface 


ax 


a ae 


Transcription to the §-plane yields 
Z d x 
2841 76. 
Combining this with (1), the motion of the vortex particle is found to be 


do = 
tg ele 


Behavior for Large o 


For large o, Foe o/3, and the motion is given by 


t ig te oF 
dt 3 
It can be noted that for increasing time, vortex particles move to the left, 
towards the spiral core. 
Integration yields 
¢7/> ¢ = const ; 
or 7 
X(t) ~ %4(0) . 
This means that as 
a(t) ->0oo, %x,(t) > x,(0) . 
But since x,(¢) ~ 73 o, it follows that ¢t > 0. 
In other words, going backwards in time, to time equal zero, the vortex 


particles on the curled-up sheet go continuously back toward their original 
positions on the uncurled slit. | 

4. Extension to Complex Potentials : 
We assume the existence of the following analytic functions over the plane, 


C=& +7268 1 
A= OLIV, y= Ori, y= Ori : 

fos 3 ; 

: 
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where 
x: 


bole 


oe be 


~ — 


We assume further that in approaching the curve from the concave side, 
D4 approaches y, +y approaches x 


where x and x are defined on the curve as the jump and mean of the potential 
% across the curve. 


5 


Figure 7 


€-plane. Analytic continuation of 7 and {, 


Since dy/df = ®; — 1 @®,, equation (1) on the curve can be satisfied by 


dy 2, es arr 
dice os eed 
Equation (1) on the curve is satisfied simply by taking Y = const on the 
curve. We can arbitrarily choose this constant zero, thereby making equations 
(1) satisfied by 
on the curve. 
Equation (2), which is also defined only on the curve, can be satisfied by 


over the whole plane. 2» 
Combination with equation for (1) yields 


a (3) 


This equation then satisfies equations (1), (1), and (2), all on the curve. This 
equation cannot yet hold over the whole plane as it stands, however, for it 
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involves complex conjugates which would violate the assumed analyticity of 
the complex potentials. 
We consider 7 as the pee oo variable with ¢ a function of us (= C(x is 


We define the function C*(y) a 
C(x) - 


A 


et) 


oy) = 
(If ¢ were expressed as a power series in y, ¢* would be equivalent to taking 
complex conjugates of all of the coefficients in the series.) 

Since Vi = i on the curve, it follows that ¢ ae ) = C(x) on the curve. 


Our equation (3) can now be rewritten 


(3*) 


Equation (3*) now not only satisfies the required relationships on the curve 
but also allows for analytic continuation over the plane, since a function 


dependent upon ¢* can be analytic whereas a function dependent upon fe 3 
throughout the plane cannot. 


Continuation 


We have assumed everywhere 
1 
R=“ ZY: 


The mean and jump of the potential across the curve have been defined as 
= iL ¥ 
Y Mahe MW Write Aol aay 6 oy 


where y and are assigned values on the curve. Substituting for ty and “y¥ 
its definitions from above which hold everywhere, even at the curve, 


f= (ets ele eel feats 


on the curve. This has been equated to *y on the curve since a y on the 
curve. Collection of terms of the right- hand equation yields 


1 
"Le ee 


Now define-<y(C) = #y(C eam. “7 1s the value that +y assumes on the 
next sheet of the Rieman surface. a (uric 
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Since this is an equation in analytic functions and approaches the equation 


be 1 
by Se Saati 0 cae 4 


~ A 


on the sheet, it follows that the equation holds everywhere. 
As in the differential equation (3*), we use y as the independent variable. 
Substituting the symbol 4 = y(¢) and further introducing the symbol 


I= Ka) = y(t), 
we have 
(0) = (a), (4) — x) =5 (A+). 


~ 


Collecting equation (3*) and the two continuation equations yields finally, 
, = , sy iL 
3x ={20%-AlK}e, yA-WW=FU4+), M= 


where A is the independent variable and prime indicates differentiation with 
respect to A. 

The three dependent variables in these three equations are x(a), ¢(A), and 
UA). It should be noted that for 4 real, ¢(A) represents the curve itself, or 
more specifically, the concave side of the curve. 


5. Behavior at Center of Spiral 


Introduce 


Zat—b, X=x-GZl 
where ¢, is the location of the center of the spiral. The equations remain 
3X’ ={2Z*-AZ*}2', X*-XY=F A+), ZQ=2H. 6) 


Our procedure now is to develop expansions that will approximate the 
solution in the neighborhood of the origin. 

From the unrolled sheet we know that A must go to zero at the origin. 
Further, from our assumption that the sheet continually curls and expands 
homogeneously in time, it follows that the spiral must be unending in any 
neighborhood including the origin. Since our curve is defined by A real, this 
suggests a representation for Z(A) of the form 


Z = (a) 
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where the functions / and g satisfy the relations 4 


f(A) = f*(A), (A) = g*(A) , 
and asymptotically, 


The equation Z(/) = Z(A) implies that an expansion for / must start with a 


linear term in A. ; 
No assumption on X(A) need be made. This variable can be essentially 
eliminated by insertion of the assumed expansion for Z(A) into the first equation 
and then integrating to obtain an expansion for X(A). This result, together with 
an assumed expansion for /(A), can then be inserted into the second equation. 
We will assume 


Z mwa de eit 
Or, more generally, we take ‘ 


{()=a4* (1+ RO), g(4)=—o* [1+ PO) ae 


where . 
o=b-12, Rlo)= Y'0,0,, Po) = 3) Po, 

v=1 y=1 q 

and both @, and 4, are real. ‘J 
A 

Evaluation of « and pb j 

? 

Substitution of the assumed form for Z(A) into the first equation of (4) 


yields 
3X’ =A Pg tt Q2tpP eS 2ff Af. 


The two lowest order terms are 
Ale Ais GA BP gees ae ieee 2D, — p, (v—1)}o +--+] 
Zg5" Go 24 a7 8 gt 10h see]. 


Then the first two terms in the differential equation are 


Nae ak 8 Be ke |-o° +2 (—e si 4] ot + +]. 


Integration gives the leading terms in X 


fear Ls 2 A2 (ed 1 29 
X= 5 a? 2 Ae) — 2a op Oa a 


i) 
R 
| 
wD 
— 
> 
e 
+ 
a) 
— 
Q 
| 
a 
+ 
: ee 
ot ee 
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We now assume an expansion for /(A). 


t= alt+ Sol, 
y=] 
Substitution into the second equation, 


X() —X(@)=-F+), 


yields 
X(l) — X(A) 
: os s : 
=3 a? BA | —1, 0? + | ko 2 & ot 1 + 2(-o +4 m)h+---| 
=-A-fmdo—-. 


Since this should be an identity in A, we find by equating both exponents and 
coefficients of these first two terms from either side of the equation 


Za—-B=1, Sa?ftbm=l, m= (B+1) x. (6) 


We follow a similar procedure in the third equation 


Taking the logarithm 


. L 1+R : 
F4{s[1+ P(s}— ofl + P(o)] }+a log=- + logs i ae aed 
where 
s=o[1+ So : 
v=1 
Again, by insertion of the expansions, the two lowest terms are 
i 
s-t = o 1 =—f2%,+ [—B 2 + Abt) a] o+-:-- 5 
aw log, == 1h MeO =": 
or one 
1B xwy+ [ary + 6B my — “FE dl Of =2n1. 
Satisfying the identity here yields’ 
1 i 
mae, ma bt Natagna. (7 


Z 
Bs 
= 


ZAMP VII/22 


—ee ee 


. 
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> 

{ 

There are two determinations now for x); one from equation (6) and one 
from equation (7). Equating them yields 


0, = 05 052 


and then f = 3 from the first equation of (6). 


Comparison with Kaden 


Our result yields, to lowest orders terms, 
DODGE Ae ee 


where for real J, this defines the spiral. 

Kapen’s result for the form of the spiral in the neighborhood of the 

origin is : 
r~c p78 


ee 


or é' 
Ve 9-23 ¢ -10 ; t 


In the neighborhood of the origin, then, our results are equivalent with the 
correspondence 


DP =A) 


We find further that the mean of the potential, x1); has a pole of order 3 
The series for /(A) gives 


Py teensy) + (57+ Gai) a+ At + 


A map of the Z-plane onto the A-plane would then show both sides of the spir 
to first cusp at the origin and then separate as shown. 


Areal 


Figure 8 
Map of the spiral onto A-plane. 


UJ 
U 
ies ei 
Vv 
0 es satis ee DS iit * 
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6. Determination of Higher Order Terms 


The procedure for determining higher order terms is clear. Substitution of 
the series into equations (5), and equating coefficients, yields two independent 
determinations of x,, in terms of 9, and ,. @, and p, are then successively 
determined. Several of these computations have been carried out and appear 
in the appendix. Q, $2, @3, fs have been computed and the spiral including 
these terms has been plotted. 

It will now be shown that successively all coefficients are uniquely deter- 
mined. > 


Successive Determination of On, Pn 


Consider first the equations obtained by equating the real parts of the two 
determinations of x,,,, where m = 1. 

In equation (5), for X’, the contributing term is — A f? g’/3. Of this, the 
contributing factor is 
. Bee = (n — 1) bah or-t. 
Then 


| 


 A[X(s) — X(0)] = ++ + enga t+ = 2 en — (=D Pa fra te dorten, 
Xn+1 2{ On (1 1) Pn} % + ae 


‘In equation Z(/) = Z(A) e?*, the contributing term is s~1 P(s) — 0-1 P(o) , 


: Se ee 13 eae be 8 an 1) ot, +s Lote; 
Nass = (W—1)pn%y+-- 

-Equating real parts of #41, gives 

; Go Ne te 


Consider now the equations obtained by equating the imaginary parts of the 
two determinations of x4» - 
In the equation for X’ there will be two contributions from 7 (2 f? g’/3) and 


one from — 4 f? g’2/3 . They yield 
I [X() — X()] 


24 24 a 
: op ed [—%¥ne2 Fo + tn to hg 2 aC ne (n — 1) pa fm + 


+ Im %{—2 [en — (% — 1) Pal + ++} onthe, 


eee ee ey RN 


- 
: 
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Substituting from real part of %n41 and imaginary part of x, in this equation 


~ 


yields 


An+2 = = fe 2{ On oe (n a 1) bn} a Ft {2 en - (n ar 1) Pu} my 


ese oat 2 On +3 (n—1) pyhmy +e. 


From equation Z(!) = Z(A) e'”* there is a contribution from each of the 
factors. Given in order, they are 


[oes $26 Ming bee bE 3 0 On ty to {—3 Mags + ++ +3 (4 — 1) pn Im me 
4 etpiottl te are 1 


Substituting for x,,, and Im x, gives 
cere ze 
na = (WN) Pa Mt Go 7 Ont + (HN Pata too 
2 3 
na = G15 On 2 (n— 1) Ba} my oo 


Equating imaginary parts of x,,5 gives 


— a aed 


(3+ 4) Qn—2 (0-1) Py =o 


Finally, then, the determination of 0, , comes from 


es 


2 — (n— 1) On| 


3n+4 =2(n—1) Pn 
For m = 1, the determinant vanishes only for 7 = 1. Therefore, successively all 
coefficients, Qn, Pn, %n, Can be uniquely determined after 0,, $1, %, %, and Im 3 
have been determined. | 

Now, Q1, %, %2 have already been determined by choosing «, B to satisfy 
equations (6) and (7). p; remains arbitrary and is seen not to influence succes- 
sive coefficients. This is expected since #, establishes a rotational orientation of 
the center; since the curve is an unending spiral there, such an orientation is 
unnecessary. | 

The only critical test, then, for the feasibility of successive determination 
of all coefficients lies in the successful determination of Im #3. Interestingly 
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enough, the two independent determinations of Im %3 agree identically; both 
yield Im x, = 11 x?/3. 


Successive determination of all other coefficients uniquely is therefore 
established. 


‘ 
&> 


v4 


Figure 9 
The rolled vortex core. 


Appendix 


Equating Two Determinations of % 


. é i Ot 11 og. 1! ear Langs 
x3 imaginary: Rl ee eae real: 2 Q2 — Cao ts 


: : 5 
“, imaginary: 5 0;— p2= | *1) 


iil 
(on = 4g *f. be = 35 74). 


aes: ee ee ; 
(c= i6g “i> Ps= Fen “4 


For the case where 


st > 4% be 


ee wey 
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the solution Z(A) is 
Ai) = ape [t+ 35 CE) + see Ye 


xe [14 (St) as + sex (3) +] 


For A real, the form of the spiral has been plotted for a range of 


0:30 <4 0-65. 


The author wishes to express his gratitude to Professor K. O. FRIEDRICHS, 
New York University, under whose guidance this work was done. 
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Zusammenfassung j 


Die vorliegenden Untersuchungen befassen sich mit dem Problem der zwei- 
dimensionalen zeitabhangigen Strémung, welche mit dem Abrollen einer semi- 
unendlichen Wirbelflache verbunden ist. Obwohl bereits vorliegende analytische 
Untersuchungen an diesem Problem Naherungslésungen fiir die Spiralform der 
aufgerollten Wirbelflache ergeben, so diirfte es doch schwierig sein, aus diesen 
Ergebnissen das zugehorige Str6mungsfeld zu erhalten. In dieser Arbeit werden 
Differentialgleichungen fiir den Mittelwert und den Sprung des komplexen 
Potentials abgeleitet, welche zur Beschreibung des Strémungsfeldes iiber die 
gesamte Ebene benutzt werden kénnen. Fiir den Spiralkern stimmt die hier 


angegebene Lésung mit der von KapDeEn iiberein. 
> 
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A Graphical Method for Determining the Coefficient of Viscosity 
of Newtonian Liquids Using an Oscillating Cylinder Viscometer 


il 
By Avr A. K. IBRAHIM and ABDEL Monem I, Kasret, Alexandria, Egypt?) 


It seems possible to represent our results and those of OLDRoyD, STRAW- 
BRIDGE and Toms [1]*) by a linear graph (which is most suitable for practical 
work) from whose slope or intercept the coefficient of viscosity 7 of the Newto- 
nian liquid may be estimated. 

Some little time ago, we have pointed out the equation (2): 


ta 2 Gy &—a@- Pot 6) \2]-12 , (b—a 
Rae. saben.» “gle [1 Gea i eal (@) 


where 9, is the angular amplitude of the suspended cylinder, gy, is the angular 
amplitude of the outer cylinder, } is the radius in cm of the outer cylinder, a is 
the radius in cm of the inner cylinder, J is the moment of inertia in g cm? of the 
suspended cylinder about its major axis, t is the wire torsionconstant in gcm2s~?, 
@ is the angular velocity in s~! of the outer cylinder and L in cm is the depth 
to which the suspended cylinder is immersed in the liquid in the outer cylinder. 

This equation represents the coefficient of viscosity 7 in g cm™} sec} of 
Newtonian liquids with considerable accuracy for all values of 7 >1 (Poise). 

Such a formula is of wide application and would be very valuable by the following 
representation. 

The manner in which it was shown that (1) represents the experimental 
results may be of interest, if (q/6))? be denoted by Y and [(I w? — t)/w]* be 
denoted by X, then the relation between Y and X is linear and may then be 
written 


2 
vm (=) be (2) 
where j . 
ee us 
K= aa = E a3 e 3 [(b a) /(b+a)] (3) 


K is an apparatus constant which may be either calculated from (3) or obtained 
by using a liquid of known coefficient of viscosity. The slope of the line gives 
(K/n)? and the intercept gives (/K)? from which the coefficient of viscosity 7 
may be deduced. If t = 0, therefore a plot of (9/8) as Y against I? w% as X 
gives a straight line. ; 
. The Figures which are drawn to scale, show our experimental results (in case 
‘of paraffin oil of density 0-875 g cm~® at 20°C) and those of OLDRoyD ¢¢ al. 


1) Physics Department, Faculty of Science, University of Alexandria. 
2) Numbers in brackets refer to References, page 345. 


re 
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100 


fe) 5 10 a) 20 


Figure 1 ; 
Variation of (09/@ )~2 with [(Iw? — t)/w]* for paraffin oil (OLDRoyD). 


(in case of liquids paraffin B. P. of density 0-877 g cm™* at 25°C) by the pro 
cess outlined above. 

The obtained graphs were very accurately linear which is distinctly superior 
in point of accuracy, the values of the slope and intercept with the X-axis were 
read directly from the graph. 

The following table shows the results for the liquid examined. 


: 7H in P 

3 ntercept with the 

K Slope (K/7) X-axis (n/K)? fenuoee from inter- 1 oe 
cept 

| 

0-6656 x 10-4 0:1427 x 108 7:041 x 10-8 1:762 1-766 1-764 | 

OLDRoyD [1] 

1-928 x 10-2 1-009 x 10-4 ise alle 1-911 1-928 1-920 

IBRAHIM and 

KaBIEL [3] 


It is seen from the above table that the coefficient of viscosity 7 is obtained 
with a distinct degree of accuracy which confirms the validity of the method 


ar shows the very accurate representation of equation (1) to the experimental | 
acts. : 


Also the figures show that 6,/y) =1 when X =0 or, 


‘ 


{ 
: 


i. e., at a frequency equal to that of free oscillation of the wire and the suepeude 
cylinder. 


oye 
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(85:2) 


Figure 2 
Variation of (89/@)~* with [(J w* — t)/@]? for paraffin oil (Iskanim and KastEt). 
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Zusammenfassung 


Ein Verfahren wird beschrieben, welches erlaubt, die Viskositat zaher Fliissig- 
keiten zu messen. Durch geeignete Wahl der dargestellten Gréssen k6nnen die Mess- 
ergebnisse der Verfasser, wie auch diejenigen von OLDROYD, in einem Graphikum 
linear wiedergegeben werden, wobei die Neigung und der Abzissenabschnitt den 

‘Viskositatskoeffizienten zu ermitteln gestattet. 


(Received: March 3, 1956.) 


Longitudinal Vibrations of a Composite Bar 


By VActav VopiéKa, Pilsen, Czechoslovakia*) 


There is little difficulty in solving vibration problems for uniform bars, but 
mathematical discussion of longitudinal vibrations of composite bars is in 
general more complicated. This paper gives an outline of such calculations. It 
‘may serve as a kind of guide in treating other related problems of technical or 


physical significance. 


1) Technische Hochschule. 
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1. Statement and Mathematical Formulation of the Problem 


A bar of length J =/,+/,+-+:-+1, and constant area g is composed of n 
parts of lengths J,, Youne’s moduli E;, and densities 9, (k = 1, 2,...,m). It is 
at rest with the end + = 0 fixed. At the time ¢ = 0 a constant axial force P is 
applied at the free end ¥ = 1. To find the subsequent longitudinal vibrations in 
the bar. 

Denoting by u = u(x, t) the displacement at the point x at any time ¢ and 


putting : 


5905 Spee lay ee ee (1) 


“n=1 


we obviously have to solve the system of differential equations 


07Uy, 07*Uz, _ — E;, ‘ 
tae i arate SRS es I Neer (Bh Zee, 70 
with the initial conditions : 
U.(%; 0) = 07, Bes y = 0) Spey VR Soe (i Lee (3.1)% 

t 

and with : 
u(0,f) = 0 he Oe (3.2.1) © 

O0U>(S, ¢ Onesie Ee { 

Up(Sps t) = Uysa(Sps t) , oe iv ek soy wt) » $>0, & sai (3.2.2) 
(Re aed Ze tie Le = Be i 

EA Foe bal Jaa 

Eye = GE? E208. (3:23) : 

2. Solution of the Problem : 

Putting 4 

cS : 

Ux(%, p) = pfePtuy(x, t) dt (k= 1,2,...,n) (4) 

F } 

3, 

. . . 

we easily obtain from (2) — (3.2.3) the following subsidiary equations: a 
| 

@U;(%, p) _ p? . 

dx? of a? Uy, P), Spax<*<s, (k= 1,2,...,n), > 

(05 p\-= On (6.1) | 

2 dU (sx P) dU sa(Siv P ; 

Ux (Sz, p) — Urs a(Sps Pp) , ax = & Ee (R = 1, 2, sie We — 1) GI 
(6.2) 
dU,(l,~) _ P | 

dx Guise 


This gives 


U,(#, p) = A, cosh 2* + B, ems (A = 1, 2.5... #) 
ko ay, 
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where the constants A,, B, must satisfy the conditions 


A,= ON (atk) 


, 


P Sx ; Ss Sp : Sh 
A,, cosh + B, sinh jl ee Ay... cosh Ly 4"By anh Se 
a a ay k+1 

k i Vea aaa 


Ss 
A,, sinh a : ih — i, (4 eet SINK oes + By. cosh p “2. : (ria) 
ay Qotn Anay 
a 
he k+a Qh+a Riel a hod ae 
. ay Of ( 
: l 
A,, sinh p 1B ,, eosnlS ee ae (7.3) 
Ay an P qd ay On 
Writing (7.2) in matrix form 
aes A, 
= M,(p) teas 102 Neh een) 
Byis B,, 
cosh P Sk cosh P Se S sinh P Sk eee 22 ; 
ay, Beit k ay, Bes4 
sinh P Sr cosh E Si cosh P Si sinh la 
k kh. k ke k+1 
M,(b) = H , “ : 4 (8) 
Seok! * Gan + erie Poe Oe : 
ay, Appa k ay Ge44 
— sinh Pain P Sk st cosh PSE P Sk 
k Dea k ke Veta 
(pce tsi, ,n— 1) 
and considering (7.1) we obtain 
A pay O 
= 0,(p) yay a ee heen 
Bry B, 


Q,(p) = M;,(p) My_s(>) ++ Ma(p) Milf) (F=1,2,...,2-1)- (8-1) 
Applying the notation 
Ox(P) = ||) || (v= 12; k=1,2,..,0-1) (9.1) 
_ we also can write 

Ana = (22) Bu Bess = PP) Br (R= 12." -— 1, (9:2) 

_and we find from (7.3) at once : 

i 1 

; ie Ce p [ato pe o + gh (p) cosh 2 | a 
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From (7), (9.2) and (10) we finally obtain 
Pp sinh as 
UP) ee ; , (1.1 

: eee p [a2 (p) sinh £5 4 22) (p) cosh 2] 

n 


n 


p __ 2(P) cosh P* + gf*(p) sinh 
J Devt D+ (11.2) 


q An, On p ee (p) sinh - BE + qa" ”) (P) cosh p pt) 


Unsal®, p) = 


MR Geko 


and our problem is then reduced to the determination of the elements in the 
second column of (8.1). Apart from exceptional cases this task is very hard. 

The expressions (11.1) and (11.2) represent the Laplace transforms of the 
wanted functions u,(¥%, ¢), which are to be determined by the well-known methods © 
of operational calculus. For few beginning values of n (e. g. for m = 2, 3) this can 
be done quite generally, but the results are complicated. 


t 
We now proceed to the treatment of the case characterized by h, = 1_ 


3. Special case 


(k = 1, 2,...,m— 1), i. e. according to (2) and (7.2) by 

EL, Gy = const; (k= eee (12)% 
Here it is possible to express the solution in closed form. j 
We first have by (8) . 
cosh o,, sinhp o;, 1 1. 

M,(p) = , on = 5x ( — ) (ss I, Zs n — 1) 
sinhp o,, coshp o, be ed 2 : 
(13) 


and from this by (8.1) 
cosh p(o,+ a2), sinhp(o, + a4) | 


Q.(p) = , 
sinhp(o,+ 0), coshp(o, + o,) 
But generally : 


v v vy+1 v+1 I 

cosh p py: 6,, sinhp oP ae cosh p pa o,, sinhp ms Oy ‘ 

x=1 x=1 | “x=1 =1 5 

My.4 (P) a 1 3 ; 
V+ y+l1 


sinh p Bye Gar coshp 37a, 
x%=1 “x=1 


| sinhp )"o,,, ee 
x=1 


QilP) = ; “oe (b= 1, 2h. 50 = 1) oo ee 


k 
sinh p ay o,, coshp oy G,, 
a=] x=1 


‘ 
H 
and so we obtain from (8.1) at once 
k k 
en Oy, sinhip os; (ahs 
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Herewith equations (11.1), (11.2) take the form 


; px 
sinh“=— 
U,(*,~) = 7 _. eet ae Ab 
q ay Qn ) Feat 
reap (2 Fu) 
ay “w=1 , 
‘ k 
sinh, : + o 
U; j Gs ~ J 
rail%, P) ales. - ee es Cte Moree on iverat pi} 
n 


and from this we find (using a table of Laplace transforms or by contour inte- 
gration) the following complete solution of our problem: 


a es 8 1 = 
U,(x, t) = a — an 6. 
q an On ay a? an, * 


“x=1 
oo : f 16.1) 
(—1)r (2r4 1) nx (Aree t) we ( 
Pa (2+ 1)? sin - or rte cos : rey , 
2a,{—-+ Ox Boom Se Dez 
On e=1 an x=1 
k j n—1 
tn % 8 l 
Ups (%, t) = + 6, — —|— + Oy 
: k 
(27+ 1) «( + “| 
xy Bey ses, Bpss » (27y+1) at (16.2) 
& (27+ 1)? f7. wal On alg ae 
2( | 2 (z : a) 
“x=1 n “n=1 
(ee 2 eee poco 


Remark. For »=1 formula (16.1) gives the well-known expression for 
longitudinal vibrations of a homogeneous bar. 


Zusammenfassung 


Mathematische Untersuchung von Langsschwingungen in zusammengesetzten 
Staben ist ein verhaltnismdssig kompliziertes, jedoch industriell wichtiges 
Problem. Die vorliegende Arbeit zeigt an einem typischen Beispiel die Art der 
dabei vorkommenden Schwierigkeiten. Sie kann auch bei rechnerischer Behand- 
lung anderer Fragen aus der Schwingungstheorie zur Orientierung dienen. 


_ (Received: December 16, 1955.) 


b anethin | 
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Uber eine Abschitzung der Amplituden in freien Schwingungs- 
problemen verinderlicher Kreisfrequenz 


Von AUREL WINTNER, Baltimore, Maryland, USA?) 


Es sei w(t) eine gegebene (reelle und stetige) Funktion der Zeit, und es sei 
x(t) eine beliebige Lésung des Schwingungsproblems 


x" + w(t) x= 0. (1) 


Die triviale Lésung +(¢) = 0 der Differentialgleichung (1) soll jedoch ausgeschlos- 
sen werden. Explizit kann man (1) meistens nicht losen. Dies ist der Fall bereits 
fiir die einfachste im Problem der Frequenzmodulierung vorkommende Koeffi- 
zientenfunktion, namlich 

w(t) = (a + b cost)”, 


wobeia und b gegebene, den Ungleichungen 0< | b |< a geniigende Konstanten [ 
bezeichnen. Jede Information ist daher von Belang, die die Kenntnis einer 
Lésung x(¢) von (1) nicht voraussetzt und dennoch eine Aussage macht tiber die 
numerischen Verhaltnisse in dem zu einer beliebigen Lésung x(¢) von (1) gehé-— 
rigen (¢, ¥)-Diagramm. Dabei kénnen solche Daten wie die (zu beliebigen posi- 
tiven Werten h gehérigen) Integrale 


h h h 
Mo = | 0(t) GES op, = ft w2(t) dt, py =# w2(t) at 
0 0 0 


als bekannt gelten, da diese bestimmten Integrale ein fiir allemal, das heisst un- 
abhangig von der Wahl der Partikularlésung +(¢) von (1) durch Reihen bestimmt 
oder numerisch tabuliert werden kénnen, sobald die Koeffizientenfunktion von 
(1) fir 0<t#<h gegeben ist. Im folgenden soll ein Resultat dieser Art her- 
geleitet werden. 

Es sei ¥ = x(t) eine Lésungskurve von (1), die fiir ein gewisses ¢ = #, die 
t-Achse schneidet [es gibt eine einfach unendliche Schar von solchen Lésungen, ~ 
da mit #(t,) = 0 nur eine der willkiirlichen Integrationskonstanten *(é,), %’(¢,)_ 
gegeben ist]. Man verfolge die Kurve 7 = x(¢) von ¢ = ¢, an, und es soll voraus-_ 
gesetzt werden, dass dann die Kurve wieder die ¢t-Achse erreicht. Es sei #, die ‘ 
erste, auf ¢, folgende Abszisse, fiir welche dies zutrifft [¢, hangt von der Wahl der 
Integrationskonstante +’(t,) ab]. Man kann den Anfangspunkt der ¢-Achse derart 
wahlen, dass 7, gleich 0 wird. Dann ist h> 0, wenn h = ¢, gesetzt wird. Danach — 
gilt x(t) = 0 fir ¢= 0 und fiir ¢= h, aber nicht fir 0<¢*<h. Da mit x(t) auch 
— x(t) eine Losung von (1) darstellt, so kann man annehmen, dass «(t) fiir 0< t<h 
positiv ist. Dann folgt aus (1), dass x"(t) fiir 0<¢t<h nirgends positiv wird. 
Bezeichnet daher J die Strecke zwischen ¢ = 0 und ¢ = h auf der t-Achse, so ist 
der Zu 0S t<h gehorige Teil des (¢, x)-Diagramms eine « Halbwelle», die aus i 
einem mit Riicksicht auf x(t) <0 konvexen Bogen K besteht, der die beiden 
Endpunkte von J verbindet und sonst oberhalb J liegt. Man beachte, dass derd 
héchste Punkt von K nicht iiber dem Mittelpunkt von J zu liegen braucht. 

Es stellt sich nun heraus, dass, wenn f den Flacheninhalt der durch J und Kk 
umgrenzten Halbwelle und m die in ihr gelegene Maximalamplitude bezeichnet, — 


— 


1) The Johns Hopkins University, Department of Mathematics. 


. 


. 
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also die positive Konstante m durch 


m = max x(t) 
O<t<h 


definiert ist, fiir den Quotienten m/f eine Abschatzung a priori gilt, namlich 


+ pagal si (3) 


Insbesondere ist also 
fy >t, (4) 


da auch f positiv ist. Dabei sind ,, mw, die durch (2) definierten Werte. Letztere 
sind positiv, da sie nicht negativ sind und nur dann gleich Null werden kénnen, 
wenn die Gleichung (1) sich auf ¥” = 0 reduziert. Dann ist aber a(t) eine lineare 
Funktion, die der Annahme einer Halbwelle widerspricht. 
Ausser (3) gilt auch 
m 4 


ae ee 5 
i Mo ©) 
wobei jw) wieder durch (2) definiert ist. Es wird sich aber zeigen, dass die Ab- 
schatzung (5) weniger besagt als (3). 

Die vorangehenden Tatsachen kénnen auf die folgenden beiden Ungleichungen 
gegriindet werden: 


h 
h <[t(h—2) w%() at (6) 
3 


und 
hm< 2f. (7) 


Die Ungleichung (6) ist eine Verscharfung eines klassischen Satzes von Lia- 
PouNorr, die Herr Dr. P. HARTMAN und der Verfasser in einer gemeinsamen 
Note bewiesen haben?). Um sich andererseits von der Richtigkeit der Unglei- 
chung (7) zu tiberzeugen, geniigt es, den Flacheninhalt des Dreiecks D mit den 
Ecken 

(é, ¥) = (0,0), (@*)=(0,h), ( *) = om) 


zu betrachten, wobei ¢, denjenigen Punkt der ¢-Strecke J bezeichnet, in welchem 
die Ordinate x(t) des konvexen Bogens K ihr Maximum m erreicht. Da D in dem 
durch K und J begrenzten konvexen Bereich enthalten ist und letzterer den 
Flacheninhalt f hat, so ist (7) klar. 

Um nun (3) zu beweisen, beachte man, dass (6) in der Gestalt h< hm, — My 
geschrieben werden kann, wenn man die Bezeichnungen (2) benutzt. Die letztere 
Ungleichung besagt aber, dass w,< h(u,— 1) gilt, woraus einerseits (4) und 
andererseits 

ae 
Hy — 1 


-folgt. Da Zahler und Nenner des Quotienten (8) positiv sind, so ist es klar, dass 
(7) und (8) zusammen genau zu der Behauptung (3) fiihren. 

Die Behauptung (5), die bereits bekannt ist*), kann als eine weniger scharfe 
Form der soeben bewiesenen Beziehung (3) betrachtet werden. Um dies einzu- 


<h (8) 


2) Ein einfacherer Beweis als der urspriingliche findet sich in einer Arbeit von Z. NEHARI, 
Amer. J. Math. 76, 690 (1954), wo weitere Literatur angegeben ist. 
3) P, HarTMAN and A. WINTNER, Quart. appl. Math. 10, 175 (1952), 
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sehen, beachte man, dass die Funktion ¢(h — t) fur O<¢#< h positiv ist und ihr 
Maximum h2/4 betragt, da es fiir ¢ = h/2 erreicht wird. Demnach folgt aus (6) 
und aus der ersten der drei Definitionen (2), dass h< py h?/4, also 

bia <h (9) 
Ho 
ist. Nun ist es aber klar, dass (9) und (7) genau so zu (5) fiihren, wie (8) und (7) 
zu (3) gefiihrt haben. Dabei ist (9) im wesentlichen der obenerwahnte Satz von 
LIAPOUNOFF. 


Summary 


Consider any differential equation of the form given under (1) below and, 
if x(t) is any solution of (1) which is real-valued and does not vanish identically, 
let C be a ‘half-wave’ a<t<b of the curve ¥ = x(Z) in a (t, ~)-plane [by this is 
meant that x(a) = 0 and #(b) = 0 but | #(t)| > Oifa< t< 6). Then an a priors 
upper limitation, involving the ‘moments’ of the coefficient function of (1) ; 
(cf. the three integrals in formula (2) of the text, where a= 0 and b =k); is 
obtained for the ratio m/f. Here m denotes the maximum amplitude occurring in > 
C, and ? the area of the region surrounded by the segment a<#<b on the 
t-axis and the graph of C in the (#, ~)-plane (this region is always convex). - 


(Eingegangen: 10. Marz 1956.) 
Varia — Miscellaneous — Divers ‘f 


Frihjahrstagung der Schweizerischen Physikalischen Gesellschaft : 
vom 5. bis 6. Mai 1956 in Zofingen 4 


Berichte iiber die Sitzungen fiir angewandte Physik und Mathematik — 


Stabilisierter NF-Verstarker und phasensensitiver Detektor, von T. 
KNELLWOLF und J. OEsER, Ziirich}). : 
Der Vergleich von Signalgréssen bei einem Kerninduktionsspektrometer er- 
fordert eine Apparatur, bei der die Empfindlichkeit tiber lange Zeiten auf min-— 
destens einige Prozent genau konstant bleibt. Dazu ist es unerlasslich, dass das — 
ganze elektronische System gut stabilisiert ist. _ 
Nach der HF-Gleichrichtung ist das Kerninduktionssignal eine periodische — 
NF-Spannung, deren Grundfrequenz gleich der Modulationsfrequenz wy ist, mit — 
welcher das statische Magnetfeld moduliert wird. } 


1. Abgestimmter NF-Verstirker 


Bei der iiblichen Anordnung (BLOEMBERGEN, Proctor und andere) kommt 


das Signal in einen mit TT-Riickkopplungsnetzwerk versehenen frequenzselek- 
tiven Verstarker, der auf die Grundfrequenz des Signals abgestimmt ist und 


1) Physikalisches Institut der Universitat. 


Vol. VII, 1956 Varia — Miscellaneous — Divers _ 353 


dessen Bandbreite einige Hertz betragt. Der Nachteil dieser Schaltung liegt 
darin, dass die Verstarkung fiir die Signalfrequenz nicht stabilisiert ist. 

Eine frequenzunabhangige Riickkopplung (Figur 1) tiber der TT-Stufe be- 
wirkt zwar diese Stabilisierung, erhéht aber gleichzeitig (bei gleicher Totalver- 
starkung) die Bandbreite. Es ist daher nétig, die Bandbreite des Verstirkers ohne 
frequenzunabhangige Riickkopplung ausserordentlich klein zu halten. 

Die totale Verstarkung ohne frequenzunabhangige Riickkopplung betragt 


___ & G 
1+ Brr(o) G, ° 
Mit frequenzunabhangiger Riickkopplung betragt sie 


Gy ie «2 Gy — —G, Gy 
1+ BGy 1+ Brr() G,+ 6 GG, © 


o+ 


Gy= G, Grr = 


Figur 1 
TT-Verstarker mit stabilisierter Verstaérkung (Prinzipschaltung). 


Die Bandbreite des TT-Verstarkers betragt 


4/3_ 


Awrr= om —G 
1 


‘Mit frequenzunabhangiger Riickkopplung wird sie 
Awr = 6G, G, Awrr = 4/3 BG, om, 


wenn die Bedingung Bf G, <1 <f G, G, erfiillt ist. Dabei bedeuten: 
G, Verstarkung der Stufen des TT-Verstarkers ohne TT-Rickkopplung, 


Grr Verstarkung des TT-riickgekoppelten Verstarkers, 
G, Verstarkung der Stufen, die dem TT-Verstarker vorangehen bzw. folgen. 
In unserer Ausfiihrung betragt G, = 1000, G,= 10. Bei einem Riickkopplungs- 


faktor von 16 ergibt sich eine Bandbreite von + 4 Hz. 


A 2. Phasensensitiver Verstarker 


_ Im phasensensitiven Verstarker wird die Bandbreite durch Homodynisierung 
‘und im nachfolgenden Gleichstromverstarker durch RC-Glieder weiter herab- 
gesetzt bis auf 1+ 1/20 Hz. Die meisten phasenempfindlichen Verstarker ver- 
‘wenden Mischréhren oder auch Pentoden mit Fanggittermodulation. Die Nach- 


ZAMP VII/23 
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j 
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teile dieser Schaltungen sind: Die Konversionssteilheit andert sich mit der Ampli- 
tude der Mischspannung und ausserdem auch mit dem Altern der Réhren. Um 
von den Variationen der Réhrengréssen unabhangig zu werden, verwenden wir 
Diodenmischung zur NF-Demodulation (Figur 2). 

Die Primdrseite des Transformators I liegt zwischen den Kathoden zweier 
Kathodenfolger, an deren Gitter die Misch- bzw. Signalspannung liegt. Am 
Transformator I liegt also die Spannung 


Vu coswmyt+ Vs cos(wyt+ @), 
aus der durch Gleichrichtung eine Spannung 
Vn + Vs cos 


entsteht. Der Transformator II tragt eine Spannung Vj, cos wy t, die entstehende 
Gleichspannung ist Vy. 

Die Differenz der beiden Gleichspannungen, die jetzt von Vy unabhangig ist, 
gibt man auf das Gitter eines stabilisierten Gleichstromverstarkers. Die in Gegen- 


+ 300 


heii tan od 


Figur 2 
Phasenempfindlicher Detektor 


> 
t 


taktschaltung verwendeten Transformatoren miissen zwei moglichst gleiche 


Sekundarwicklungen besitzen. ; 
Ergebnisse : 


1. Linearitat 


Der phasensensitive Verstarker ist i 
als auf 1% 


linear. 


: n einem Arbeitsbereich von 20 mA pessell 
inear. Im totalen Arbeitsbereich von. 40 mA ist er besser als auf 2% 
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2. Unabhangigheit von der Refeves zspannung 


Variationen der Mischspannung von + 30% bewirken eine Verschiebung des 
Nullpunktes um 0,1 mA (ganzer Bereich des Registrierinstrumentes 20 mA) 


3. Die Empfindlichkeit 


Uber das Impulsverhalten von Flachentransistoren, von E. BALDINGER 
W. Czaja, M. Nicoret, Basel?). 

Zur Beschreibung der Wirkungsweise von Flachentransistoren ist das in 
Figur 1 dargestellte Ersatzschema fiir Impulsschaltungen besonders geeignet?). Die 
Parameter F,, R,, S und C lassen sich aus dem ebenen Transistormodell®) be- 
rechnen. Das Ergebnis ist ebenfalls der Figur 1 zu entnehmen. 


Figur 1 


Ersatzschema fiir Flachentransistoren. 


Da bisher Angaben dariiber fehlten, wie gut dieses Modell das Verhalten 
eines wirklichen Transistors beschreibt, sind hier an zwei méglichst verschiedenen 


1) Physikalisches Institut der Universitat, Abteilung fur angewandte Physik. 

2) In Impulsschaltungen ist der Belastungswiderstand oft klein; deshalb kénnen der Kollektor- 
Emitter-Widerstand (R,) und der Kollektor-Basis-Widerstand (Rg) veriachliassigt werden (gestri- 
chelt eingezeichnet). Fiir den Philips-Transistor OC 71 betragen zum Beispiel bei J, = 1mA: 
R, ~ 7 MQ und Ry ~ 100k2. 

3) A. W. Lo et al., Transistor Electronics (Prentice-Hall, 1955), Kapitel 8. 

C. W. MuELterR and J. I. Panxove, A p-n-p-Triode Alloy-Junction Transistor for Radio- 
Frequency Amplification, Proc. I. R. E. 42, 386 (1954). 
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th-Kurve Philips 0C71 


~<—Philco $B4100 
80 


60 


40 


20 
10 : 
SSS a nn a, ee 
0 1 2 3 mA 
Figur 2 


Steilheit als Funktion des Emitterstromes. 


Transistoren (Philips OC 71 und Philco SB-100) Messungen ausgefiihrt worden, 
die zeigen, dass die modellmdssig zu erwartenden Abhangigkeiten vom Emitter- 
strom I, (vgl. Figur 1) befriedigend erfiillt sind. Die Resultate sind den Figuren 2 
bis 4 zu entnehmen. 

Es ergibt sich, dass innerhalb der fiir Anwendungen erforderlichen Genauig- 
keit das Ersatzschema von Figur 1 das Verhalten eines Flachentransistors tiblicher 
Konstruktion richtig wiedergibt. Insbesondere ist demnach die Steilheit aller 
Transistoren annahernd gleich, da sie ausser vom Stromaufteilungsfaktor « 
(x = 1) und dem Arbeitspunkt I, nur von der Einheitsladung q, der Boltzmann- 
Konstanten k und der absoluten Temperatur T abhangt. Dagegen erscheint die 
Kapazitat C als das bestimmende Element fiir das Frequenzverhalten des Transi- 
stors; ihr Wert hangt neben der Diffusionskonstanten D, der Minderheitsladungs- 
trager in der Basis auch noch von der Dicke W der Basisschicht ab und ist deshalb 
eine stark konstruktionsbedingte Grdésse. 

Ferner sind die Ausdriicke C/S und R,C in guter Naherung unabhangig vom 
Arbeitspunkt und deshalb fiir den Transistor charakteristisch. Die gemessenen 
Werte sind in der folgenden Tabelle zusammengestellt: 
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| Philco SB-100 


~1,5-10-9 


18+ 20-5 ~4 +1078 
g empmersfa Ko}. / Bias 08 


Philips OC 71 


C/S [s] 
RC [s] 
Auflésungsvermégen4) [Hz] 


| 


Zum Vergleich sei das Auflésungsvermégen von Zweifachuntersetzern ange- 
fiihrt, das mit denselben Transistortypen erzielt wurde‘). Es lasst sich daraus 
folgern, dass das mit SB-100-Transistoren erreichte Auflésungsvermogen yon 
5 MHz noch wesentlich verbessert werden kénnte. 

Wir sprechen der Stiftung Hasler-Werke, durch deren finanzielle Unterstiit- 
zung die vorliegende Arbeit zustande kam, unseren Dank aus. 


pF 
100 


Philco SB8400 
Up= -3V 
Tx 22°C 


Phileo SB100 


Up=-3V oi 
T~ 22 
2 
0 1 2 3 mA 
—Se 
k2 
2 


Philips 0C71 
Up=-6V 2:10 
T~ 22°C 


CRy+R) 


Philips OC 71 
Up= -6V 
T ~ 22°C 


a5 i 2mA* 0 i 2 Lp 3 mA 
mA~ 4 2 | Q7 = 05 e 
Figur 3 Figur 4 
Eingangswiderstand als Funktion Kapazitat C als Hyon, 
des Emitterstromes. des Emitterstromes?), 


4) E. Batprncer und M. Nico.et, Untersetzer mit Fldachentransistoren, ZAMP 6, 503 (1955). 
5) Anmerkung bei der Korrektur: Der ausgepragte Knick in der Kapazitatskurve des Typs 
SB-100 trennt die Gebiete kleiner und grosser Injektionsdichte der Minoritatsladungstrager in der 
Basis. Vergleiche zum Beispiel das inzwischen erschienene Werk der RCA Laboratories Princeton 
Transistors I (1956), S. 48. 
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Gedanken zum heutigen Stand der Reaktortheorie, von P. Scumip, Ziirich*). 

Ahnlich wie in der Thermodynamik kénnen wir in der Lehre iiber neutronen- 
vervielfachende Systeme zwei sich erganzende Auffassungen unterscheiden, nam- 
lich die modellmassig-mechanistische und die phanomenologische Theorie. Wah- 
rend im Zeitalter der Dampfmaschine die Praxis der Theorie oft vorausging, 
wiahrend sich eine phanomenologische Theorie an ihrem Objekt selbst entwickeln 
konnte, steht im Zeitalter der « Atomenergie» die mechanistische Theorie chrono- 
logisch an erster Stelle, und sie ist wohl bis heute auch dominierend. Viele Jahre 
nach dem Bau des ersten Kernreaktors begegnen wir den Anfangen einer sich von 
speziellen Modellen lésenden Auffassung der Reaktortheorie, so zum Beispiel 
E. P. WicNERS Stérungstheorie?), A.M. WerInBERGS allgemeiner Theorie der 
multiplizierenden Medien*), G. GorrtzELs Arbeiten tiber nichtuniforme Brenn- 
stoffverteilung*) und S. M. FEInBERGS heterogenen Methoden zur Reaktorberech- 
nung’). Diese Arbeiten zeichnen sich durch mathematische Eleganz aus, was 
natiirlich damit zusammenhdngt, dass von komplizierten Modellen abstrahiert 
wird. Darin liegt gerade die Kraft solcher Formalismen, und zwar in zweierlei 
Hinsicht ; Einmal lassen sich allgemeingiiltige Gesetzmassigkeiten leichter erken- 
nen und besser formulieren, zum andern kommt man der Wirklichkeit mit allge- 
meinen mathematischen Formen, zusammen mit guten Experimenten, bedeu- 
tend naher als mit primar modellgebundenen Theorien. 

Die vielleicht zu wenig beachtete phanomenologische Reaktortheorie soll nun 
etwas in den Vordergrund geriickt werden. Ihre Aufgabe k6nnen wir wie folgt 
umreissen: Nach der Internationalen Konferenz iiber die friedlichen Anwendungen 
dev Atomenergie steht uns viel experimentelles Erfahrungsmaterial zur Verftigung, 
das schlecht ausgenutzt wird, wenn man es nur als Priifstein fiir diese oder jene 
mechanistische Theorie verwendet. Hier scheint nun vom Standpunkt der Technik 
aus der Einsatz phanomenologischer Theorien lohnender zu sein als die Verfei- 
nerung der Modelle, obwohl das letztere von einem wissenschaftlichen Standpunkt 
aus befriedigender ware. Zuverlassige Experimente sollen also durch mathema- 
tische Formen beschrieben werden. Die Ergebnisse sind auf systematische Zu- 
sammenhange und allgemeine Prinzipien hin zu untersuchen. Im Hinblick auf 
neuartige Konstruktionen sind Experimente vorzuschlagen, welche im Rahmen 
solcher Theorien eine pragnante Formulierung der an die Natur zu stellenden 
Frage verkérpern. Zur Verfolgung dieser Richtung der Reaktortheorie miissen 
gelegentlich bedeutende Rechenhilfsmittel, wie Lochkartenkalkulatoren, elek- 
tronische « digital computers», Analogierechengerate usw., eingesetzt werden. 

Fiir die Schweiz erhebt sich nun die Frage, was auf dem Gebiet der Reaktor- 
theorie in Zukunft getan werden soll und was im Rahmen der Méglichkeiten 
eines kleinen Landes getan werden kann. Nach der persénlichen Auffassung des 
Vortragenden sind folgende Massnahmen ins Auge zu fassen: 


1) Reaktor AG. ; 


*) S. GLAssToNE and M. C. Epiunp, The Elements of Nuclear Reactor Theory (D. van Nostrand 
Co. Inc., New York 1952), p. 372—383. 
3) A. M. WEINBERG, Theory of Oscillating Absorber in a Chain Reactor, AECD-2044, 
A. ue poe Current Status of Nuclear Reactor Theory, Amer. J. Phys. 20, 401 (1952). 
as re igse ken and W. A. Logs, Nonuniform Fuel Distribution in Nuclear Reactors, a 
© GOERTZEL, Minimum Critical Mass and Flat Flux, J. Nuclear En. 2, 193 (1956). 4 
) S. M. FEINBERG, Heterogeneous Methods for Calculating Reactors. Survey of Results and 
Comparison with Experiments, A/Conf. 8/P/669, Genf, August 1955 (Original russisch). 
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4 Auf den Hochschulen werden junge Leute fiir Reaktortheorie interessiert 
and darin zu selbstandigen Arbeiten ermuntert. Diese Disziplin kénnte gewiss 
manchen angewandten Mathematiker und theoretischen Physiker befriedigen 
sind doch bedeutende Wissenschafter, wie E. P. WIGNER oder G. PLACZEK, auth 
nach dem Kriege an diesem und verwandten Gebieten interessiert geblieben. 

2. Der verhaltnismassig bescheidene Stab von Experimentalphysikern und 
Ingenieuren bei der Reaktor AG. sollte durch einige Reaktortheoretiker erganzt 
werden, die zum Teil fiir gréssere Rechenarbeiten im Zusammenhang mit dem 
Versuchsreaktor eingesetzt werden kénnen, zum Teil F achpublikationen verar- 
beiten, gute Experimente analysieren und fortschrittliche theoretische Methoden 
entwickeln. 

Zum Schluss sei noch die Frage aufgeworfen, ob wir nicht heute schon Leute 
zum Studium von schnellen Reaktoren und der gesteuerten Verwendung der 
thermonuklearen Reaktion einsetzen sollten. Jedenfalls liegen Anzeichen dafiir 
vor, dass andernorts solche Probleme zum Teil mit riesigem Aufwand bearbeitet 
werden. 


Gitterberechnungen fiir (D,O—- Unat)-Systeme, von G. DE Couton}), 
R. MEIER’), W. WINKLER®) und W. Ztnt!14), Ziirich. 

Das materielle Buckling B? ist die fiir die Beurteilung eines multiplizierenden 
Mediums am besten geeignete physikalische Grésse. In der fiir thermische Sy- 
steme tiblichen Zweigruppendarstellung wird B? durch die Lésung mit positiver 
Wurzel der Gleichung 
(1+ B’Li) (1+ BLi)=k 
definiert. 

Darin bedeuten L,, die thermische Diffusionslange, L, die Bremslange und k 
den Multiplikationsfaktor der Neutronen im unendlichen Medium'). 

Das Buckling steht in direktem Zusammenhang mit den Abmessungen des 
unreflektierten, kritischen Reaktors. Im Falle der Kugelgeometrie gilt beispiels- ° 
weise B = a/R, wobei R den kritischen Radius des multiplizierenden Mediums 
bedeutet. Im Falle einer allgemeineren Geometrie bleibt die Proportionalitat 
zwischen Lineardimension D und Buckling erhalten: B~ 1/D. 

Die Diffusionsparameter L,, und L, und der Multiplikationsfaktor k lassen 
sich bei gegebener mikroskopischer Struktur aus den Grundkonstanten (Wir- 
kungsquerschnitte, Dichten) der Moderator- und Spaltsubstanzen berechnen, 
womit auch das Buckling bestimmt ist. Dieser rein rechnerische Weg fiihrt hin- 
gegen bei den heutigen Kenntnissen der Grundkonstanten noch nicht auf eine 
Genauigkeit, die fiir eine zuverlassige Dimensionierung eines Reaktors notwendig 
ware. Dies ist aus Tabelle 1 ersichtlich, wo fiir ein quadratisches Gitter von kreis- 
zylindrischen Staben natiirlichen Urans in einem Schwerwassermoderator die 
Fehler 4B2, herriihrend von einzelnen, wesentlichen Grundkonstanten X + AX, 
zusammengestellt sind. Der Fehler im kritischen Volumen ist noch um den 
Faktor 3/2 grésser, ausserdem wurden Fehler der Theorie nicht beriicksichtigt. 
Direkte Messungen des Buckling sind daher vorderhand unerlasslich. 

Fiir (D,O — Unat)-Systeme sind aus Exponentialexperimenten kiirzlich in 
zwei Laboratorien Bucklingwerte fiir eine Reihe von Uranstabradien und Git- 


» 1), 2), 3), 4) in Reaktor AG.; ) auch Gebriider Sulzer AG., Winterthur; und 3) auch Escher 
Wyss AG., Ziirich. 

: 5) S. Grasstone and M. C. Eptunp, The E lements of Nuclear Reactor Theory (D. van Nostrand 
Co., New York 1952). 


ve 


i 
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Tabelle 1 ‘ 

Buckling-Empfindlichkeit (S) auf Fehler von verschiedenen Grundkonstanten (X) fiir ein” 
quadratisches (D,O — Unat)-Gitter mit Stabradius { 


a=1,27 cm und Gitterkonstante b = 12,3 cm 5 


Experimentelle Resultierende 
Fehler S Fehler 
AX/X in % AB?) B2 ime, 


1:3ipet 0,025.43 
y 2,47 + 0,10 
Chass BOO RE 45.be Tl .ae 2,6 2,6 6,8 
Cait LOT: HAS. bie tenes 4,7 — 0,62 eG 
OTys20 2,76+0,04b... 1,5 O44 ENS4 
%,0 2 4 02 be: 11 “— 0,08 2.99 
A 2,50: 6,05 cm 4 2,0 0, 50iaon 1,0 
D,O , 
2 } 
Lap.0 ....10200 500 ent? 5 10,50 25. 
Linio 125 brome? -dgaqi! 4,8 0,50 2,4 


9416 


> 
terkonstanten ermittelt worden. In einem schwedischen Versuch (ZEBRAS) 
wurde ein Tank von 1m Durchmesser und eine Ra-Be-Quelle verwendet. Ein 
amerikanisches Experiment (NAA’)) bestand aus einem Tank von 1,52 m Durch- 
messer und einem Water-Boiler-Reaktor als Neutronenquelle. j 
Diese Resultate gaben uns Gelegenheit, die Gittertheorie zu priifen und die 
mit Hilfe direkter Messungen nicht mit geniigender Genauigkeit bestimmbaren 
Grundkonstanten anzupassen. Eine von MumMaRy§) erstmals fiir eraphierien al 
rierte Systeme angewendete Methode ist dazu besonders geeignet, weil sie die 
in der Rechnung unsichersten Werte von 7 und dem Resonanzintegral als Para- 
meter frei lasst, wahrenddem L,, Ly, f und e, die weniger empfindliche Ts 
enthalten, berechnet werden. 


Dem effektiven Resonanzintegral legen wir die experimentell gut oe 
Form 


(/oete = (8+ ee) | 


mit den Bezeichnungen S fiir die Oberflache one M fiir die Masse pro Zentimete 
Uranstab zugrunde. 


6) R. Persson, E. Bromsj6, M. Bustraan and R. MEIER (wird demnachst im Journal of 
Nuclear Energy publiziert). 


a K, Conen, Geneva Paper No. P/605 (1955). 
8) W. Mummary, Geneva Paper No. P/429. 
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Die Analyse der Versuchsdaten fiihrt dann auf die Werte von Tabelle 2. 


Tabelle 2 


A lu 


ZEBRA‘) 
6,4 NAA’) 


Mit den Mittelwerten der Daten von Tabelle 2 wurde die Bucklingflache in 
Abhangigkeit vom Stabradius a und der Gitterkonstanten 6 berechnet. Das 
Relief dieser Flache ist in Figur 1 dargestellt. 

Figur 2 zeigt einen Schnitt durch die berechnete Bucklingflache mit dem 
Stabradius a= 1,27 cm. Vergleichsweise sind die von R. PERSSON e¢ al.6) und 
K. COHEN’) gemessenen Werte aufgetragen. Die experimentellen Punkte konnen 
innerhalb der Messgenauigkeit angepasst werden. 


Figur 1 ~ 
‘Buckling-Flache (Parameter B? m~®) fiir quadratisches Gitter (b cm) kreiszylindrischer Uranstibe 
(Radius a cm) in schwerem Wasser (99,76%)- Aluminiumhiille 0,1 cm dick. 
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Figur 2 


Buckling B? fiir (D,O — Unat)-Gitter, Stabradius a = 1,27 cm. ——— berechnet mit y = 1,302 
ue Ai he 
und ( f o(E£) dE|E) eff = 5,6 + 32 S/M; a NAA-Exponentialexperiment; ° ZEBRA-Exponentiai-— 


experiment. 


Der thermische Ausgleichsvorgang in einem Medium mit internen 
Warmequellen, von K. HINTERMANN, Ziirich!). 


Einfihrung 


In einem Reaktor bestehen in allen Teilen kontinuierlich verteilte innere 
Warmequelien. Im stationdéren Fall miissen die dadurch hervorgerufenen Tem- 
peraturfelder fiir die Konstruktion der Kiihlsysteme erfasst werden k6nnen. 
Wegen des Einflusses der Temperatur auf die Reaktivitaét eines Reaktors miissen, © 
im Hinblick auf die Regulierung, die Temperaturen im Reaktorkern jedoch auch — 
bei zeitlichen Anderungen betrachtet werden. Ein wesentlicher Teil dieses Pro- 
blems ist der thermische Ausgleichsvorgang in den Brennstoffelementen selbst, 
in welchen tiber 90% der gesamten Warme entstehen. 


Bevrechnung : . 

Die Warmeleitungsgleichung mit einem Quellenglied lautet: ; | 
oT 

AT — 22S —=O6¢-% 

a - (7, t) (1) 


1) Reaktor AG., Ziirich, und AG. Brown, Boveri & Cie., Baden. 
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oder eindimensional fiir eine Platte, nach entsprechender Wahl der Einheiten: 


oT oT “aml 
“aed ore) toe (2) 


Hier erscheinen die Langen mit der Dimension //Zeit (Einheit: //s). Die Einheiten 
stellen nichts anderes dar, als die Langen in Zentimeter, dividiert durch die 
Quadratwurzel aus der Temperaturleitzahl a. 

Die Lésung der Gleichung (2) ist mittels Laplace-Transformation?) méglich. 
Falls der Nullpunkt der Temperaturskala so gewahlt wird, dass er mit der Rand- 
temperatur zusammenfallt, lautet die Lésung: 


l 1 t 
T(x, t) = [Tw &;t) T(E, 0) dé + fas fre, £;t— 1) OE, 2) dt, (3) 
0 0 Gas 


ee : 
I(x, €; 2) = wr De ” sinn + x sinn + & 


(I Plattendicke in //s). 


Trotzdem am Rande der Heizelemente keine konstante Temperatur herrscht, 
sondern ein Warmeiibergang an einen Warmetrager stattfindet, st6rt hier die 
Voraussetzung konstanter Randtemperaturen nicht, da in jedem Falle der Warme- 
iibergang durch Warmeleitung ersetzt werden kann. Es braucht dazu nur ein 
quellenfreies Gebiet ausserhalb des Randes in die Lésung einbezogen zu werden. 


Diskussion der Lésung 


Zur Diskussion sei die Lésung fiir raumlich konstante, zur Zeit ¢= 0 mit 


_konstanter Starke einsetzende Warmeproduktion explizit angeschrieben: 


1 
eZ oe —n tll _- 4 - A 
Dita by D>? sinn | xf T(é, 0) sinn sa & dé 


n=1 
: (4) 
4]? _ JI —n? x3 t/l? ; 4 
= aa @ ps ere he ) sinn = *. 
n=2m—-1 


Das erste Glied stellt zur Zeit t = 0 die Fourier-Reihe fiir den Anfangstem- 


‘peraturverlauf T(x, 0) dar und wird mit ¢ > oo gegen 0 abgebaut, und zwar in 


erster Naherung mit der Zeitkonstanten /?/z?. 
Das zweite Glied beginnt bei ¢= 0 mit 0 und geht fiir ¢ > oo mit derselben 


-Zeitkonstanten [2/m® in diejenige Fourier-Reihe iiber, die denjenigen Parabel- 


bogen darstellt, den man auch aus der stationaren Warmeleitungsgleichung er- 


halt, wie es auch sein muss. 


Ausgleichszeiten, Z ylindergeometrie 


Die Zeitkonstante des Ausgleichsvorganges (« Ausgleichszeit ») bei einer Platte 
betragt 72/2? fiir konstante Randtemperaturen. . 


2) G. Dorrscu, Theorie und Anwendung der Laplace-Transformation (Springer, Berlin 1937). 


eee yh Oe ee 
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Setzen wir statt J in /s din Zentimeter ein, so ist: 


i a d*co 


(a Temperaturleitzahl, c spezifische Warme, @ Dichte, k Warmeleitkoeffizient). 


‘ 


Die Ausgleichszeiten sind also proportional dem Quadrat der linearen Ausdeh-— 
nung eines Kérpers und umgekehrt proportional der Temperaturleitzahl. Sie sind 


unabhangig von der Quellstarke. 


Der Zahler des Ausdruckes fiir die Ausgleichszeit (5) ist proportional dem Ver-— 
haltnis Warmeinhalt zu erzeugter Warmemenge. Dass dieser «spezifische » Warme-— 
inhalt proportional zu d? ist, sieht man ein, wenn man den stationaren Fall 
betrachtet: Ausser der MaBstabsvergrésserung mit zunehmendem d miissen noch ~ 


die Temperaturgefalle linear zunehmen, um die mehr erzeugte Warme abzufuhren. 


25 


Nach oo § 
(sfationdr) 


Oo QUIS ORF a OF ieee ORD GO fe £07 
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Temperaturverlauf in einer ebenen Platte nach Einsetzen interner Warmequellen. T Temperatu- / 


ren; Ty Randtemperatur; Ty.» Temperatur in Plattenmitte nach cos; d?/a = 57,38; d Platten- 
dicke; a Temperaturleitzahl. 


Die Flache unter der Temperaturparabel ist also proportional zu d%, wahrend die 
durch die Quellenverteilung erzeugte Warme linear mit d verlauit. 

Wird nun die Plattendicke als Durchmesser in die Zylindergeometrie iiber- 
gefiihrt, so zeigt sich, dass die Parabel des stationdren Zustandes eine Parabel 
bleibt, jedoch ihre Scheitelhdhe nur mehr die Halfte derjenigen der Plattengeo- 
metrie betragt (Umfang zu Flache wird doppelt). Ferner ist die mittlere Tempe- 


} 
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ratur einer parabolischen Verteilung im Zylinder 1/2 der Scheiteltemperatur 
.wahrend Sie in einer Platte 2/3 betragt. Der Warmeinhalt in Zylindergeometrie 
und damit auch die Ausgleichszeiten betragen also 


4) gee re 3 
a a Ht ee tpl: 


Resultate 


Figur 1 zeigt den Temperaturverlauf in einer Platte mit Anfangstemperatur 
T = 0 bei raumlich konstanter, zur Zeit ¢ = 0 mit konstanter Starke einsetzender 
Quellenverteilung zu verschiedenen Zeiten. Man sieht deutlich, wie sich zu 
Beginn der Randeinfluss noch nicht voll bis zur Mitte auswirken kann und daher 
die Kurven anfanglich etwas abgeplattet sind. 


0 5 
0 Zz 4 6 8 70 Platengeomerrieé 
0 ] 2 3 Zylindergeomeirie 

Figur 2 


Zeitlicher Verlauf von Temperaturen in einem Medium mit internen Warmequellen. T Tempera- 

_turen; Ty) Randtemperatur; T,. Temperatur fiir >00; ¢ Zeit seit Einsetzen der Quellen. 1 mittlere 

Temperatur, gleichmdssige Verhdltnisse; 2 mittlere Temperatur, gemittelt iiber Plattengeometrie; 

3 zeitlicher Verlauf des Temperaturmaximums. d?/a = 57,38; d Plattendicke, Zylinderdurch- 
messer; a Temperaturleitzahl [cm? sp: 


Die angegebenen Zeiten beziehen sich auf eine Platte von 2,6 cm Uran. Die 
_Kurven werden jedoch allgemein giiltig, wenn man die Zeiten statt in s in Ein- 
_heiten von d?/a angibt. 
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(Zeitkonstante /?/12). 

Wenn man also bei variierenden Verhaltnissen mit den gleichbleibenden Be- 
dingungen, wie sie der Kurve 1 zugrunde liegen, rechnet, so ist der grésstmégliche 
Fehler (bei idealem Warmeiibergang am Rand) 18%. Mit endlichen Warmetiber- 
gangszahlen am Rand wird der Fehler kleiner. 

. 
Untersuchung iiber Vakuumiiberschlage an Beschleunigungsrohren, von 
F. BuMILLER, Ziirich!), } 

Die Isolationsfestigkeit des Beschleunigungsrohres des Ziircher Bandgenera-_ 
tors?) sank im Verlauf von etwa 300 Betriebsstunden von 1,8 MV auf 1,2 MV. — 
‘Nach der Zerlegung des Rohres liess sich auf den Isolatoren eine starke Verkoh-: 
lung feststellen, die mit ziemlicher Sicherheit auf das zum Kleben verwendete 
Araldit zuriickzufiihren war, gegen dessen Eindringen in den Vakuumraum 


keine besonderen Vorsichtsmassregeln beobachtet wurden (Figur 1). 7 
; 

Araldit e 
Stah/ ‘ 

; 

i 

: 

' 

#s 

oF 

Porzellan r 
glasiert i 
4 


: 
} 


Figur 1 : 

Auf Grund von neueren Erfahrungen, wie sie insbesondere von der Van de 
Graaff-Gruppe in Saclay (Paris)*) gesammelt wurden, wurde ein neues Beschleu- 
nigungsrohr mit glatten, unglasierten Porzellanisolatoren und Aluminiumelek-§ 
troden gebaut. Obwohl dieses Rohr zu Beginn 1,9 MV bis 1,8 MV isolierte, sani 
die Spannungsfestigkeit im Verlauf einer Woche auf 1,2 MV. Die einzelnen Iso- 
latoren zeigten kraftige Aluminiumbahnen, deren Widerstand zum Teil unter 
1 MQ lag. Besonders aufschlussreich waren Aluminiumfusspunkte, welche an den © 
Stellen am Isolator entstanden, wo der Elektrodenring das Porzellan beriihrte. . 
Sie sind leicht als der Beginn der Aluminiumbahnen zu erkennen (Figur 2). 

Um den Bandgenerator betriebsbereit zu halten, wurde aus den gereinigten 
Isolatoren des ersten und den Elektroden des zweiten Rohres ein Ubergangsrohr 
gebaut. Dabei wurden Erfahrungen, die aus der Zerst6rung der ersten Rohre 


1) Physikalisches Institut der Universitat. 
2) F, Bumiier, V. Meyer, H. H. Straus und H. 


WINKLER, Z. Naturf. 10 i : 
8) Rapport C. E. A. No. 189, 1953 (Centre d’ uae 


Etudes Nucléaires Saclay, Paris). 
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Silberanstrich 


Aluminium 


Figur 2 


gezogen werden konnten, bereits beriicksichtigt. So wurde durch Versilbern der 
Isolator-Metalltrennflache und Einlegen einer Aluminiumfolie ein besserer Kon- 
takt ermdéglicht; durch eine Nute im Zwischenblech wurde das Vordringen des 
Araldits in den Vakuumraum verhindert (Figur 3). Nach einer anfanglichen 
Isolationsfestigkeit von 1,8 MV halt dieses Rohr nach etwa 500 Betriebsstunden 
immerhin noch 1,5 MV. 

Al-folie i 
Aluminium 
(von) 


giasiert 777 7 
(ORL) oe Dae 7 


In der Zwischenzeit wurden Versuche fiir die beste Isolator-Elektroden-An- 
ordnung durchgefiihrt. Uber eine Giessharzausfiihrung wurde eine Spannung von 
maximal 250 kV vom Bandgenerator an das Versuchsgefass angelegt. Zur Beur- 
teilung der Messungen wurde die Beziehung von CRANBERG‘) 


V2 
VE = Coder ee = C 
fiir Parallelanordnung herangezogen. V ist die angelegte Spannung, E die Feld- 
starke, d der Elektrodenabstand. C ist eine vom Material abhangige Konstante 
(iiberschreitet VE ihren Wert, so erfolgt Durchbruch). Der Wert von C betragt 


fiir Elektrodendurchschlage etwa 
j 7 (MV)? 


, 


4) L. CRANBERG, J. appl. Phys. 23, 518 (1952). 


é 
.. 


> 
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fiir Isolatoriiberschlage liegt er allgemein zwischen 0,1 und 0,3 (MV)?/m und kann | 
fiir sehr gute Isolatoren auf rund 0,8 (MV)?/m ansteigen. Es wurden folgende 
Messungen durchgefiihrt: “ ; % 

1. Die Festigkeit gegen Elektrodendurchschlag wurde gepruft. Sie lag iiber 
der von der Anlage erreichten Spannung von 250kV. Dies entspricht einem 
Wert von C w» 7 (MV)?/m. Bae 

2. Es wurde vermutet, dass die Aluminiumfusspunkte, wie sie im Rohr II 
auftraten, vom Kontakt des Elektrodenringes mit dem Isolator herriihrten. Um 
eine gréssere Zahl von Uberschlagen zu erhalten, wurde eine Anordnung mit ver- 
kiirztem Ringabstand aufgebaut, wobei im ersten Modell die Ringe satt am Por- 
zellan lagen, im zweiten der Abstand jedoch 2 mm betrug. 


Ag eingebrannt 


a Sahl 


= " verchromt 
guarzsandgestrahit 


Neves Rohr 


Figur 4 


Es gelang mit dem ersten Aufbau wiederum Fusspunkte und beginnende | 


Bahnen zu erzeugen. Die Isolationsfestigkeit betrug 120 kV bei einer Funkenfolge ~ 
von etwa 1/s. ; 

Die zweite Anordnung hielt 170 kV bei einer Funkenfolge von etwa 1/min. — 
Es waren keinerlei Spuren am Isolator festzustellen. 


3. Auf Grund dieses giinstigen Ergebnisses wurde eine Elektrodenanordnung 


aus verchromtem Stahl angefertigt, wie sie aus Figur 4 hervorgeht. Mit diesen 
Elektroden wurden folgende Isolatoren ausprobiert: J 


Isolator Viibersepltue . (6; 
a) Glatter Isolator, quarzsandgestrahlt .. . 160 kV 0,43 : 
b), Gewellter Isolator, plasiert) “2.2.4. 2. 190 kV 0,60 (0,47) - 
c) Gewellter Isolator, quarzsandgestrahlt . . 215kV 0,77 (0,61) 


Die C-Werte in Klammern wiirden einer Beriicksichtigung der Wellung ent- : 
sprechen. Mikroskopische Untersuchung von Funkenaufnahmen ergab jedoch, © 
dass die meisten Uberschlage geraden Bahnen folgen. Die gewellte Oberflache 
scheint eher eine schnelle Verschmelzung zu verhindern. Nach diesen Messungen 


und einigen Zwischenuntersuchungen diirften fiir ein Beschleunigungsrohr fol- 
gende Bauregeln zu empfehlen sein: 


ee 
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Porzellan ist, zumindest im Vakuum, der beste [solator, wobei eine rauhe, 
unglasierte Oberflache vorzuziehen ist. Trotz héherer Sekundaremission siud 
hoherverdampfende Schwermetalle den Aluminiumelektroden vorzuziehen. Das 
Eindringen von organischen Materialien in den Vakuumiraum ist unbedingt zu 
vermeiden. 

Absolute Sauberkeit der Isolatoren-(Fingerabdriicke!) und gute Politur der 
Elektroden haben einen wesentlichen Einfluss auf die Isolationsfestigkeit. 

Fur guten Kontakt zwischen Metall und Porzellan ist zu sorgen. Freie Isola- 
toren- und Seitenflachen sollen sich unter rechtem Winkel schneiden. 

Nach diesen Erfahrungen wurde ein neues Beschleunigungsrohr vom Typ 3c 
gebaut, welches in den nachsten Wochen in Betrieb genommen wird. 


Knickung von Schraubenfedern unter Druck und konservativer Torsion, 
von Hex. Maaa, Ziirich!). 

Es wurde das Knickproblem einer Schraubenfeder untersucht, die durch eine 
axiale Druckkraft und ein semi-, quasi- oder pseudotangentiales Torsionsmoment 
belastet ist. Diese Momente ergeben sich nach H. ZIEGLER?) aus einfachen Krafte- 
verteilungen, wahrend das iibliche axiale Torsionsmoment, mit dem vielfach 
gerechnet wird, sich ausser an Einspannungen praktisch kaum realisieren lasst. 
Dazu sind diese Momente — wieder im Gegensatz zum axialen Torsionsmoment — ~ 
konservatiy, so dass man die Stabilitatsuntersuchung rein statisch durchfiihren 
darf; also die kritischen Lasten auf Grund des erstmaligen Auftretens von Gleich- 
gewicht bei einer ausgebogenen elastischen Linie ermitteln kann. 

Da sich eine Schraubenfeder mit geringer Steigung elastisch als weicher Stab 
behandeln lasst, stellt diese Arbeit eine Erweiterung der Untersuchung von 
M. Becxk®) iiber das Knickproblem bei schlanken Staben unter Druck und kon- 
servativer Torsion dar. Die Erweiterung besteht darin, dass im Gegensatz zum 
schlanken Stab bei der Feder auch der Langenanderung und vor allem der 
Schiebung Rechnung getragen werden muss. Mit der fiir solche Eigenwertpro- 
bleme typischen Methode werden die Differentialgleichungen der elastischen 
Linie unter Beriicksichtigung der Randbedingungen fiir verschiedene Lagerungen 
der Feder gelést und die Knickgleichungen als Bedingungen fiir nichtidentisches 
Verschwinden dieser Lésungen berechnet. Diese Knickgleichungen stellen im all- 
gemeinen transzendente Beziehungen zwischen der kritischen Druckkraft, dem 
kritischen Torsionsmoment und Form- und Materialkonstanten der Feder dar. 
In der Mehrzahl der untersuchten Falle wurden diese Gleichungen in Diagrammen 
graphisch dargestellt und der Einfluss des Schlankheitsgrades der F eder und der 
Querschnittsform des Federdrahtes auf die kritischen Lasten aufgezeigt. Eine aus- 
fiihrliche Veréffentlichung ist in Vorbereitung. 


Uber die Anwendung des elektrisch leitenden «Teledeltos»-Papiers fur 
die Ausmessung von Potentialfeldern, von E. B. Bas, E, BERNHARD und 
M. NIEDERHAUSER, Ziirich?). 

- Die Arbeit behandelt den Einfluss des anisotropischen Charakters der elek- 
trischen Leitfahigkeit von « Teledeltos »-Papier L 39 bei der Anwendung dieses 


1) Assistenz fiir technische Mechanik an der ETH. 

2) H. Zrecter, ZAMP 3, 96 (1952). 

3) M. Beck, Ing.-Arch. 23, 231 (1955). 
1) ETH, Abteilung fiir industrielle Forschung des Institutes fiir technische Physik. 
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Papiers fiir Potentialfeldmessungen. Es wird gezeigt, dass die elektrische Leit- 
fahigkeit in Walzrichtung des Papiers um rund 6% grodsser ist als quer dazu. 
Diese Leitfahigkeitsanisotropie kann durch eine affine Transformation beritick- 
sichtigt werden. An zwei theoretisch errechenbaren Beispielen wird der Einfluss 
der Anisotropie studiert. Das erste Beispiel behandelt die Potentialverteilung 
entlang der Achse einer elektronenoptischen zweidimensionalen Zweielektroden- 
Kathodenlinse. Bei Nichtberiicksichtigung der Anisotropie wurde eine maximale 
Abweichung des gemessenen Potentialverlaufes von dem gerechneten Verlauf 
von etwa 2% gefunden (bezogen auf die Elektrodenspannung). Dieser Fehler hat 
entlang der Achse denselben Verlauf wie die errechnete Funktion 0°D/d%*. Be- 
riicksichtigt man die Anisotropie, so reduziert sich die maximale Abweichung 
auf rund 5/9). Das zweite Beispiel behandelt die Potentialverteilung in einem 
hyperbolischen Quadropol: y = u (¥?— y?). Es wird gezeigt, dass hier das Feld 
in der Nahe des Sattelpunktes auf die Leitfahigkeitsanisotropie sehr empfindlich 
ist und dass auch hier die Beriicksichtigung der Anisotropie durch affine Koordi- 
natentransformation weitgehende Verbesserungen erwirkt. Aus diesen Untersu- 
chungen kann das Resultat gezogen werden, dass man besonders bei der Ausmes- 
sung von stark inhomogenen Feldern mit Hilfe des « Teledeltos »-Papiers vorsichtig 
sein soll. Es ist empfehlenswert, die Anisotropie stets zu korrigieren, und zwar 
mit einem Anisotropieverhaltnis des gerade zur Messung verwendeten Papier- 
stiickes. Zur leichten Bestimmung des Anisotropieverhdltnisses (o;/o;;) wird in 
Anlehnung an die erwahnten Messungen am Quadropol ein kleines Gerat vorge- 
schlagen, welches einen Quadropol aus zwei gekreuzten Plattenkondensatoren als 
Elektroden darstellt. Eine der vier Platten wird mit Hilfe einer Mikrometer- 
schraube verschiebbar angeordnet. Das Mittelpunktpotential des auf diese Weise 
affin verzerrbaren Quadropols wird gemessen und durch Verschiebung der be- 
weglichen Platte auf das richtige Sattelpunktpotential eingestellt. Die bendtigte 
affine Verzerrung des Modells kann dann direkt vom Mikrometer abgelesen 
werden, Mit einem solchen Gerat kann das verwendete Papier in allen Bezirken 
auf die Anisotropie kontrolliert und ein beziiglich Anisotropie homogenes Stiick 
ausgesucht werden. Ein ausfiihrlicher Bericht erscheint demnachst in der ZAMP. 


Anwendung von kolloidalen radioaktiven Substanzen in der Medizin. 
Messtechnik. Von G. Poretti und E. Wi1LpBo1z, Bern). . 
Durch Einspritzung von radioaktiven Substanzen in das Tumorgewebe kann 
eine einwandfreie therapeutische Wirkung erreicht werden, wenn: 
_1. die Tumoren von aussen oder chirurgisch erreichbar sind und ihr Volumen 
mit guter N&herung bestimmt werden kann; : 
2. die eingespritzte Substanz physiologisch neutral ist und lange Zeit an Ort 
und Stelle bleibt. 


_ Ist das Tumorvolumen bekannt und wird die Bedingung 2 erfillt, so kann 
die Strahlendosis bestimmt werden. q 
Schon seit einigen Jahren ist es bekannt, dass die kolloidalen Substanzen die 
Bedingung 2 mehr oder weniger erfiillen. Es wurden &8ZnS, 52MnO,, 18Au_ 
metallisch, **FeSO,, Cr?2PO,?) und ZrOH*?PO, angewandt. Wir haben haupt- 
sachlich die therapeutische Wirkung des 22P i yl- 
: peutische Wirkung ‘des *P im Chromphosphat und Zirkonyl-_ 

phosphat bei Blasenkarzinom erprobt und dabei versucht, fiir das CrPO, die Giil- 


tigkeit des Punktes 2 am Menschen zu beweisen. CrPO, wurde statt mit °2P mit 

1) Inselspital. 

*) M. E. Morton, Nucleonics 10, No. 11, 92 (1952). . 
2 
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51Crmarkiert, welches mit hoher spezifischer AktivitAt [(n, y)-R Bo 
aus Harwell erhaltlich ist. §1Cr zerfallt nach pete cebeeienimesiial: 
beute etwa 26%) und emittiert eine Gammalinie (3%) von 0 323 MeV, welche mit 
Hilfe eines Nal(T1)-Szintillationszahlers gut messbar ist. 

Die in unseren Versuchen von den eingespritzten 1,2 mc Chrom gelieferte 
Dosis wurde in 100 rep total errechnet. Der Geometriefaktor (Kreiszylinder) fiir 
die Gammadosis wurde an Hand der Berechnung von PLEscH’) bestimmt. 


Figur 1 


Die Chrommarkierung erlaubt, das an der Harnblasenwand angelagerte 
‘CrPO, von aussen nachzuweisen. Die Auffindung der aktiven Stelle erfolgte mit 
Hilfe eines Scinti-Scanners*) (Figur 1), einer an der Universitat Bern konstruierten 
Apparatur, welche die automatische Abtastung von Korperteilen mittels eines 
‘Szintillationszahlers erlaubt. Die aktiven Stellen oder Organe werden graphisch 
registriert (Schilddriise: 1°11; Leber: 1®*Au usw.). 

Die Stelle, an der sich das CrPO, angereichert hat, ist in Figur 2 ersichtlich, 
‘und daraus entnimmt man, dass die Substanz effektiv 27 Tage lang (Chromhalb- 
-wertzeit) und mehr in loco geblieben ist. Eine analoge Messung an einem zweiten 


‘Patienten lieferte gleiche Resultate. 


3) R. Prescn, Strahlentherapie 97, 277 (1955). 
4) W. Botiicer und G. PorertTI, Experientia 12, 116 (1956). 


Sees. 


372 Varia — Miscellaneous — Divers ZAMP 


Patient S.A. ’91 (14. 12. 1955) ee 
= emp ea FPR eli hot fa a 
| i Len Wve Peo thei Akay 
Pe yt 
Lill 


Intra dt a ; 
aT f eg 


iu 
| 
| 


Durch diese Ergebnisse wird bestatigt, dass 

1. die effektive Halbwertzeit des *2P in der mit Cr®*PO, infiltrierten Stelle 
14 Tage betragt; 

2. die von den 32P-f-Strahlen in dem infiltrierten Volumen abgegebene Energie 
(Dosis) mit Hilfe der Formel?) ; 


Dee z. Ty, C [rads] 
(E, = mittlere Energie des f- Bed in MeV, T,,, = effektive Halbwertzeitl 
in Tagen, C = *2P-Konzentration in Mikrocurie pro Gramm Gewebe) é 


errechnet werden kann. 

Fiir eine Einspritzung von 4 mc 32P in 10 g Tumorgewebe betragt deshalb die 
Dosis 290 000 rads. Die von der 1,7 MeV Betastrahlung im Gewebe befreite Brems-— 5 
strahlungsenergie betragt weniger als 1% der Ionisationsenergie®). 


Experimentelle Untersuchungen tiber den Einfluss der Eintrittsgrenz-_ 
schicht auf die Druckumsetzung in geraden und gekrimmten Diffusoren’), 
von H. SPRENGER, Ziirich?). 
Diffusoren sind Kanale mit zunehmendem Strémungsquerschnitt. Sie dienen 
bei tropfbaren Fliissigkeiten und bei Gasen (im Unterschallbereich) zur teilweisen 
Umsetzung von kinetischer Energie in statischen Druck. Die Energieumsetzung 
in Diffusoren wird von zahlreichen Parametern beeinflusst, beispielsweise durch 
das Flachenverhaltnis, den Erweiterungswinkel, die Kanalkriimmung, die 
Unterschiede der Form aufeinanderfolgender Querschnitte, durch den Turbulenz- | 
grad des Str6mungsmediums, sowie in besonders hohem Masse durch die Dicke 
und Symmetrie der Reibungsgrenzschicht beim Eintritt in die divergenten Ka- 
nale. An 14 verschiedenen Diffusoren mit D, = 10 cm Eintrittsdurchmesser und 
dem maximalen Flachenverhaltnis 1:4 wurde versucht, experimentell einen Teil 
dieser Abhangigkeiten zu ermitteln. Die Untersuchungen sind mit Luft im in- 


q 
’) L. D. MarrNexii, Radiology 62, 656 (1954). 
8) S. J. Warp, Meleouics TGR No. 7, 44 (1955). 
1) Ein ausfiihrlicher Bericht erscheint in den « Mitteilungen aus dem Institut fiir Aerodynamik 
an der ETH », herausgegeben von Prof. Dr. J. ACKERET, 4 
2) Tstitut fiir Aerodynamik der ETH. : 

‘ 

. 


Vol. VII, 1956 Varia — Miscellaneous — Divers 373 


kompressiblen Bereich ausgefiihrt worden bei mittleren Eintrittsgeschwindig- 
keiten %, = 600 bis 12000 cm/s; diese ergeben Reynolds-Zahlen (bezogen auf D,) 
zwischen 30000 und 800000. Zur Beurteilung der Reibungsgrenzschicht dient 
wie iiblich die Verdrangungsdicke 6*; sie konnte durch Vorschalten zylindrischer 
Rohre variiert werden zwischen 0,3 und 6%, bezogen auf den Rohrradius. 

Die Messungen brachten im wesentlichen folgende Ergebnisse: 

Bei geraden 8°-Diffusoren und diinnen turbulenten Grenzschichten am Eintritt 
wird ein Druckanstieg von 50% des Eintrittsstaudrucks praktisch verlustlos er- 
reicht. Gekriimmte Diffusoren mit gleichen Eintrittsbedingungen haben allge- 
mein schlechtere Wirkungsgrade, kénnen aber durch strémungstechnische Mass- 
nahmen, wie Abplattung in der Biegeebene, etwas verbessert werden; Flachen- 
verhaltnisse tiber 1:3 bringen nur noch bescheidene Gewinne. Der Einfluss dicker 
Grenzschichten erwies sich als sehr erheblich, der Druckriickgewinn sinkt bei 


257De = 1,15 Yo 
Re =520 000 


L/De=4,85 


287/De=1,75% 
Re = 520000 


L/De=6.85 


2. 5°/De=2,3% 
Re = 520000 


¥ = 31% { 
shee - \ = 0,803) 


Finfluss der Eintrittsgrenzschicht (Verdrangungsdicke 6*) auf den Druckriickgewinn (7) fir den 
geraden Kreiskegeldiffusor (F,:Fg = 1:4) bei konstanter Reynoldscher Zahl (Re). Die Aus- 
trittsenergieverteilungen sind mit einem relativ tragen Registriermanometer aufgenommen ; die 
Schwankungen vermitteln nur ein qualitatives Bild der Turbulenz. (Die beiden Definitionen von 
n = (bq — Pe)l(Ie — Iq) Unterscheiden sich durch die Art der Mittelwertbildung der Staudrucke. 
Es ist q = 0/2 u? der Staudruck der mittleren Geschwindigkeit, wahrend drei Striche one 
dass durch Integration der Verteilung von u® ein Staudruck bestimmt wurde, der dem Wert der 
. effektiv eintretenden Energie entspricht.) 
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geraden Diffusoren um etwa 20% und bei Hochkantkriimmern bis 40%. Die 
Resultate vermitteln Einblick in die Stromungsverhaltnisse bei technischen Difs® 
fusoren, auf deren Bedeutung fiir den Wirkungsgrad von Turbokompressoren — 
und Niederdruckturbinen abschliessend hingewiesen sei. « 
Das Bild zeigt instruktiv einige Ergebnisse der Messungen am geraden Kreis-_ 
kegeldiffusor. Nach einem neuen, von Herrn Prof. Dr. J. ACKERET vorgeschla- “ 
genen Rechnungsgang gelang es, den Verlauf der Grenzschicht und somit auch 
die Druckumsetzung dieses Diffusors mathematisch zu erfassen, wenn gewisse 
empirisch ermittelte Konstanten mitverwendet werden. Das Verfahren beruht 
auf dem Einbau der Kontinuitatsgleichung in differentieller Form in die Differen- 
tialgleichung der Grenzschicht (Impuls- oder Energiesatz), wodurch Iterationen 
wegfallen und gute Konvergenz erreicht wird. Es ist anzunehmen, dass die mit- 
benutzten, empirisch gefundenen Werte sich auf Grund noch auszufiihrender 
Turbulenzmessungen im Diffusorinnern auch physikalisch interpretieren lassen. 


= so 
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Kleines Handbuch technischer Regelvorgange. Von WINFRIED OPPELT 
(Verlag Chemie, Weinheim 1954). 448 S., 353 Fig., 110 Taf. und 3 Tab.; Ganzleinen” 
DM 31.80. P i 
Dieses «kleine» Handbuch ist aus friiheren Biichern des Verfassers entstanden 
und hat einen ganz betrachtlichen Umfang angenommen. Es fiillt eine Liicke im 
deutschen Schrifttum und erganzt, ja tibertrifft die ausgezeichneten amerikani- 
schen Biicher der Servotechnik von CHESTNUT and MAYER, NIXON usw. Die Rege- 
_ lungstechnik, die unter dem Namen «Servomechanisms and regulating systems» 
wahrend und nach dem Zweiten Weltkrieg in den USA. einen enormen Auf-_ 
schwung genommen hat und nun auch in Europa iiberall Eingang findet, kann 
einerseits in die allen Systemen gemeinsame «Regelungstheorie» und anderseits- 
in die spezielle Kenntnis des Aufbaues der Servomechanismen, in die «Gerdte- 
technik», aufgeteilt werden. Um eine Vorstellung von den heutigen Méglichkeiten — 
zu geben, behandelt das Buch aus der Geratetechnik mechanische, hydraulische, 
pneumatische, kalorische und elektrische Beispiele. Eine systematische Einrei-_ 
hung erfahren die Messwerke (Geber, Fiihler) sowie die Verstarker, darunter tie. 
magnetischen Verstarker, die gesteuerten Stromrichter und die Maschinen yeas 
starker (zum Beispiel die Amplidynmaschinen). > 
___ Der Hauptteil des Buches befasst sich mit der Regelungslehre. Die komplexen 
Ubertragungsfunktionen werden mittels ungedampfter sinusformiger Variatio-_ 
nen, unter Vermeidung der Laplace-Transformation, erklart und durch komplexe 
Frequenzgangkurven in der GauBschen Ebene dargestellt. Die von BLocH und 
MICHAEL entwickelte Ortskurvengeometrie in der komplexen Zahlenebene ve 


Phasenganges nach BopE wenig gepflegt. Eingehend besprochen werden die Sta- 
bilitatskriterien nach Roura-Hurwirz, insbesondere die Stabilitatsgrenzen be: 
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Differentialgleichungen dritter und vierter Ordnung, ferner die Kriterien von 
CREMER-LEONHARD und Nyquist und das besondere Verhalten der Proportional- 
und Integralregler sowie ihrer Kombinationen. Die nichtlinearen Regelkreise, 
die Kreise mit Reibung, Anschlagen und Spiel («Lose»), die «Ein-Aus»-Regler, 
die Modellregelkreise und die Blockschaltbilder von elektronischen Analogie- 
rechengeraten beschliessen das Buch..In-sprachlicher Hinsicht enthalt es ver- 
schiedene Neuworter. Der Text ist klar gefasst, die Ausstattung ist einwandfrei. 
Die besonderen Kennzeichen des Buches sind die Griindlichkeit sowie die enorme 
Mannigfaltigkeit der Figuren, Tabellen und der im Text verarbeiteten Literatur- 
stellen. 

Im Friihjahr 1956 erschien bereits eine zweite, neu bearbeitete und erwei- 
terte Auflage mit 555 Seiten, 471 Abbildungen, 126 Tafeln und 3 Tabellen. 
Neu aufgenommen wurde die Wurzelartmethode. Die nichtlinearen Vorgange 
werden auf 30 Seiten unter Verwendung der Beschreibungsfunktion und der 
Phasenebene zusammenfassend behandelt. Das Buch kann allen Theoretikern 
und Praktikern der Regelungstechnik sehr empfohlen werden. Ed. Gerecke 


Advances in Electronics. Bd. 5, von L. Marton (Academic Press Inc., 
New York 1953). XI + 420 S., 177 Fig., 1 Tafel; $9.50. 

Der fiinfte Band der von L. Marton herausgegebenen Buchreihe ‘“‘ Advances 
in Electronics’’ umfasst die folgenden Beitrage: 

R.C. Jones, Performance of Detectors for Visible and Infrared Radiation; 
R. W. Haywarp, Beta-Ray Spectvometers; F. E. Wittiams, Solid-State Lumines- 
cence; W. E. DANFORTH, Thorium Oxide and Electronics; H. C. WEINGARTNER 
and S. W. KENNEDY, A Review of Modern Vacuum Pumps in Electronics Manu- 
facturing; R. Q. Twiss, On the Steady-State Theory of the Magnetron; CH. J. 
Hirscu, A Review of Recent Work in Color Television; J. S. SCHAFFNER, Junction 
Transistor Applications. 


Aus dem dem Band 5 beigegebenen Gesamtregister fiir die bisher erschie- 
nenen fiinf Bande der Reihe geht hervor, dass in diesen alle Gebiete der Elek- - 


tronik, sowohl die grundlegenden physikalischen Erscheinungen als auch deren 
technische Anwendungen, in gleichem Masse zur Darstellung gekommen sind. 
Um dies in Zukunft noch mehr zu betonen, werden die folgenden Bande den Titel 
«Advances in Electronics and Electron Physics”’ tragen. Es ist dem Werk auch 


unter dem neuen Namen eine weite Verbreitung zu wiinschen. A.A. Rusterholz. 


_ Résistance des Matériaux théorique et expérimentale. Par RoBERT 
L’Hermite. Tome I (Dunod, Paris 1954). 860 p., 417 fig.; relié toile, fr.fr. 8400.—. 
_ L’ouvrage en quatre volumes que M. Rosert L’HERmiITE, directeur des 
Laboratoires du batiment et des travaux publics, a Paris, est en train d’écrire sur 
la Résistance des matériaux, est destiné aux ingénieurs et aux physiciens qui dé- 
sirent approfondir et élargir leurs connaissances en cette matiére. L’auteur 
cherche en quelque sorte a rédiger une encyclopédie de cette partie de la science, 
en dressant un tableau aussi clair et aussi complet que possible des lois et des 
expériences acquises dans ce domaine. II] y réussira certainement, si l’on en juge 
par le tome premier, remarquablement écrit, qui vient de paraitre. 

Le chapitre I de ce volume se rapporte aux lois et aux équations générales de la 
théorie de l’élasticité. Les chapitres II a IX traitent l’élasticité plane, la flexion 
et la torsion des corps a fibre moyenne, la théorie des structures isostatiques, 
celle des structures hyperstatiques, les plaques et dalles fléchies et les coques. 
Un dernier chapitre (X) est consacré a des critiques et 4 des remarques générales. 
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Ce volume est donc une ceuvre théorique. Il est destiné 4 donner une vue 
d’ensemble sur la théorie de l’élasticité. Il constitue également un recueil des” 
différentes méthodes de calcul et des cas d’intégration connus, mais un recueil” 
résumé, comme le dit l’auteur dans la préface. Il a été.en effet obligé, pour éviter 
une extension exagérée du texte, de condenser au maximum les exposés des diffé-— 
tentes théories. Il se borne méme, dans la plupart des cas, a en énoncer les prin-— 
cipes et les résultats essentiels. Toutefois, il a eu l’heureuse idée de placer une. 
bibliographie en cours de texte, et une autre en fin de chapitre, pour permettr q 
au lecteur de puiser facilement aux sources, s’il le désire. i. 

Une autre qualité de ce livre est que chaque chapitre et méme de nombreux 
paragraphes peuvent étre lus indépendamment, qualité que de nombreux lecteurs~ 
apprécieront. a 4 

Il est prévu que le second volume décrive les théories de la stabilité, du flam-_ 
bement et de la plasticité; le troisiéme traitera de la mécanique physique des” 
matériaux; enfin, l’étude des phénoménes dynamiques, chocs et vibrations, des 
associations et assemblages de matériaux, des statistiques et de la sécurité, 
formera la matiére du quatriéme volume. = 
Nous recommandons vivement l’ouvrage de M. L’HERMITE non seulement 
aux ingénieurs et aux physiciens de l'industrie, mais 4 toute personne se consa-— 
crant a des recherches dans le domaine de la résistance des matériaux, car cet 
ouvrage, facile et agréable 4 consulter, est une source inépuisable de renseigne-_ 
ments. H. Favre 


Sa Ne 


Limit Distributions for Sums of Independent Random variables. Von 
B. V. GNEDENKO und A. N. Kormocoroy, aus dem Russischen iibersetzt und 
_mit Anmerkungen versehen von K. L, Cuune (Addison-Wesley Publishing Com-— 
pany, Cambridge, Mass., 1954). 261 Seiten; $8.50. _ 
Jeder Mathematiker, der sich mit modernen Grenzwertsatzen der Wahr-_ 
_ scheinlichkeitsrechnung beschaftigt, wird sich mit dem obigen, ins Englische 
iibersetzten Buch der bedeutenden russischen Mathematiker GNEDENKO und 
KoLtmoGorov befassen miissen. Denn es enthalt eine klare, axiomatisch saubere 
Darstellung klassischer und moderner Grenzwertsatze der Wahrscheinlichkeits- 
rechnung und der damit zusammenhdngenden Begriffe, wie chardiccesigaeua 
Funktionen, unendlich-teilbare Verteilungen usw. Diese Grenzwertsatze bilden 
zum Teil die theoretische Basis der modernen mathematischen Statistik, es sei_ 
nur an das klassische Gesetz der grossen Zahlen erinnert. a 
__ Das vorliegende Buch setzt die einfacheren Erkenntnisse der Wahrschein- 
lichkeitsrechnung als bekannt voraus. Es gibt vor allem eine Ubersicht der fun-_ 
damentalen russischen Beitrage an die Wahrscheinlichkeitsrechnung im Gebiete 
der Grenzwertsatze und kann als wohl bedeutsamste Publikation iiber diesen 
Gegenstand in den letzten Jahren bezeichnet werden. = W. Saxer 


i Theorie ideal plastischer Kérper. Von W. PRAGER und P.G. Hopes, jra 
ins Deutsche iibertragen von F. CHMELKA (Springer-Verlag, Wien 1954). 274 S. 
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G. DOETSCH 
Flandbuch der Laplace-Transformation 
Band IIT 


Anwendungen der Laplace-Transformation 
2. Abteilung 
300 Seiten mit 44 Figuren 
Ganzleinen Fr./DM 40.-, broschiert Fr./DM 36.— 


Aus dem Inhalt: Partielle Differentialgleichungen — Differenzengleichungen — Integralgleichungen 

und Integralrelationen — Ganze Funktionen vom Exponentialtypus — Nachtrage zu Band I — Lite- 

farische und historische Nachweise — Biicher iiber die Laplace-Transformation — Literaturverzeich- 
nis — Sachregister — Berichtigungen zu Band II. 


Das «Handbuch der Laplace-Transformation» ist mit dem III. Band, der die zweite Abteilung 
der «Anwendungen» bringt, abgeschlossen. Der III. Band behandelt in der Hauptsache die Funk- 
tionalgleichungen (ausser den bereits im II. Band dargestellten gewohnlichen Differentialgleichun- 
gen), also partielle Differentialgleichungen, Differenzengleichungen und Integralgleichungen. Diese 
Gebiete sind sowohl fiir den reinen Mathematiker als auch fiir den Physiker und Ingenieur von 
Bedeutung. Insbesondere die bisher noch nirgends in Buchform dargestellte Behandlung der Dif- 
ferenzengleichungen mit Laplace-Transformation wird wegen des in letzter Zeit in der Technik 
neuerwachten Interesses an den Differenzengleichungen (Theorie der Kettenleiter, Schrittregler 
usw.) begriisst werden. 

Der Schlussteil beschaftigt sich mit der Theorie der ganzen Funktionen vom Exponentialtypus, 
die ein unerschépfliches Anwendungsgebiet der Laplace-Transformation abgeben. Hier sind noch 
viele Resultate zu erwarten, und die vorliegende Darstellung méchte zu weiterem Eindringen in 
dieses Gebiet anregen. 

Das Literaturverzeichnis enthalt die im II. und III. Band zitierten Arbeiten sowie die seit 1950 
erschienene Literatur aus dem Stoffgebiet des I. Bandes. 
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